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CHAPTER I 
DfERODUCTIOH 



1.1 Object and Scope 

This thesis is a study of the effectiveness of the surrovindlng 
soil in restraining the lateral movement of shedlow-buried arches subjected 
to blast overpressures from a nuclear explosion. An understanding of this 
effect is Important in the safe and economical proportioning of such 
structures. 

The theory is developed for parabolic arches with circular or 
barrel arches included as a special ceise. It can be applied equally veil 
to arches oriented with their longitudinal axes perpendiculeir to or parallel 
to the direction of travel of the blast wave. 

The arch is considered to be buried near the surface of a semi- 
infinite soil mass. The term "shallow-buried** as used herein refers to a 
configuration and amoxmt of earth cover such that reflected pressures and 
dynamic pressures need not be considered in the loading of the arch and yet 
not so great as to provide significant free-field attenuation of air-blast- 
induced overpressures in the soil. This is a convenient and meaningful area 
of investigation. It is convenient because it eliminates the dynamic and 
reflected pressure loading components at one extreme and the free-field 
attenuation of overpressure in soil at the other. Consideration of either 
of these extremes would serve further to complicate the problem \diile con- 
tributing little to the significance of the results. It is meaningful in 
that an economiceO. balance may be reached in prewitlcal installations by 
locating structures near the ground surface in order to minimize excavation 
requirements and facilitate access ^ile at the same time using earth cover 
to eliminate dynamic and reflected pressure loading components from the blast 
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loading yAilch the structviral system must withstand. In this study the depth 
of cover over the crown is limited to one -half of the eurch span. Ibis 
restricts the location of the structure to a region adjacent to the earth's 
surface in ^rtiich there appears to be no measurable attenuation of free-field 
overpressure in the soil and in which the restraint offered to soil movement 
by the overpressiire acting on the ground surface is significant. 

The two -hinged rather than the fixed arch is selected as repre- 
senting more closely the shallow-buried arch structures used in full-scale 
nuclear tests to date.^^^^^^^ While these structures were not actual ly pinned 
at the arch abutments, their behavior \onder blast loading appears to the 
writer to resemble more closely the pinned them the fixed case. Further, such 
techniques as separating the arch footings from the floor slab in order to 
reduce accelerations imposed on the floor slab, and keeping the arch footings 
small to encourage downward movement under blast loading in order to reduce 
direct compression stresses in the arch, tend to produce a more nearly pinned 
than fixed condition at the abutments. The greater structural flexibility 
of the two-hinged airch is an additional favorable feature in enabling under- 
ground structures to adjust to Imposed loading systems. 

The principed variables studied in the numerical examples include: 

1. Depth of cover 

2. Overpressure loading 

3 . Arch material 

4. Arch speua 

1.2 Method of Approach 

The arch is replaced by a model consisting of four inextensible 

(4) 

bars connected at moment-resistant joints. This model is developed 



(1) Numbers in parenthesis refer to the List of References. 
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herein in Section 5«1* Bie arch is thus reduced to a two-degree-of-freedom 
system. Tills greatly facilitates the dynamic analysis while at the same 
time permitting meaningful assumptions to he made regarding the loading and 
resistance parameters at the two haunches and at the crown. 

The two Independent variables chosen to describe the system are 
the horizontal displacement components at the haunches. While any two in- 
dependent displacements could have been chosen, these are most convenient 
since they provide a direct plctiure of the prlmeiry mode of response of the 
model and are used directly in the determination of the state of stress within 
the soil mass at each side of the arch. 

The two Independent modes of response of the model are referred to 
herein as the symmetric mode, in •vdilch the haunches move simultaneously in 
opposite directions, l.e. both inward or both outward, and the asymmetric 
mode, in \diich the haunches move simultaneously in the same direction, l.e. 
both toward the rl^t or both toward the left. 

Solutions for the response of this system are obtained using the 
Newmark p-methods. These are step-by-step mmierlcal integration 

techniques in which the accelerations, velocities and displacements of the 
dynamic system are determined at each time interval from the time of arrival 
of the shock wave until the significant response of the system has occurred. 
These techniques are readily adaptable for use in digital computers. One of 
the University of Ill inois Structural Research Laboratory Codes was modified 
and used for the solution of the numerical examples presented in this thesis. 

The blast overpressure loading is concentrated as discussed in 
Section 2.k and Chapter VI and is placed directly on the arch. Modification 
of the blast loading due to its transmission Ihro ugh the soil mass from the 
ground sxirface to the arch rib is handled separately as "soil resistance". 

This is discussed further in the following sections. 
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In general, the resistance of the arch-soil "strxicture" is handled 
by: 

1. Considering the bending resistance of the arch as it deforms, 

2. Considering the frictional and shearing stresses developed 
within the soil as the arch moves into or away from it, 

5. Combining these arch and soil resisteinces vectorially as they 
act upon the model, and 

4. Considering the mass effects of the arch and soil. 

While these resistances are closely interrelated and act in combination during 
the response of the system, it is convenient and admissible to handle them 
separately in developing the equations of motion for the system. 

The soil resistance is developed by considering the states of stress 
•vhich exist within the soil surrovinding the arch as the arch-soil "structure" 
deforms laterally. 

Direct compression in the arch is handled separately. This is an 
approximation but a necessary one with this model, since to Include it direct- 
ly would violate one of the assumptions upon which the model is based, i.e. 
the bars are inextensible. Note that it would also change the model from 
a two-degree-of -freedom system to a six-degree-of-freedom system. However, 
the resjjonse studies presented remain vedid with respect to lateral movements 
as long as the arch section is sufficient to withstand direct compression 
and buckling. A sepeirate "rib- shortening correction" (cf. Section 5* 5) is 
applied to certain of the resxilts to give an approximate picttire of its 
significance. This is discussed further in Chapter IX. 

The arch supports are considered to be unyielding. While this is 
a common assvraiption in structural problems, it is an approximation. Again, 
the assumption is necessary in order to limit the system to two degrees of 



freedom; however, this does offer an Interesting possibility for further 
study as stated in Chapter IX. 
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Throughout this thesis displacements, velocities and accelerations 
to the rig^t and upward are considered positive, and conversely. Loads are 
considered ixjsitive \dien they act on the structure to the right or upward. 
Resistances are considered as negative l0€ids. 

Ihe method of approach described above is approximate. The neces- 
sity of approximating actual conditions is discussed in Chapter II eind is 
also pointed out as occasions arise in the development of the loading and 
resistance expressions for the system. Ihe approach used gives calculated 
meignltudes of response which differ from values measvured in field tests. 
However, as will be seen, the results do describe the behavior of the actual 
structure and permit meaningful comparisons to be made of arch behavior 
under varying combinations of structure geometry and loading conditions 
(cf. Chapter IX ). 

1.3 Notation 

Each symbol used in the text is fully explained vhen it is first 
introduced. A summary of the most Important notation is presented here for 
the convenience of the reader. 

vertlcsil projection of a lower bar of the dynamic model 
(cf. Figure 2.l) 

cross-sectional area of the arch rib 

Influence coefficients for total joint rotation (cf. Equa- 
tion (3.2.11)) 

horizontal projection of a lower bar of the dynamic model 
(cf. Figure 2. l) 

span of the erch (cf. Figure 2.l) 

vertical projection of an iipper bar of the dynamic model 
(cf. Figure 2.l) 



a = 

A 

b 

B 

c = 
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c 

s 


seismic velocity of the soil 


C 

s 


apparent vertical component of the velocity of the shock 
wave in the soil 


d 


horizontal projection of an upper bar of the dynamic model 
(cf. Figure 2.l) 


d 

c 


depth of soil cover over the crown of the arch 


11 


vertical dimensions of the soil prisms used in the mass 
calculations (cf. Figures 5*1 and. 5*2) 


d 

o 


total depth of soil cover from the ground surface to the 
arch abutments 


E 


modulus of elasticity of the arch material 


J^ra 


equivalent mass factor (cf. Eqviations (7-l.^)) 


g 


acceleration of gravity, 52.2 ft/sec^ 


G . = 


equivalent meiss factor (cf. Equations (7.1.4)) 


h 

a 


thickness of the eirch rib 


II 

*o 


displacement 


11 


velocity 


• • 

“j - 


acceleration 


“j = 


effective external load 


H. 

J 


equivalent external load 


I 


effective moment of inertia of the arch rib cross section 


J 


as a subscript, unless otherwise specified, refers to the 
horizontal component at joint j 


J+5 


as a subscript, unless otherwise specified, refers to the 
vertical component at joint j 




influence coefficients for the equiveilent load on the dynamic 
model (cf. Eqxiations (6.5*2)) 


11 

"4*4 


Influence coefficients for equivalent soil stiffness 
(cf. Equation (4.4.10)) 


k 

o 


observed vedue of the ratio between the horizontal and 
vertical soil pressiires for the condition of earth pressure 
at rest 


11 


hydrostatic pressxire ratio for the "active" state of stress 
in the soil 
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Sa 






bending stiffness of the arch rib as concentrated at a 
joint of the dynamic model 

ratio of the soil pressure in the direction of P! to the 
vertical soil pressure for the condition of earth pressure 
at rest 

ratio of the soil pressure in the direction of P' to the 
vertices, soil pressure for the condition of earth pressvure 
at rest 

hydrostatic pressure ratio for the "passive” state of stress 
in the soil 




ratio of the horizontal to the vertical overpressure induced 
in the soil by the shock wave in air 




soil stiffness dxjring the development of an "eu:tive" state 
of stress in the soil adjacent to haunch j 




= soil stiffness during the development of a "passive" state 
of stress in the soil adjacent to haimch j 



L 

m 



j 

j 




m 



6j 



m 









M 



length of b6U j of the dynamic model (cf. Figure 2.l) 

concentrated mass at a joint of the dynamic model 

arch contribution to m . 

J 

soil contribution to m . 

J 

equivalent mass 

moment resistance of the arch-rib cross section 

maximum moment -res is ting capacity of the arch rib 
cross section 



n = as a subscript, refers to time 

n+1 = as a subscript, refers to time t = 

p = overpressure in the shock wave in air 



P 



so 



peak value of the overpressure in the shock wave in air 



P 



a 



vertical component of the average overpressure acting on 
the arch 



= uniform radial component of the average overpressure acting 
on the arch 



= average overpressure acting at the groxind surface on 
loading region j 
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^Sj 


average oveirressiire acting at the grxjund svirface on the 
potential "active" failure vedge in the soil adjacent to 
haunch j 




average overpressure acting at the ground surface on the 
potential "passive" failure wedge in the soil adjacent to 
havinch j 


P' 


force on the side of the arch due to the weight of the soil 


p" 


force on the side of the arch due to the "svu*charge" of the 
blast oveirpressvire acting on the ground surface 




force acting on the side of the arch due to the weight of 
the soil and corresponding to a fully developed "active" 
state of stress in the soil mass 


'■I " 


force acting on the side of the arch due to the blast over- 
pressure and corresponding to a fully developed "active" 
state of stress in the soil mass 


^OA 


component in the direction of P' of the force acting on the 
side of the arch due to the weight of the soil and correspond 
ing to the condition of earth pressure at rest 


P” = 

OA 


component in the direction of P^ of the force acting on the 
side of the arch due to the blast overpressure and corre- 
sponding to the condition of earth pressxire at rest 


P' = 

OP 


component in the direction of Pp of the force acting on 
the side of the arch due to the wel^t of the soil and 
corresponding to the condition of earth pressure at rest 


"op = 


component in the direction of Pp of the force acting on 
the side of the arch due to the blast overpressure and 
corresponding to the condition of earth pressure at rest 


"p ■ 


force acting on the side of the arch due to the weight of 
the soil and corresponding to a fully developed "passive" 
state of stress in the soil mass 


p" = 

p 


force acting on the side of the eirch due to the blast over- 
pressure and corresponding to a fully developod "passive" 
state of stress in the soil mass 




influence coefficients for equivalent arch resistance 
(cf. Equation (5.4.?^ 


= 


equivalent resistance due to deformation of the eirch rib 




effective resistance at an "interior" haunch due to deforma- 
tion in the soil mass 




effective resistance at an "exterior" haunch due to deforma- 
tion in the soil mass 
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equivalent resistance due to deformation in the soil mass 
radius of a circular eirch 

maximum value of soil resistance (based upon a consideration 
of gravity forces) for a fully developed "active” state of 
stress in the soil 

maximum value of soil resistance (based upon a consideration 
of overpressure forces) for a fully developed "active” state 
of stress in the soil 

maximum veilue of soil resistance (based upon a consideration 
of gravity forces) for a fully developed "passive" state of 
stress in the soil. 

maximum value of soil resistance (based upon a consideration 
of overpressure forces) for a fully developed "passive" state 
of stress in the soil 

shear in beir J 

time 

elapsed time following passage of the shock front at a point 
of interest 

duration of the positive phase of the shock wave 

time of arrival of the shock front at the windward boundary 
of loading region j 

time of arrival of the shock front at the leeward boundary 
of loading region j 

thrust in the arch rib 

natural period of vibration of the system 

velocity of the shock wave in air 

effective inertial force 

equivalent inertial force 

element of the constraint loadings 

one -half the angle of opening of a circular arch 

weight per vinit volume of the arch materisJ. 

wei^t per unit volume of the soil 
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At 



friction angle corresponding to the friction developed 
between the surface of the arch and the soil 

"rib -shortening" correction 

change in length of bar J due to thrust in the eirch rib 
rotation of beir j 

time interval used for numerical integration 



T) = inclination with respect to the horizontal of the shock 

frt)nt in soil 



\ 

m 









V 










T 



J 






inclination of bar J with respect to the horizontaO. 

multiplying factor for the constraint loadings 

influence coefficient for equivalent mass (cf. E]uation(5.2.6)) 

influence coefficient for equivalent mass (cf. Equation (%2.6)) 

rotation of the lower bar of the model at haunch / 

rotation of the lower bar of the model corresponding to full 
development of am "active" state of stress in the soil 
adjacent to haunch t 

rotation of the lower bar of the model corresponding to 
full development of a "paissive" state of stress in the 
soil aidjacent to haiuich t 

influence coefficient for eiuiveilent arch resistance 
(cf. Equation (7.2.5)) 

angle of internal friction of the soil 



. = rotation of joint j 

J 

ijr = joint rotation required to produce yielding of the joint 

y 

w = inclination of the lower bar of the model with respect to 

the vertical 



b) s 

m 



circular natural frequency of the system 
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CHAPTER II 

GENERAL DISCUSSION OF ^ VARIABLES 
2.1 Arch Resistance 

The inertial resistance of the arch is discussed separately in 
Section 2.3* The discussion here considers the arch resistance in terms of 
the forces required to pull a raffssless, deformable structure into the deflec- 
tion configuration defined by the two independent displacements. 

For the purpose of the dynamic analysis the arch, which has an 

infinite number of degrees of freedom as it responds to dynamic loading, is 

replaced by a mathematical model. This model consists of four rigid, inex- 

tensible bars connected together at three moment-resistant joints and 

connected to the eirch abutments at two frictionless hinges as shown in 

Figure 2.1. Ihls model has two degrees of freedom. Note that it permits 

assumptions to be made regarding conditions at the two haunches and at the 

( 8 ) 

crown. It has been demonstrated' ' that these are the critical sections of 
arch behavior, and thus the model provides a relatively simple though useful 
tool for the study of the dynamic response of the arch Itself. 

In this thesis, an arch of unit width is considered. Also it is 
assumed that the arch section under consideration is sufficiently far removed 
from endwalls or intermediate stiffeners that they have no effect on its 
dynamic response. This, in effect, reduces the problem of one of two 
dimensions. 

Hie flexural resistance of the erch is determined herein by first 
finding the moment-resisting capacity of the arch cross section and then 
concentrating this moment resistance at the three moment-reslsteint Joints 
of the model. Hie joint resistance is then represented by an idealized 
moment-rotation relationship. Ihls is believed to be a reasonable 
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representation of the flexural resistance of the arch, since it will be shown 
(cf. Chapter IX ) that the flexural resistance of the arch contributes rela- 
tively little of the total resistance of the soil-arch ''structxxre” to lateral 
deformation vinder blast loading. Thus, modification of the flexural resistance 
to allow for dead, and live loads acting on the arch prior to the anrival of 
the shock wave or for thrust in the arch due to the blast loading does not 
appeair warranted and no modification is made herein. 

Since the model used to represent the arch is comixjsed of inex- 
tensible bsirs, it is necessary to handle the direct compression resistance 
of the arch separately. This is an approximation; however, its effect upon 
the lateral deformation response of the soil-arch "structure" is not considered 
serious as long as the 8U*ch section is sufficient to resist the maximvan thrust 
imposed on the arch by the blast overpressure without general yielding or 
buckling inward. Well compacted soil surrovinding the arch shoiild restrain 
any tendency of the arch to buckle outward into the soil mass. It is important 
to note at this point also that the arch section must be sufficient to with- 
stand imposed dead and live loadings during construction and by normal service 
conditions prior to the nuclear explosion. Also worthy of note is the probable 
loading condition on the eu*ch after the nuclear blast wave heis passed and 
dynamic response has ceased. The importeince of attention to these conditions 
in the design of any buried structure is obvious. 

An approximate "rib “shortening correction" is applied to certain 
of the nvimeriC6il results obtained in this study. While this is not directly 
a part of a study of lateral deformation, it does give some indication of the 
total deformation which may be expected under the given conditions. It also 
serves to indicate the effect of direct compression in the arch rib on arch 



response. 
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2.2 Soil Resistance 

The present state of knowledge regarding the interaction of soil 
and underground structures contains many uncertainties. This is true for 
"static" behavior and is even more true for "dynamic" behavior such as that 
now under consideration. The development in this thesis is considered 
reasonable but is necesseirily approximate, conseiuently a hi^ degree of 
refinement is not attempted. 

As with the arch discussion in the preceding section, the inertieJ. 
effects of the soil eire covered separately in Section 2.5. The discussion 
here considers the shear and frictional forces developed within the soil 
surrounding the arch as the soil-arch "structure" deforms under blast loading. 
It is assumed throughout that prior to the arrival of the blast wave the soil 
surrounding the arch is in the "at rest" state described below. 

The concept of soil resistance used in this thesis derives from a 
consideration of the states of stress which eure developed within the soil 
mass 8\irrounding the eu:ch as the eurch xindergoes lateral deformation. Ihe 
confining action of the blast overpressure on the soil is also considered. 

The states of equilibrium within the soil mass which have been 

(9) 

used in developing this concept are: 

1. Active State of Plastic Equilibrium - the state of stress 
within a soil meuss which has undergone a "stretching," i.e. the soil mass 
has been permitted to expand laterally, \intil the ratio of the horizonteJ. 
to the vertical press\ire within the soil mass has reached a minimum value, 
called the "coefficient of active earth pressxire." 

2. Passive State of Plastic Equilibrium - the state of stress 



within a soil mass which has undergone a lateral compression until the ratio 
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of horizontal to vertical pressure has reached a maximum vaJLue, called the 
"coefficient of passive earth pressure." 

3 . Earth Pressure at Rest - this is a state of stress which is 
intermediate between the two states of plastic equilibrium described above 
for the soil mass and exists after the soil has been deposited by natural 
and artificial processes in the vicinity of the eurch. It is a state of 
elastic equilibrium. Since this definition includes "artificial" processes 
for deposit of soil, it will be restricted herein to require careful com- 
paction for optimum density of all fill materials used during construction 
and backfilling of the arch. 

The hydrostatic pressure ratios, or ratios of horizontal to 
vertical soil pressure, referred to in the foregoing definitions differ 
because of the differences in magnitude and direction of the shear and 
frictional forces in the soil mass. In general, these stresses are developed 
within the surrounding soil as the arch-soil "structiare" deforms under blast 
loading. The extent of development of these stresses at any particular time 
depends ui)on: 

1. The state of stress within the soil mass, and 

2. The restraint imposed upon the deformation of the soil mass 
by the confining action of the blast overpressure. 

Since the blast loading is applied to the eirch rib throu^ the soil, these 
stresses will assist the arch in "resisting" the loading. 

The degree of restraint imposed upon the deformation of the soil 
mass by the blast overpressure will, of course, vary with time as the blast 
overpressure varies. 

Hie primary variables determining the state of stress within the 



soil mass are the direction and amount of deformation to which the soil mass 



I 
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has been subjected. These will, of course, vary with time, but, at any 
given time, will also vary from one region to another around the arch. In 
this thesis, three regions are considered, viz. the regions adjacent to each 
of the haunches and the region adjacent to the crown. 

As noted above, prior to any deformation it is assumed th- t the 
"at rest" stress state exists in ihe soil surrounding the eirch. As a haunch 
moves inward, an "active state of plastic equilibrium" tends to develop. 

This does not develop fully immediately but requires a certain amount of 
deformation within the soil mass as discussed below. As a haunch moves out- 
ward into the soil mass, a "passive state of plastic equilibrium" tends to 
develop. This, too, does not develop fully until sufficient deformation 
has occurred. 

Ihe eunount of deformation re-iuired for the full develofment of an 
active state of plastic equilibrium within the soil mass was studied by 
Dr. Keirl Terzaghi. In a series of tests on a large-scale retaining wall 

Dr. Terzaghi demonstrated that the lateral force on a wall varies as the 
wall undergoes lateral movement. The relationship between the force and 
the movement is shown qualitatively in Figure 2.2.^^^^ The resistance- 
displacement fimctions developed in this thesis are essentially idealizations 
of this figure. 

While Figure 2.2 is qualitative only, Dr. Terzaghi and others 
have estimated the amount of wall movement re uired to reduce the eeu*th 
pressxire to the active state for retaining walls and specific soil conditions. 
These recommendations form the basis for the quantitative assumptions of 
deformation required to develop fully the active state of stress which are 
used in this thesis. There is little guidance in tlie literature regarding 
the amount of deformation required to develop fully the passive st? te of 
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stress within the soil mass. Limited investigation^^*^^ ^ heis indicated, 
however, that it is greater than that required to develop the active 
state of stress. Accordingly, in this study it is assumed to be such that 
the ideeilized resistance -displacement function for the soil has a constant 
slope in the active and passive "elastic" regions for a given value of over- 
pressure restraint (cf. Section 4.4). 

In determining numerical values of the hydrostatic pressure ratios 

corresponding to the limiting states of stress described above the classiceil 

(14) 

Co\ilomb assumptions have been used. It is assumed that: 

1. The soil slides along a solid pressure surface, i.e. the soil 
failure wedge is bounded on one side by the pressure sxirface and on the 
other side by the surface of rupture. 

2. The obliquity of the resultant pressure on the pressure surface 
has a vedue less than or equal to the friction angle of the soil. 

5. Sufficient movement has occurred within the soil mass to 
develop the state of stress being considered, i.e. active or passive. These 
assumptions yield directly the pressure acting on the pressure surface. A 
corapeirison between results obtained from the Coulomb Equation and those 
obtained by the Trial Wedge Method is given in Api>endix A. The equation 
gives results which are considered sufficiently accurate for this study. 

An additional assumption re ^uired to define the soil failiire wedge 
is the character of the surface of rupture within +he soil mass. A plane 
surface of rupture is assumed throughout. 

The pressure surface is taken as the lower bar of the model on 
each side and is considered to continue upward in the same direction from 
the haunch to the ground surface. This facilitates the development of tlie 
soil resistance concept used in this study and, for the relatively shallow 
depths of burial considered, does not appear unreasonable. 
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The soil resisteince at the crovn is limited primarily to the 
inertial resistance of the soil at and adjacent to the crown. While there 
must obviously be shearing forces developed in this soil mass as the arch- 
soil "structure" deforms, the relatively small volume of this region for 
the shaJ-low depths of buried considered, eis compared with the volumes of 
the soil wedges influenced by lateral movements of the haunches, indicates 
that these shearing forces will be relatively small. Also, it is possible 
that they will act, at least in part, only to carry the load to another 
portion of the arch. The extension upward to the ground surface of the 
pressiire surfaces at the sides of the arch as described in the preceding 
paragraph is an attempt to include that portion of these shearing forces 
which will be effective in reducing the load brought to the arch. No load 
reduction is made, i.e. no "deformation resistance" is developed, for the 
remainder of the soil in this region. While this approach is somewhat 
arbitrary, it is not considered to have a significant influence on the 
response of the system as a \rtiole. 

Mass Effects 

The inertied resistance of the arch-soil "structure" has ein 
important effect on both the magnitude and time history of its dynamic 
response. As noted previously in Section 2 . 1 , the model used in this study 
for the dynamic analysis permits appreciable latitude in the choice of the 
masses which are considered to move with the system at any given time. This 
property of the model will be used to advantage in adapting the model to fit 
changing conditions dviring the response history of the actual structure. 



The mass of the arch is concentrated at the joints connecting the 
inextensible, massless bars, except for the mass of the lower i>ortion of the 
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arch on each side. This portion of the mass of liie arch is concentrated 
at the abutments, i.e. is not considered to have a significant effect upon 
the response of the system. This technique is commonly used in dynamic 
analyses and gives reasonable results. 

The mass of soil \diich is assinned to move with the eirch is 
selected according to: 

1. The location of the joint, 

2. The direction of movement of the joint, and 

5. The depth of burial of the eirch. 

Varying the mass assumption with the direction of movement of the 
joint has the effect of giving the dynamic system a variable mass. While 
this increases the complexity of the eunalysis somewhat, it may be handled 
readily by the nvunerical integration procedure, and it appears to give a 
better approximation to the actual field condition than the use of a constant 
mass assumption. Logic indicates, for example, that the quantity of soil 
influenced by the inward movement of a haunch is different from the quantity 
influenced by the outward movement of the haunch into the soil. 

Limits are placed on the dimensions of the soil masses considered. 
In no case is a dimension considered which is greater than the corresponding 
speui of the arch.^^^^ 

2.h Blast Pressure Loading 

Since the object of this study is to investigate the response of 
a shallow-buried arch to blast loading from a nucleeir explosion, the blast 
overpressure is the primary weapon effect considered. In applying the 
results of this study to any field situation, other weapon effects obviously 
must be considered as well. Further, the fact that the limiting cases of 
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the numerical examples studied herein would not prove feasible in the field 
because of other-than-blast effects shoxild not detract from their value in 
developing an understanding of structural response to the blast overpressure 
itself. 

The results of this study are considered applicable to stnictures 

located outside of the crater, i.e. outside both the zone of rupture and 

^ l6) 

the zone of plastic deformation. ' ' It is doubtful that structures can be 

built to withstand weapons effects levels at closer ranges. 

Direct ground shock is not considered. For structures located near 

(l‘ 

the gix)\ind surface and outside the crater this omission appears reasonable. 
Much of the energy of the direct ground shock is dissipated in disruptive 
effects in and immediately adjacent to the crater. Fvurther, as the direct 
ground shock wave travels outward from the point of the explosion, its rise 
time increases so that it becomes even less effective in causing dameige. 

Thus, the air -induced ground shock is considered the major stinic- 
tural damaging effect. It is assvnned that if the structure can withstand 
the air-induced ground shock it can resist both the direct ground shock and 
differential soil movements. This latter consideration may become more 
imixjrtant with larger structures than those considered herein. However, no 
truly "large” underground arch structures have been tested to date, and the 
possible importance of differential soil movements in this case has not 
been established. ^ 

The shock wave in air is assumed to have a single front, i.e. the 
structure is assumed to be in the region of Mach reflection. This assvmiption 
is not necessary but is convenient when defining the load on the structure. 
The study may be extended readily to cover other types of blast waves if 



desired. 
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No attenuation of free-field overpreBsures with depth is assumed. 
For the relatively shallow depths considered in this study, this is entirely 
resLSonahle (cf. Section l.l).^^®^ 

The decay of fr«e-fleld overpressure with time in the soil mass is 

assumed to be the same as that of the shock wave in air. Again, for the 

relatively shallow depths considered in this study, this is a reasonable 

(18) 

assumption. 

The free-field horizontal -to -vertical overpressure ratio in the 

soil is assumed as l/4. This is the lowest of the values normally recom- 
(l5 17 19) 

mended ’ ^ y > gjjjj -^,^8 chosen to accentuate any tendency of the enrch to 
move latereilly. A low value of this ratio acts, of course, to decrease the 
horizontal component of the pressure acting on the arch during passage of the 
blast wave and thus to decrease the restraint vAiich the load itself offers 
to lateral movement. 

A plane shock front is assumed for the shock wave in soil. A study 
of the inclination of this shock front is presented in Appendix B. The rise 
time of the overpressure at this shock front is assumed instanteuieous. Ihe 
fact that studies and field measurements do indicate a finite and increasing 
rise time as the shock wave travels downward from the ground svirface' ' is 
partially allowed for in developing the loading on the arch by considering 
the average overpressure over a region of Interest rather than the 
instantaneous overpressure at any particular i>oint. 




Figure 2.1. Model Considered in Analysis 




Movement of Woll 



Figure 2.2. Effect of Movement of a Wall on 
the Lateral Tlirust^^^^ 
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CHAPTER III 
ARCH RESISTANCE 

In developing an expression for the resistance offered by the arch 
to lateral deformation under nuclear blast loading one must consider the manner 
in which the eLrch will respond and the relative importance of the veirious re- 
sistance peirameters. In general, as the ea:ch responds to a dynamic loading, 
the arch will mobilize both inertial and deformation resistance. This chapter 
is limited to a discussion of resistance due to the deformation of the arch 
rib. Inertial effects are discussed separately in Chapter V. 

As the arch moves laterally in response to a dynamic, unsymmetrical 
loading, significant deformations occ\ir at the haunches and at the crown. It 
has been shown' ^ that these are ailso the critical sections for bending moment 
and that the thrust is nearly uniform throughout the arch. 

3.1 Development of Model 

Since the arch Itself is continuous, it has an infinite number of 
degrees of freedom as it responds to a dynamic loading. Thus, it is incon- 
venient to aneilyze the eirch directly. In this study, for the purpose of the 

dynamic aneilysls, the arch is replaced by the four-beir model shown in 

(U) 

Figure 2.1. ' ' 

The model consists of four inextensible bars connected at three 
interior Joints. At the interior Joints are concentrated both the moment 
resistance of the eirch rib and the mass of the arch. These Joints also act eis 
node points. Ihe model permits the introduction of load and soil resistances 
into the analysis at each Joint location. As mentioned in Section 2.1, the 
four -bar model is connected to the arch abutments by frictionless hinges. 

This is permissible for the two-hinged arches studied in this thesis. 



23 



If the bars were capable of deformation, the model vould be a six- 
degree -of -freedom system. Since the bars Eire inextensible, one can write four 
"constraint" relationships among the joint displacements, one "constraint" 
relationship for each bar. This reduces the nvonber of degrees of freedom to 
two. These "constraint" relationships are based upon the condition that the 
deformations of the model at the joints correspond to no change in the length 
of the members. These relationships are developed for the model in Section 3-2. 

In this study the constraint relations will be derived using four 
loading patterns which produce only axial stresses in the bars. They thus 
cause no movement of the model and, of course, do no work when applied to the 
model. The principle of virtual work and the fact that these force pattenis 
cause no movement Eire then used in deriving relationships among the displace- 
ments smd among any general set of loading components. These relationships 
are then consistent with the two -degree -of -freedom system. 

A limitation in the use of this model is that direct compression in 
the arch cEinnot be handled directly. This follows from the inextensible nature 
of the bars. This limitation would be serious if this study were of the total 
stability of the arch, since the Eirch rib must be able to withstand the thrust 
loading Imposed upon it without failing by general yielding or buckling. 

However, if the arch section is sufficient to carry the imposed thrusts, it may 
then be studied from the standpoint of latereil deformation response without 
simultaneous consideration of direct compression load components. 

The prime reason for accepting this limitation in the model is, of 
course, to retain the convenience of the two-degree-of-freedora system. If 
the bars were extensible, the model would have not two but six degrees of 
freedom and the analysis would become correspondingly more involved. The nature 
of the behavioral assumptions upon which this study is based is approximate 
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enou^ that a more refined analysis is not considered j'ustified at the present 
time. 

The independent displacement components chosen to define the con- 
figuration of the model at any instant are the horizontal displacements of 
the haunches. While other pairs of independent displacements could be used, 
those chosen are particularly convenient in that they: (l) indicate directly 

the primary mode of response of the system (cf. Article 2.l) and (2) simplify 
the expression of the soil resistance (cf. Chapter IV ). 

5.2 Constraint Relationships 

The constraint relationships relate the various eirch parameters and 
serve to reduce the continuous arch to a two-degree-of-freedom system. These 
relationships are developed herein by considering four independent loading 
patterns, each of which produces only axial forces in the bars of the model. 
Thus, since the bars are inextensible, no deformation of the arch can take 
place under any of these loading patterns. 

The four loeiding patterns chosen are shown in Figure 5*1* For con- 
venience these loading patterns will be referred to herein as the "constraint 
loadings". Other loading patterns could be developed by combining these; 
however, these possess certain orthogonal properties which are convenient in 
deriving the following relationships. The constraint loadings may be tabxilated 
conveniently as shown in Table 5»1* 

5.2.1 Loading Relationships 

The load concentrations which will be applied to the model are 
shown in Figure 5.2. Relationships among these load components are now derived 
using the constraint loadings shown in Figure 5.I and the fact that these con- 
straint loadings, when applied to the model, produce no movement of the model. 
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Now any general system of loads can be replaced by a set of loading 
patterns consisting of some proportion of each of the four constraint loadings 
plus a loading pattern \rtiich does produce movement of the model. This latter 
pattern is that portion of the original system of loads which is effective in 
producing dynamic response of the system. This concept is expressed in 
Equation (3.2.1). 



H 



J 



k 



m= 1 



X 

m 




+ H, 



( 3 . 2 . 1 ) 



where Hj = einy system of loads on the model (cf. Figure 3.2) 

= multiplying factors for the constraint loadings 

a. = elements of the constraint loadings 

J,m 

H. = that portion of the loading system which is effective 
in producing motion of the model 

j = load component index; j = 1,2,... 6 

ra = constraint loading index; m = 1,2,3>^. 

In developing the loeid relationships from Equation (3.2.I) the 
constraint loadings serve to reduce the number of equations from six to two. 
This is accomplished algebraically by seeking four values of the multiplying 
factor, X^, such that four of the quantities Hj are zero. Since the horizontal 
displacements at the haunches have been chosen as the independent variables, 
the corresponding load components and will be retained here. Thus, in 

the following derivation multiplying factors are sought such that = 

2 = s 0. 

Using either Figure 3.I or Table 3«1 and Equation (3.2.I), one may 



write: 



H, 






- a X^ 



+ a Xi 






+ 2X^ 



0 

0 

0 



- a X^ 



(3.2.2) 



- a X|^ = 0 
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Note that the constraint loadings were chosen so that Equations 
( 5 . 2 . 2 ) are orthogonal, i.e. ^ = 0, for m ^ n and J = 2,i^,5,6. 

J 

Thus, the raviltiplying factors, X^, may be obtained from: 



X = 
m 



') H, 

J 



, J = 2,4,5,6 



(3.2.5) 



Equation ( 5 . 2 . 5 ) gives directly: 

\ 



^ i O*!. ^ ®6> 

c 



(3.2.4) 



S ~ 2 d ^® 2 ^ 

\ i (% - %) 

Substituting Equations (5.2.4) into Equation ( 5 . 2 .I) and solving for 
and yield the load relationships for the two -degree -of -freedom system. 
They express the "equivalent" load on the system, i.e. that portion of the 
total load which is effective in causing movement of the structure in the 
horizontal direction at the haxmches. 



“1 - * I (1 - - <X - |(f - 

% - 4 * -S - * <i>"6 

These load relationships will be vised later in developing the 
equations of motion for the model (cf. Chapter VIl). 



( 3 . 2 . 5 ) 



5 . 2.2 Deflection Relationships 

In deriving relationships among the displacement components of the 
model use is again made of the constraint loadings of Figure 5 .I together with 
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the principle of virtued work. Since none of the loading patterns produces 
any displacement of the model, the product of the loadings in the pattern 
times the displacements in any consistent set of displacements of the struc- 
ture in the direction of the loadings in the pattern must he zero. This is 
expressed algebraically in Equation (5-2.6). 



z 

J=1 



h, a = 0 



( 3 . 2 . 6 ) 



where 



h . = displacement components 

J 

ot. = elements of the constraint loadings 
J = displacement component index; J = l,r, ...6 
m = constraint loading index; m = 1, 

Again referring to the choice of the horizontal displacement com- 
ponents of the haunches as the Independent variables, the other displacement 
components may be written in terms of them using Equation (5.2.6). 

Solving Equation (5-2.6) explicitly for each of the deflection 
comjKjnents yields: 

= |(1 - bc/ad)(h3^ + h^) 
h|^ = - (b/a)h^ 



hj = |)(h, - hj) 

h^ = (b/a)h^ 

where a,b,c,d = generalized arch dimensions (cf. Figure 2.l) 
hj^ = horizontal displacement component of Joint 1 
h,^ = horizontal displacement component of joint 2 
h^ = horizontal displacement component of Joint 5 
hj^ = vertical displacement component of Joint 1 
h^ = vertical displacement component of Joint 2 
hg = vertical displacement comijonent of Joint 5 



(5.2.7) 
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5.2.3 Joint Rotation Relationships 

It is necessary to know the relationships between the rotations of 
the joints and the independent displacement variables in order to determine 
the moment resistance developed in the joints of the model for any particular 
deflected shape. Since the horizontal displacement components of the haunches 
have already been chosen as the Independent variables, they will be used as 
such in this derivation also. 

Referring to Figure 3*3^ it can be seen that for small displacements 
the rotations of the bars may be expressed as: 

(h . , " h .) CSC 9 , 






A 



( 3 . 2 . 8 ) 



where = rotation of bar j 



h., h. ^ = horizontal displacement components at the ends of bar j 
j J “ J- 



L. = length of bar j 
J 

6 . = inclination of bar j 
J 



Substituting into Equation ( 5 . 2. 8) and rearranging terms give the 
following expressions for the bar rotations of the model in terms of the in- 
dependent displacement coordinates and the geometry of the sirch: 

60 ^ Z - h^/ a 

b0» - + ^(l/c + b/ad)h^ - ^l/c - b/ad)h, 

2 ^ ^ ^ (3.2.9) 

60^ = - ~(l/c - b/ad)h^ + |(l/c + b/ad)h^ 

60 j^ r - h^/a 



The total rotation at each of the joints is simply the algebraic 
difference of the rotations of the bars connecting at the joint. This may be 
expressed as : 



^. = &0 
J 



J+1 




(3.2.10) 
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where = to ted rotation of joint j 

= rotation of bar j. 

Substituting Equations (5-2.9) into Equation (5.2.10) and intro- 
ducing simplified notation give the desired joint rotation relationships. 

^^1-V5 

= -Aghj^ + A^h^ (5-2.11) 

where A^ = ^( 2/a + l/c + b/ad ) 

= l/c 

A^ = |(l/c - b/ad) 



5.5 Moment -Rotation Relationships for Joints 

As the arch displaces laterally under load, bending moments are 
developed in the arch which act to resist the displacement. Likewise in the 
model, as a change occurs in the angle between the bars meeting at a joint, 
a resisting moment is developed in the joint. In Section 5-2.5 expressions 
are derived relating joint rotations and model displacements. It is the pur- 
pose of this section to derive relationships between these joint rotations and 
the corresponding moment resistance developed by the joints. 

If bending in the arch is assumed to be viniform over a region such 
as that shown in Figure 5-^^ "the moment -rotation relationship for the joint 
may be derived as follows: 

Let M. = moment resistance of the arch rib cross section 

u 

E. = modulus of elasticity of the material, from vdilch 
^ the arch rib is made 



I . = moment of inertia of the arch rib cross section 
J 

L. = length of bar (assimied of equal length on either 
side of joint) 



^ = joint rotation 

J 
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Then, from Figure 5.4: 





2 



Hence : 



M 



J 




For an arch of constant cross section, E and I. may be considered 

J J 

constant. may be chosen as a constant ^en setting up the model. For these 
conditions it is convenient to simplify the notation as: 



Mj = Kj (5.3,1) 

J 

where K_ = E I./l. = bending stiffness of the arch rib as concentrated 
j J J J at joint J of the model. 

Strictly speaking. Equation (3.5*l) is valid only in the elastic 
range. Since it is intended to permit plastic deformation of the arch rib, 
the moment -TO tat ion relationship is extended as shown in idealized form in 
Figure 5.5. Ihis follows directly from Equation (5 . 3 • l) and the requirement 
that the arch rib possess sufficient ductility to continue to rotate without 
rupture or decreasing moment resistance once the maximum moment resistance is 
developed. This can be achieved readily in the design of the arch rib. 

The limits of this moment-rotation diagram for any particular case 
will thus be determined by the strength, ductility and geometry of the arch 
rib cross section as well as by the geometry of the arch itself. As mentioned 
previously (cf. Section 2.l), no modification of bending resistance is made 
for thrust in the arch caused by the blast loading, or for dead and live 
loads existing on the arch prior to the arrival of the shock wave. In general, 
these effects will tend to reduce the stiffness of the arch. Again, while 
the arch must be capable of withstanding these loads, their net effect on the 
total resistance of the arch-soil "struct\ire"to lateral deformation is 
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relatively small and such modification is not considered warranted for this 
study. 



Development of Equivalent Arch Resistance 

In order that the moment resisteince of the Joints may be inserted 
conveniently into the equations of motion for the model, relationships are 
needed expressing the resistance of the arch to lateral deformation. These 
are derived in this section first in general terms end then specialized for 
the choice of independent displacement variables used in this thesis, viz. 
the horizontal displacements of the haiinches. 

In the sign convention used for this derivation clockwise shears are 
considered positive. Bending moments are considered positive when they tend 
to produce compression in the outer fibers of the arch rib. Bar rotations 
are measured at the left ends of bars, counterclockwise from a horizontal 
axis directed toward the right. As stated in Section 1.2, loads and displace- 
ments are taken as positive horizontaJJ.y to the rl^t and verticailly upward. 

Referring to Figure ^,6 and considering the equilibrium of either 
one of the bars, it can be seen that: 



S 






M. 



A 




( 5 .‘». 1 ) 



where Sj = shear in bar J 

M. = resisting moment in joint J 
J 

Lj = length of beir j 

Svumning horizontal forces acting on the free body shown in 
Figure 5*7 gives: 



( 3 . ‘♦. 2 ) 
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where 



Qj = horizontal component of arch resistance at joint J 
Sj = shear in bar j 

= inclination of bar j 

60j= rotation of bar j xmder nucleeir blast loading. 

Slimming vertical forces acting on the free body shown in Figure 5«7 



gives : 

where ^j+5 “ vertical component of the arch resistance of Joint j. 

In deriving the equivalent arch resistance at the haunches, it is 

convenient to express the inclination of the sheeiring forces in the bars in 

terms of the arch geometry. This is readily accomplished using the familiar 

trigonometric identity stated in Equation 

sln(0, + 60.) = sin 0. cos 60 + cos 0. sin 60. (5-^«*+) 

J J J J J J 



Using the notation of Figure 2.1, and limiting consideration to small rotations, 
one can write: 



sin(0^ + 60^) = ^ [a - (“)h^) 

^ [c + |(d/c + b/a)h^ - ^(d/c - b/a)h^] 
^ [-C - |(d/c - b/a)h^ + |(d/c + b/a)h^] 



sin(02 + 60g) 



( 5 . *+. 5 ) 



sin(0^ + 60^) 



sin(0j^ + 60j^) = ^ [-a - (b/a)h^J 



It is now possible to write expressions for the equivalent arch 
resistance in the horizontal direction at the haunches by substituting 
Equations (5.4.l) and (5.4.5) into Equation (j.4.2). Ohus: 



I 
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- + -)h- 
c a 1 



2^c 





- + ^)h, 

2 c a 1 2'c a^ 5 







0 . 4 . 6 ) 



vhere and Q_ = horizonted components of arch resistance at joints 1 

^ and j3 respectively 

M|, M , M, * resisting moments developed in joints 1, 2 and 3 
^ respectively 



a,b,c,d = generalized eirch dimensions as given in Figure 2.1 

h-. = horizontal components of displacement at joints 1 
^ and 3 respectively. 

By limiting consideration to equal bar lengths, rearranging terms 
and simplifying notation, one can write the equivalent arch resistances as: 



^ = Mi(q^ + %\ - 

% = + MgCq^ + q^h^ - q^h^) 

\rtiere q = (l/L^)(a + c) 

s. 

q-jj = (l/2L^)(d/c - b/a) 

q^ = (c/L^) 

q^ = (l/2L^)(d/c + b/a) 

3.5 Rib Shortening 



The results of the dynamic analysis developed else\diere in this 
thesis may be modified xising the "rib-shortening" correction derived in this 
section in order to give an indication of the significance of the effect of 
direct compression in the eurch rib on the total displacement of the arch. Hie 
necessity of considering direct compression effects separately was discussed 
earlier (cf. Sections 2.1 and 3*l)- Essentleilly, the model ceuanot accept axieLL 



changes in the length of the bars and retain its characteristics as a two- 
degree -of -freedom system. While this is a definite limitation, it is not 
serious in a study of lateral deformations as long as the arch rib has suf- 
ficient capacity to withstand the thrust in the arch without general yielding 
or buckling. 

The approach used in this section is to consider the response of 
the arch to a uniform compression loading. This response is then superimposed 
on the lateral response in order to approximate the total response of the arch. 
This is an approximation in that these effects are in fact interrelated euid 
cannot be entirely separated, althou^ if one considers the relative compres- 
sibilities of the arch and the soil, direct compression effects are not as 
important as might othen#ise be expected (cf. Chapter IX ). 

Note that the soil provides a high degree of damping to the response 
of the arch in this mode. This effect cannot be neglected as it tends to pro- 
duce an over-critically damped response. Further, the extremely short natural 
period of vibration of the arch-soil "structure" in this mode makes the rise 
time of the shock wave in soil significant. This also tends to reduce the 
magnitude of the response. Therefore, for the purpose of the approximate 
correction derived in this section, a dynamic magnification factor of unity 
is assumed for this motion. 

The arch loading does contain a significant "uniform compression" 
component because of the rapid envelopement of the arch by the shock wave 
6uid the relatively slow decay of overpressure with time behind the shock 
front. Therefore, it appears reasonable to consider this component of the 
loading in computing the thrust in the arch rib. In this study, the uniform 
compression loading on the arch is taken as a proportion of the average vertical 
ovearpressure acting on the arch (cf. Chapter Vl). This proportion depends 
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principally upon the ratio of the horizontal to the vertical overpressure 
in the soil, althou^ the fact that the crown is subjected to a hi^er value of 
overpressure than the average radial overpressvure on the arch must also be 
considered. 

Two additional assumptions are made in order to simplify the compu- 
tation of the rib shortening itself. First, the average thrust in the arch is 
assumed to be equal to that in a circuleir ring under uniform compression. 
Second, the arch is assumed to maintain its same general shape as it deforas 
under the uniform compressive loading. This is approximate, but then the 
nature of the correction derived in this section does not warrant a rigorous 
solution for rib shortening in the arch. As will be seen, these assumptions 
lead to simple expressions which give reasonable results and are therefore 
justified on this basis. 

Consider the arch to be loaded as shown in Figure 5 •8* The thrust 
in the arch, T, may then be determined by equating the sum of the vertical 
forces acting on the arch to zero. 

2Tsin 3 = 2Rsin ^ p' 

a 

T = p' R 
a 

t = |f;b ( 5 . 5 . 1 ) 

where p' = "uniform compression" component of the average overpressvtre 
loading the arch (cf. Chapter VI ). 

The change in length of any one of the bars of the model caused by 

the thrust is given by Equation (5.5-2). 

6L. = TL /aE 
J J 

5 L^ = |p^BLyAE (5.5.2) 

tL. = change in length of bar j due to thrust 
J 

L. = length of bar j 
J 



where 
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A = cross-sectional area of the arch rib (eissxDned constant) 

E = modulus of elasticity of the arch material. 

The model is shovn in its undeformed and deformed configurations in 
Figinre 3* 9* llie arch abutments are considered immovable (cf. Section 1.2); 
therefore, joint 0 remains stationary. For the model to maintain the same 
general shape as the bars shorten, let bar 2 rotate about point O' such that 
joint 1 moves in a direction perpendicular to the original direction of bar 2. 
This is shown in Figure '^.10. Now, from symmetry joint 2 must move vertically 
downward. Tlierefore, joint 2 may be considered first to perform the same 
motion as joint 1 and then to rotate and shorten until it is vertically below 
its original position. This is illustrated in Figvtre 

The displacement coordinates of joint 1, i.e. the "rib -shortening” 
correction for joint 1, may be written directly from Figure 3- 10* 

6^ = 6L^/s in ( 0^-02 ) 

= 6L^/sln 0^cos 6^-coe 0^sin 0^ 

= 6L^(L)^/ad-bc ( 3 . 5 . 3 ) 

6h^ = - 5^3 in 0^ 

= - (bI^)Lc/ad-bc (3.5.*^) 

6hj^ = + 6^cos 0g 

= + (6L^)Ld/ad-bc (3.5.5) 

The displacement coordinates for the crown may be written directly from 
Figure 5 . 11 . 

bhg = 0 ( 3 . 5 . 6 ) 

6h^ = bL^/sin 0^ 

= ( 5 L 2 )l/c (3.5.7) 

The "rib-shortening" corrections for the model in terras of the 

arch geometry and the uniform overpressure which is effective on the arch 



W'^ 
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may now be obtained by substituting Equation (3»5«2) into Equations (3 •5*^) 
through (3»5'7) and considering symmetry. 



5h 



1 



1 BL^c 1 

2 ^a [_AE(ad-bc)J 



bh^ = 0 




oh 



4 ■ 




5h 



5 



+ 



1 -I I BL^ j 

2 ^a |_cAE J 



oh6 



bh 



k 



(3.5.8) 



2 

Note that if the overpressure, p', is expressed in units of fo rc e/ (length) , 
the cross-sectional area of the arch rib, A, must be expressed in units of 
( length )^/ length . 

An alternate approach to the derivation of the "rib -shortening" 
corrections is to allow the havinches and the crown to move such that, if the 
elevation of the abutments is used as a reference, the ratio of the elevation 
of the haunches to the elevation of the crown remains constant. Expressions 
were derived on this basis and found to give results, for the nvunerical 
examples studied herein (cf. Chapter VIIl), which do not differ significantly 
from those given by Equations (3*5.8). Since the alternate expressions are 
somewhat more involved, they were discarded in favor of Equations (3- 5*8). 
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TABLE 3.1 

CONSTRAINT LOACn^GS 



m 



1 

3 

3 

4 



2,m 




0 



-2(f) 



0 





5,m o,m 




0 

+a 

0 



-a 



-2 0 

0 +a 

0 0 

0 +a 
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Figure 5.1. Constraint Loadings 
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Figure 3*'^» Load Concentrations on Model 





Figure 30* Diagram Used in Determining 
Joint Rotation Relationships 
















.'^v 



■n:i 



1 
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Figure 5*4. Diagram Used in Determining Moment-Rotation 
Relationships for Joints 




Figure 5* 5* Idealized Moment -Rotation Diagram for Joints 



I 



k2 




Figure 5*6. Free -Body Diagram of Typiceil Bars and Joint 




Figure 5»7* Free-Body Diagram of Typical Joint 



i 
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Figure 5 -S' Uniform Radial Load on Arch 




Figure 3* 9* Deformation of Model Due to Rib Shortening 




Figure 5 . 10 . Diagram Used in Determining Rib- 
Shortening Correction at Haxmch 







Figure 5 •11. Diagram Used in Determining Rib- 
Shortening Correction at Crown 
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CHAPTER IV 
SOIL RESISTANCE 

In the case of an underground arch tlie soil contributes significant- 
ly to the total resistance to nuclear bleist loading of the arch-soil "structure". 
As with the arch, this resistance can be conveniently considered in two peirts, 
viz. inertial resistance, or that due to the mass of the soil which moves with 
the arch, and deformation resistance, or that due to the presence of shear and 
frictional forces within that portion of the soil which deforms with the arch 
under blast loading. Deformation resistance is discussed here. Inertial or 
mass effects are considered separately in Chapter V. 

As discussed in Section 2.2, the concept developed in this study 
regarding soil resistance due to frictional and shearing forces in the soil is 
based upon consideration of the state of stress existing in the soil at any 
particular instant during the response of the structure. Specifically, in the 
derivations which follow, the limiting stress conditions of the "Active State 
of Plastic Equilibrium" and the "Passive State of Plastic Equilibrium, " as 
defined in Section 2 . 2 , are considered as well as the intermediate condition 
of "Earth Pressxrre at Rest. " In general, differences in the gravity forces 
acting on the arch for different states of stress in the soil, together with 
the difference between the blast loading exerted on the ground surface and 
that transmitted to the arch by the soil, are considered as due to "soil 
resistance." A linear interpolation is used between the "at rest" and the 
"plastic equilibri\mi" states, and eui unchanging soil resistance assumed for 
deformations greater them those required for "plastic equilibrium." This in 
effect assumes em elasto -plastic resistance function for the soil vdiich is 
consistent with the present state of knowledge of soil -structure behavior. 
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Since the model developed for the dynamic analysis (cf. Section 
contains three interior hinges, three regions, one adjacent to each hinge, are 
considered separately from the standpoint of soil deformation and resistance. 
Conditions are different at the haunches from those at the crovn. Further, 
since conditions at a ha\inch axe different depending upon vdiether the haunch 
has moved invard into the arch or outward into the soil, these cases are treated 
separately. For convenience, a haunch ^ich has moved inward into the arch with 
respect to its original position prior to the arrival of the shock wave, or with 
respect to its original position as modified by any "permement set" in the soil 
adjacent to the arch, is referred to herein as an "interior" haunch. A haunch 
which has moved outward into the soil with respect to the same reference posi- 
tion is referred to as an "exterior" haunch. This classification is particular- 
ly usefvil when studying "rebound, " or response of the system after the first 
maximum response has occurred (cf. Section 4.4). 

4.1 Interior Haunch 

As a haunch moves inward into the arch, the soil tends to follow it 
under the influence of the overpressure acting on the ground surface. Ihls 
tends to create an "active" state of stress within the soil surrounding the 
haunch. Further, since the blast loading is applied to the arch through the 
soil, the shearing and frictional stresses developed within the soil mass will 
assist the structure in resisting the load. 

To illustrate how this "soil resistance" is developed, consider the 
loading "sensed" by the arch as the arch deforms. As mentioned in Section 2.2, 
the pressure sixrface, or surface on \diich the soil pressure acts, at the side 
of the arch is assumed to coincide with the lower bar of the model and in 
addition to extend upward in the same direction from the haunch to the ground 
sxirface. In Figxu*e 4.1 P' represents the force on the side of the arch due to 
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the weight of the soil. It is assumed to act at the lower third point of the 
pressure svurface. P" represents the force on the side of the arch due to the 
"surcharge" of the blast overpressure acting on the ground sxirface. It is 
assumed to act at the midheight of the pressvure surface. Since classical 
earth pressure theory (Coulomb eissumptions, cf. Section 2.2)is being used to 
define the limiting stress conditions for the calculation of hydrostatic pres- 
sure ratios for the soil stresses, P' and P" are assumed to act at an angle to 
the direction of the pressure surface. This angle, 6, is the friction angle 
corresponding to the friction developed between the surface of the arch and 
the soil. 

Now consider the hinges of the model to be "locked" so that when the 
blast loading is applied, no deformation can take place. If no movement takes 
place, the state of stress in the soil cannot change, and it therefore remains 
in the "at rest" condition. For this condition of "no deformation" one can 
write expressions for the forces acting on the sides of the arch as: 



^OA ■ 2 ’^s‘^0 ^OA 



(4.1.1) 



where 



and: 

vhere 



^OA ~ ^ direction of P^ of the force acting on the 

side of tlie arch due to the weight of the soil and corre- 
sponding to the condition of earth pressure at rest 



7^ = weight of soil per unit volvnne 



d 



o 



total depth from tlie ground siirface to the arch abutment 






ratio of the soil pressure in the direction of P^ to the 
vertical soil pressvtre for the condition of earth pressvure 
at rest. 




- p ,d K^. 
^sj o OA 



(4.1.2) 



Pq^ = component in the direction of Pj^ of the force acting on the 

side of the arch due to the blast overpressure euid correspond- 
ing to the condition of earth pressijre at rest 



k& 



Pg . = average overpressvire acting at the ground surface on the 

ixatentieil "active" failure wedge in the soil (defined 

below). Note that since p . acts downward, it is considered 

sj ' 

negative in accordance with the sign convention used herein 
(cf. Section 1.2). 

Now consider the hinges of the model to be "xinlocked" and the 
model permitted to deform until a fully developed "Active State of Plastic 
Equilibrium" is reached. At this point frictional and shearing forces exist 
within the soil mass and the expressions for the forces acting on the side of 



the arch are: 







d^ K. 
o A 






where = force acting on the side of the arch due to the wei^t of the 

soil and corresponding to a fully developed "active" state of 
stress within the soil mass. 

K. = hydrostatic pressxire ratio for the "active" state of stress 
in the soil, 



and; P” = - (I..1.4) 

where P^ = force acting on the side of the arch due to blast overpressure 
and corresponding to a fully developed "active" state of stress 
within the soil mass. 

The results of the preceding discussion are presented graphically 
in Figure 4.2. The "no deformation" loewi components. Equations (4.1.1) and 
(4.1.2), are plotted along the vertical axis, i.e. corresponding to no rota- 
tion of the pressure sxirface. The "fully active" load components, Equations 
( 4 . 1 . 5 ) and (4.1.4), are plotted at a deformation corresponding to that 
amount of rotation of the pressvire surface required to reach the "Active 
State of Plastic Equilibrium" within the soil mass. A discxission of the 
amount of deformation required is contained in Sections 2.2 and 4.4. 

Because of the forces developed within the soil during deformation, 
Pj^ + P^ is less than Pq^ + Pq^. Hiis difference is defined herein as the 
"soil resistance" corresponding to an inward movement of a haunch sufficient 
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to develop fvilly an "active" state of stress in the adj8icent soil. An 
expression for the soil resistance corresponding to this "fully active" condi- 
tion can now be derived by taking the difference of each of the load components 
cons ider ed , i . e . : 

Let: and = maximum values of soil resistance corresponding to a 

fully developed "active" state of stress in the soil. 

Then, using Equations (4.1.1) and (4.1.3): 



K = PAa - Pa 
A OA A 






(4.1.5) 



and, using Equations (4.1.2) and (4.1.4): 



R. = PAa - P, 



‘A ‘OA ‘A 

^sj^o^A ^ Psj*^o^A 

"a = - p.A%a - ‘'a) 



(A. 1.6) 



The nature of the variation of the soil resistance between the 
"at rest" and the "fully active" states can only be inferred from field 
observations of test structures, laboratory model studies and a general con- 
sideration of the nature of stresses in masses of soil. No exact relationship 
is known to the writer. One possible relationship is shown qualitatively in 
Figure 4.5. This corresponds, in general, to Figure 2.2. For the purposes 
of this study the relationships between deformation and soil resistance will 
be approximated using linear functions. This is discussed in detail in 
Section 4.4. 



While Figures 4.2 and 4.5 are helpful in visualizing the development 
of the soil resistance concept used in this study, they are somewhat misleading 
in that the two components of soil resistance as expressed by Equations (4.1.5) 
and (4.1.6) cannot be added directly to one einother for insertion in the 
equations of motion for the model because they act at different points. In 
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order to combine these quantities into what will be termed herein the 

"effective" soil resistance at the haunch, the lower bar is considered as a 

beam simply supported at the haunch and at the abutment and "loaded" with the 

two soil resisteuices and R^. The "reaction" at the havinch is then taken 

as the soil resistance which is "effective" at the haunch. This concept is 

illxxstrated in Figure U.4. While this neglects moments at the haunch, it 

is sufficiently accvirste for this study. 

Also shown in Figure h.U is the soil wedge corresponding to the 

development of the "active" state of stress in the soil. Ihis wedge is 

booinded by the ground surface, the pressure surface and a plane passing upward 

n 

through the arch abutment at an angle of + g) with respect to the horizontal, 

where "P is the angle of internal friction of the soil. This concept follows 

directly from the use of classical earth pressure theory to define the "active" 

state of stress in the soil plus the assumption of a plane failvtre surface for 

the determination of K^. As with the arch itself, this wedge is considered 

herein to have a unit thickness in the direction of the longitudinal arch axis. 

The overpressure term, p ., in Equations (U.1.2), (4.1. U) and (4.1.6) is taken 

s J 

as the average overpressure acting at the ground surface on this failure wedge. 

Keeping in mind that R^ and R^ as given by Equations (4.1.5) and 
(4.1.6) are maximvmi values, the maximum effective soil resistance for an 
interior haunch may be obtained directly from Figure 4.4 as: 



Qgj = i cos(w + 6) [R;^(dy5) + 
^sj = ^ ^ 5 R"a) cos (u, + 6) 



(4.1.7) 



bd "I 



and: 



' ■ b ['*A 

«sj.3 ' - ^ ^''a) 



(4.1.8) 



51 



I where 

I I 



Q . = maximum value of the horizontal component of the soil 
^ resistance \rfiich may be considered effective at an 
interior haunch 



I 



A 



Q. 



sj+3 



maximian value of the vertical component of the soil 
resistance vdiich may be considered effective at an interior 
haunch 



a^ b ^ c ^ d 



generalized arch dimensions (cf. Figure 2.1 ) 



I 



inclination of the lower bar of the model with respect to 
the vertical 



6 = friction angle corresponding to the friction developed 
between the surface of the arch and the soil. 



d^ = depth from the ground sxurface to the arch abutment 

/ " maximum values of soil resistance components as defined by 

Equations (4.1.5) and (4.1.6) 



Substituting Equations (4.1.5) and $+.1.6) into Equations (4.1.7) 
and (4.1.8): 

«sj * SI <Sa - ■'aX'A - <‘‘•^•9) 

«SJ.5 ' - 51 ‘■'OA - 

These equations will be vised in Section 4.4 in developing expressions for the 
total soil resistance of the system. 

One further point to be considered here involves the hydrostatic 
pressure ratios for this case. Methods for determining are compared in 
Appendix A. Values for Kq^ cannot be calcvilated using classical earth pressure 
theory, althou^ values have been measured in tests. In the numerical, 

solutions of this study an observed value of the ratio of horizontal to 
vertical soil pressures coi*responding to the "at rest" condition is assvmed. 

It is then necessary to "resolve" this ratio from that corresponding to hori- 
zontal pressures to that corresponding to pressures paredlel to the line of 
action of P'. A diagram of the relationships involved is given in Figure 4.5. 
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Referring to this figure, one may write: 



1 2 
^ Ka 
2 s o A 



2 V/ 



= cos 



1 + k 



j^(w 



+ 6 ) 



arctan f 
0-1 



From which: 



=\/~1 + k^ cos |^(w + 6) - arctan (U.l.ll) 

\rtiere k = observed value of the ratio between the horizontal and 

vertical soil pressures for the condition of earth press\rre 
at rest. 



4. 2 Exterior Ha\inch 

Ihe soil resistance developed at an exterior haunch may be derived 
in a manner analogous to that used in the previous section for the interior 
haunch. The essential difference is that movement of the arch outward against 
the soil tends to produce a "Passive State of Plastic Equilibrixm" (cf. Sec- 
tion 2.2) within the soil. As in the previous case, the "soil resistance" 
will be taken as the differences in 'Uie loads transmitted to the arch by the 
soil in the "locked" and "unlocked" conditions. The forces acting on the arch 
for this case are shown in Figiare 4.6. 

If the arch is now considered to be "locked" so that no deformation 
can take place, hence no frictional forces develop in the soil, and the bleist 
loading applied, one can write: 



\diere 



1 p 

pf ^ ZL y ir 

OP " 2 's o OP 



(4.2.1) 



PqP = component in the direction of Pp of the force acting on the 

side of the arch due to the wei^t of the soil and correspond- 
ing to the condition of earth pressure at rest, 

Kqp = ratio of the soil pressure in the direction of Pp to the 

vertical soil pressure for the condition of earth pressure 
at rest, 
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and: 

where 



OP 



p;j%^op 



(4.2.2) 



PqP = component in the direction of Pp of the force acting on the 
side of the eirch due to the blast overpressure and corre- 
sponding to the condition of earth pressure at rest, 



p' = average overpressure acting at the ground surface on the 
potential "passive" soil failxire wedge (defined below). 

Note that since p' . acts downward, it is considered negative 
in accordance with the sign convention used herein (cf. 

Section I. 2 ). 

If the structure is now "unlocked" and permitted to deform until the 
"Passive State of Plastic Equilibrium" is reached, frictional forces will 
develop within the soil mass and the forces acting directly on the side of the 
arch will be: 






(4.2.3) 



where Pp = force acting on the side of the arch due to the wei^t of the 
soil and corresponding to a fully developed "passive" state 
of stress within the soil, 

Kp = hydrostatic pressure ratio for the "passive" state of stress in 
in the soil, 

and: P^ = - P^jd^Kp (4.2.4) 

where Pp = force acting on the side of the arch due to the blast 

overpressure and corresponding to a fully developed "passive" 
state of stress within the soil. 

These load components eire plotted in Figvire 4.7. The "locked" 
condition is represented by = 0, i.e. no rotation of the pressure surface. 

The "unlocked" condition is represented by i.e. sufficient rotation of 

the pressure surface to develop fully the "Passive State of Plastic Equilibrivim" 
within the soil (cf. Section 2.2). 

As in the case of the interior haunch, the differences between the 
"locked" and "unlocked" conditions are taken as the "soil resistance." 

Referring to Figure 4.7^ these differences may be expressed by subtracting 
Equations (4.2.1) and (4.2.2) from Equations (4.2.5) and (4.2.4). 
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Thus, if: and R^ = maximuin values of soil resistance corresponding to a 

fully developed "passive" state of stress in the soil, 

Then: 

= I - I 

And: R^ = Pp - Pqp 

^ - ^ 0 ?'^ ('^• 2 . 6 ) 

Since the load components P' and P" are assumed to act at different 
points on the pressure surface (cf. Figure 4 . 6 ), the corresponding soil re- 
sistances R^ and R^ must be combined in a similar manner to that used in the 
case of the interior haunch. Again considering the lower bar of the model as 
a simply supported beam acted upon by the soil resistances R^ and R^, the 
effective soil resistance components at the haunch are (cf. Figure 4 . 8 ): 







i cos(t. - [r^ (^) + r; (^)j 

d 

(2R^ + 3R^) cos(6 - w) 



( 4 . 2 . 7 ) 



And: 






- r bd bd -j 

- i =ln(a - -) [r^ (3;^) . R; (^)J 

d 

- (2R^ + 3R^) sin(6 - w) 



( 4 . 2 . 8 ) 



\diere 




maximxmi value of the horizontal component of the soil 
resistance which may be considered effective at an 
exterior haionch. 






maximum value of the vertical component of the soil 
resistance which may be considered effective at an 
exterior haunch. 



a,b, c,d 



generalized arch dimensions (cf. Figure 2.l), 
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w = inclination of the lower bar of the model with respect 
to the vertical, 



5 = friction angle corresponding to the friction developed 
between the surface of the arch and the soil. 



Substituting Equations (4.2.5) and (4.2.6) into Equations (4.2.7) 
and (4.2.8): 



Note that the soil failure wedge now considered (cf. Figure 4.8) is 
bounded by the ground surface, the pressvire siirface and a plane throti^ the 



where is the angle of intemal friction of the soil. A unit thickness of 
wedge parallel to the horizontal arch axis is considered in this thesis. The 
assumptions leading to this wedge are the same as those leading to the defini- 
tion of the "active" soil failure wedge described in Section 4.1, i.e. the 
iise of classical earth pressure theory (Coulomb's assumptions) and the assump- 
tion of a plane failure svirface in the determination of the hydrostatic pressure 
ratio corresponding to full development of the "passive" state of stress within 
the soil (cf. Appendix A). 

Equations (4.2.9) and (4.2.10) are used in Section 4.4 in deriving 
expressions for the total soil resistance of the system. In order to use 
these equations, one must be able to calculate the hydrostatic pressure ratios 
for the soil pressures corresponding to the "at rest" and "passive" states of 
stress. Methods of finding Kp are compared in Appendix A. K^p is found by 
"resolving" obseirved hydrostatic pressure ratios relating horizontal to 
vertical soil pressures to a ratio relating soil pressures acting parallel to 
the line of action of Pp to vertic£il soil pressures. The technique used is 




(4.2.10) 



(4.2.9) 



arch abutment inclined 




with respect to the horizonted. 
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analogous to that described in Section 4.1 for finding A diagram of the 
relationships involved is given in Figure 4.9. Referring to this figure, one 
may vrite: 



From which: 
•vrtiere k 

o 



1 7 d^ , 1 + 

2 s o V o 



= cos l^arctan 



i t (6 - u) j 

o j 



K, 



OP 



1 + k cos 
o 



arctan 

o 



.)j 



(4.2.11) 



observed value of the ratio between the horizontal and vertical 
soil pressvires for the condition of earth pressure at rest. 



4.5 Crown 

As the soil-arch "structure” deforms in response to a blast loading, 
shear and frictional forces eire developed in the soil near the crown. Ihe 
definition of these forces is not an easy matter. In evaluating the relative 
importance of these forces, consider that for the sheillow depths of cover 
included in this study (cf. Section l.l) the volume of soil at and adjacent to 
the crown is relatively much smaller theua the vol\ane of the soil wedges defined 
in the preceding sections, affected by lateral displacements of the sides of 
the arch. Thus, the maximum resultant value of these soil forces will 
necessarily be relatively smaller at the crown than at the sides of the arch 
because of the smaller volume of soil involved. Further, some of these 
resisting forces serve only to distribute the applied load away from the 
crown to other regions of the arch, thus contributing no net resistance to the 
system. A portion of these forces can, however, be assumed logically to act 
to transmit some of the applied load away from the arch. This latter portion 
of the soil forces does assist the arch in resisting an applied load but is 
not readily defined explicitly. 
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An attempt to include at least part of this "effective" portion of 
the soil resistance near the crown in the total soil resistance of the system 
is made by continuing upward the pressiure siarfaces at the sides of the arch, 
in the same direction as the lower bars of the model, until they intersect 
the ground surface. This is not a rigorous approach, but it is believed to 
be sxifficiently accurate for the present. Since the forces involved are 
relatively small, it would seem to follow that their influence on the response 
of the system as a ^ole should be relatively small as well. 

Die effect of the inertial resistance of the soil near the crown 
is discussed in Chapter V. 

As a result of the foregoing discussion, the effective soil resistance 
to be placed on the model at the crown is zero. Daus, the following expressions 
may be written: 

Q32 = 0 (J+.5.1) 

where = horizontal component of the effective soil resistance at 

the crown. 



and: 



= 0 
s? 



(4.5.2) 



where ^ " vertical component of the effective soil resistance at the 

^ crown. 



4. 4 Development of Equiveilent Soil Resistance 

Since the model is a two -degree -of -freedom system (cf. Section 5*l)^ 
the effective soil resistances derived in the preceding sections must be 
combined such that their net, or "equivalent, " effect on the model is defined 
in the direction of the independent displacement coordinates. This is done 
in this section for the set of independent coordinates already selected, viz. 
the horizontal displacements of the haunches. Thus, the object here is to 




I».t a 
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define the resultant resistance offered by the soil to horizontal movements 
of the haunches. Note that this is not the same as the soil resistance 
components derived in Sections U.l and k.2 since now the entire system is 
considered simultaneously. 

The technique used herein is simply to apply the effective soil 

resistance components to the model as a system of loads. Then, by applying 

the "Load Relationships" derived in Section 5*2.1 (Equations (5*2. 5))^ the 

equivalent soil resistances in the horizontal direction at the haunches are 

determined. This may be expressed algebraically by substituting the effective 

soil resistance components, Q ., for the load components, H., in Equations 

s j j 

(3.2.5)* Thus: 



«sl - «sl ^ I ^ |(f - !> «S5 



(i.li.l) 



\diere = equivalent soil resistance at the windward haunch, 

Qgj = equivalent soil resistance at the leeward haunch. 

Since the expressions for the components of the effective soil 
resistance at the joints were derived in the preceding sections in terms of 
the limiting or maximum values, it is convenient here to write meiximum values 
for the equivalent soil resistance. Further, since the haunches may be 
displaced either inward or outward, it will be convenient to have two sets of 
equations, one for interior haxmches and one for exterior havinches. Ihese 
equations may then be used in various combinations for dynamic analysis as 
the configuration of the system dictates. 

For the case of both haunches displaced inward, the maximum equiva- 



lent soil resistance may be foxmd by substituting Eqviations (4.1.9), (4.1.10), 
(4.3.1) and (4.3.2) into Equations (4.4.1). This gives: 
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%1 = SI ^ (|)6in(u H- b)j 

W + 6)J 



2 

%i ‘ si ■ 5Pb3* 



+ 6) + (^)sin( 

Q. 



Note that the preceding expressions differ only in the overpressure term. 
Thus, one may write, in general; 

A d^ I V, 1 

= SI ^^OA • ‘^a)^’'s‘^o - 5Psj) + &)J 



A 



(4.4.2) 



where Q . = maximum value of the horizontal component of the e>^uivalent 
soil resistance at haunch j when the haunch is displaced 
inweird, 

p . = average overpressure acting at the ground surface on the 

"active" soil failure wedge adjacent to joint j (cf. Figure 4.4). 

Similar expressions may he written in the case ^ere hoth haunches 
are displaced outweird. Substituting E-juatlons (4.2.9)> (4.2.10), (4.3.1) and 
(4.3.<) into Equations (4.4.l): 

A d^ 

^8l = SI ■ ^OP^^^s'^o ■ 5 Ps 1^ (|)sin(<^ - &)j 

A d^ 

^S3 = SI (|)6ln(« - 6)J 



Note that these expressions also differ only in the overpressure term. Thus, 
in general: 

.2 



«aj = ^ <>S> - ''oP^Vo - [=“=<“ - 



(It.'l.j) 



where Q' . = maximum value of the horizontsil component of the equivalent 
^ soil resistance at haxmch j when the haunch is displaced 
outward, 

p' . = average overpressure acting at the ground surface on the 
"passive” soil failure wedge adjacent to joint j 
(cf. Figure 4.8). 



Note that both Equation (k.k.2) and Equation (4.4.;5) express the 
equivalent soil resisteince at a haunch solely in terms of the components of 
'' the effective soil resistance acting at that haunch. This permits the haunches 
I to be considered separately and reduces the number of cases which must be con- 
I sidered in the remainder of this discussion to two^ viz. the interior haunch 
I and the exterior haunch. The resulting expressions can then be applied to any 
conceivable configuration of the system. 

An impKjrtant aspect of the soil resistance developed eis the system 
responds to a dynamic load involves the amounts of deformation required to 
1 develop fully the "active” and "passive" states of plastic equilibrium. In 
Section 2. 2 reference was made to a number of estimates of these 

parameters. There is general agreement that in the case of retaining wedls 
backfilled with cohesionless soil, a rotation of the wall about its base 
corresponding to a movement away from the backfill at the top on the order 
of O.I 5 & of the wall hei^t is sufficient to reduce the state of stress in the 
soil behind the wall from the "at rest" to the "fully active" state. In the 
case of the development of the "passive" state of stress, quantitative esti- 
mates are generally lacking, eilthoxi^ it appears certain that the required 
movements are greater than for the "active" case (cf. Figure 2.2). A 
technique for estimating this required movement is outlined below. Note that 
the pressure surfeices assumed in deriving the effective soil resisteinces 
(cf. Section ^.1 eind k.2) are similar to the pressvure surfaces behind retaining 
walls. Therefore, the significant deformations of the model for the develop- 
ment of soil resistance are the rotations of the lower bars, since these define 
the positions of the pressure surfaces. 

In this study llie idealized soil resistance -displacement function is 
obtained by assuming a linear buildup of soil resistance with increasing dis- 
placement until a maximxroi value is reached. The soil resistance is then 
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assumed to remain at this maximum value for further displacements in the same 
direction. This results in a so-called "elasto -plastic" resistance diagram 
as illustrated in Figure 4.10. 

Thus, the soil resistance-displacement function for an interior 
haunch may be defined by using Equation (4.4.2) to find the maxim\mi value of 
the soil resistance, and the assumption that this veilue of soil resistance 
corresponds to an inward rotation of the lower beir, of 0.001. Ibis is 
illustrated in Figxire 4.11. The range of possible values for the soil resist- 
ance corresponding to each value of displacement results from the overpressvure 

tenn, p ., in Equation (4.4.2). 
s j 

The maximum value of soil resistance for an exterior haunch is given 
by Equation (4. 4. 5 ). Since quantitative estimates of the corresponding 
deformations are lacking, the assumption is made that the slope of the idealized 
soil resistance -displacement function in the elastic range is constant for any 
given value of overpressure acting on the soil failure wedges. For the specific 
examples studied herein, this deformation corresponds to a rotation of the 
pressure surface on the order of O.OO5, which appears to be in reasonable a^ee- 
ment with existing theory and experimental evidence. ^ Note, however, 
that this value has not yet been determined definitely. An exi>erimental in- 
vestigation of the deformation corresponding to full development of the 
"passive" state of stress in the soil is certainly needed (cf. Chapter IX ). A 
complete, idealized soil resistance -displacement dieigram is shown in 
Figure 4.12. 

Expressions for the "yield" deformations of the soil may be derived 
as follows: 

Let: = rotation of the lower bar of the model as the 

structure responds to a blast loading. 
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i = rotation of the lower bar of the model corresponding 
^ to full development of the "active" state of stress 
within the soil adjacent to the haimch^ 

I' = rotation of the lower bar of the model corresponding 
^ to full development of the "passive" state of stress 
within the soil adjacent to the haunch. 

Then, from Figure 4.12: 

A 

^Sj 

Substituting Equations (4.4.2) and (4.4.5) into Equation (4.4.4): 






V ~ i 
^ty ^ty 



si - "opXVo - 5 p^j) f=o»(“ - '>) * - ^)] 

- r V T 

si <’'0A - "a)' V o ■ * ®)J 



(4.4.5) 

Now, the difference between p^^ and p^^ is very small except for 

the brief Interval dvu*ing which the shock front traverses the structure and 

the adjeu:ent soil wedges. This is due to the relatively small dimensions of 

the structure as compared with the spatial dimensions and slow dece^r with time 

of the nuclear blast wave. P\irther, the large differences which do exist 

between p . and p' during the shock front transit, ■vdxen pressure buildup 
6 J s J 

occurs, are not significant here because dxiring this internal the structure 
does not have time to respond sufficiently to develop significant soil 
resistances. Thus, Equation (4.4.5) niay be reduced to: 




r"s> - ^op- 




cos(w - &) + (b/a)sin(w - 5) 


1 

1 

< 

uP 
1 




cos(w + b) + (b/a)sin(w + 6) 



(4.4.6) 



Now, note that: 

cos(u) - o) + (b/a) sin((*> - 6^ 1 _ cos (u) - 6) + tan <*) sin(<<J - b) _ . 
cos(w + 6) + (b/a) sin(w + 6)J “ cos (u + b) + tan w sln(w + b) ~ 

(4.4.7) 
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Therefore ; 



'iy 



) 



(4.4.8) 



Alternatively, the idealized soil resistance -displacement function. 
Figure 4.12, may he defined using either the maximum soil resistances or the 
"yield" displacements for the soil masses as derived above, and the slope of 
the "elastic" line, or "soil stiffness." For an interior haunch the "soil 
stiffness" may be expressed as: 



ic . = ^ 






where K , = stiffness" of the soil during the development of the 
"active" state of stress. 

Substituting Eqviat ion (4.4.2) into Equation (4.4.9): 



K 



(K^ 



6j '‘"OA 



K^)(ygdo - 5Pgj) cos(w + 6) + (^)sin(w + 6)j 



or: 



K . = k' - k"p . 

sj 



(4.4.10) 



where k' = " (K, 



7 d-' 
s o 

6& i , ^*^0A 

ty 



- K^) j^cos(w 



+ 6) + (f)sin( 

cx 



w + 5)1 



k" . (^) k- 

s o 

While for an exterior haunch, the"soil resistance" may be expressed as: 



Q' . 

ic' = -4^ 



(4.4.11) 



K' . = "stiffness" of the soil dtxring development of the "passive" 
^ state of stress. 



^dlere 






L 
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Substituting Equations (4. 4 . 3)7 (4.4.7) (4.4.8) into 

Equation (4.4.10) gives: 



or; 




SI 0^, - KopM Vo 




,'S> 



K, 



OP. 



^OA - K/ 




- 3Pg^) [cos(u) - &) + (^)sin(fa> - 6)] 

cos(t»> - S) + (b/a)sin(to - b) 
cos(tJ + b) + (b/a)sin(w + o) 




(4.4.12) 



Note that the development of the tvo expressions for "soil stiffness, " 
K . and K' (Equations (4.4.10) and (4.4.12)), is consistent with the assumption 

8 j s J 

of a constant slope for the "elastic" line of the soil resistance -displacement 

function for a given value of average overpressure, i.e. for p^j = p^j . 

Equations (4.4.10) and (4.4.12) express the soil stiffness in terms 

of a " Q .vs i. " diagram (cf. Figure 4.12). It is sometimes more convenient 
s J X 

to deal directly in terms of the horizontal displacements of the haxxnches 
rather than with the rotations of the lower bars. Now, from Equations (5.2.9): 

= 60^ = - - h^a. (4.4.13) 

Therefore, the soil stiffnesses may be expressed as: 








(4.4.14) 



and: 



(-) K' , 
a sj 





(4.4.15) 



At this point one should consider the resistance offered by the 
soil as the structure "rebounds," or reverses its direction of motion, during 
response to a blast load. Keeping in mind the fact that the overpressure 
tends to force the soil into intimate contact with the arch, the variation in 
soil resistance dviring rebound is assumed to be as shown in Figure 4.13. 
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It is now possible to examine more closely the classification of 
haunches according to their displacement into "interior" and "exterior" 
haianches eis was done at the beginning of this chapter. Referring to Figure 
4 . 15 , the original position of the haunch is represented by the origin of the 
coordinate system. The effect of the "permanent set" is to transfer the 
origin of tlie coordinate system to allow for any "plastic" deformation which 
may take place in the soil during the response. This is, of coxurse, an 
approximation of the behavior of the soil, but it is consistent with the 
behavioral assumptions made previously in this chapter. Note that a reversal 
in the direction of haunch movement does not cause an immediate reversal In 
the direction of the soil resistance. There is some evidence to indicate^^^^ 
that the reversal of the direction of the soil resistance will occur more 
rapidly with reversal of the direction of movement than Figure 4.15 indicates; 
however, further investigation of this behavior was not included in this 
study (cf. Chapter IX ). 

Consistent with the sign conventions xised elsewhere in this thesis, 
soil resistances acting to the left are positive, and conversely. 



Load on Arch 
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Figure 4.1. Forces Acting on Pressure Surface for "interior" Haunch 




Rotation of Lower Bar 



Figure 4.2. Load -Displacement Diagram for "interior" Haunch 



6 ? 




Figure Soil Resistance -Displacement Diagram for 

"interior" Haunch 




Figure U.4. Effective Soil Resistance for "interior" Haunch 
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Figure 4 . 5 . Diagram Used in Determining Kq^ 



» 




Figure 4.6. Forces Acting on Pressure Surface for "Exterior” Haunch 



Load on Arch 
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Figure 4 . 7 . Load -Displacement Diagram for "Exterior" Haunch 




Figure 4.8. Effective Soil Resistance for "Exterior" Haunch 



( 




Equivalent Soil Resistance 



TO 




Figure U. 9 . Diagram Used in Determining K^p 




Rotation of Lower Bar 



Figure 4.10. Partial, Idealized, Soil Resistance - 
Displacement Diagram 



Equivalent Soil Resistance 
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Figure 4.11. Partial, Idealized, Soil Resistance -Displacement 
Diagram Showing Effect of Overpressxire 




Figure 4.12. Complete, Idealized, Soil Resistance -Displacement 
Diagram 
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Figure 4.1j5. Assumed Variation in Soil Resistance During Rebound 
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CHAPTER V 
MASS EFFECTS 

The mass of a system s\xbjected to dynamic loading directly eiffects 
the inertial resistance of the system and thus has eui important effect on 
both the magnitude and the time history of the response of the system. It 
is the pvirpose of this chapter to list and discuss the eis sumptions made 
regarding the mass of the arch-soil "structure” and to derive the relationships 
necessary to incorporate inertial resistance into the eqiiations of motion of 
the system (cf. Chapter VII ). 

In general, the approach used here is to consider the meusses of arch 
and soil vhich may logically be assumed to move as the system responds. The 
inertial resistances of these masses are then calculated and applied to the 
model as a system of loads. The resultant effect of these ineirtial forces 
is taken as the equivalent inertial resistance of the model. As in the case 
of the equivalent deformation resistances determined earlier, the equivalent 
inertial resistance is determined in agreement with the choice of independent 
displacements for the model. 

5.1 Maas Concentrations on the Model 

Since the model has three interior joints, masses may be placed at 
each joint in such a manner as to approximate in the most realistic possible 
fashion the inertial resistance of the actual system in the region adjacent 
to that joint. Die model itself imposes no inherent restrictions on the choice 
of the masses which may be placed at the joints. Hie mass of the arch may be 
handled quite readily by concentrating at each joint the mass of the arch rib 
for a one-half bar length to either side of the joint. Note that this results 
in the concentration of the mass of the lower half of each of the lower bars 



I 












74 



at the arch abutments. Since the abutments are assumed immovable (cf. Sec- 
tion 1 . 2 ), this portion of the mass of the arch is thus assumed to have no 
significant effect upon the response of the system. As pointed out in 
Section 2.5, such assumptions are not uncommon in dynamic analyses and give 
reasonable results. 

The mass of soil \diich may logically be assumed to move with the 
8LTch is more difficult to define. In general, the mass of soil moving with 
a joint depends upon: 

1. The location of the joint, e.g. the mass of soil moving with 
the crown will be different from the mass of soil moving with one of the 
haunches. 

2. The direction of movement of the Joint, e.g. the mass of soil 
eiffected by the inweird movement of a ha\inch is different from that affected 
by its outward movement. 

5 . The depth of burial of the arch, e.g. the mass of soil moving 
with the crown will tend to increase as the depth of cover increases. 

Further discussion of the latter two considerations is weirranted 
here. With regard to the direction of movement of the joint, note that this 
approach is somewhat different from that used in Chapter IV in developing the 
soil resistance due to deformation. The effect of both approaches is the 
same during initial response. However, once rebound occurs, the "deformation" 
resistance of the soil is assvuned to decrease lineeirly to zero before building 
up in the opposite direction (cf. Section 4.4). On the other hand, by con- 
sidering only the direction of movement when choosing the effective soil mass, 
the start of rebovind causes an immediate change in the magnitude of the mass 
of soil considered to move with the system. This apparent contradiction is 
not as seriovis as it might first appear, since the start of rebound corresponds 
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to a maximum displacement, hence zero velocity, of one of the two independent 
motions of the system, so that there is actvially no momentum 'vdien the mass 
assumption is changed. The expressions derived below permit the mass assump- 
tions to be changed independently for each degree of freedom of the system, 
since in genereil the independent motions will not reverse their direction of 
movement, i.e. "rebound," at the same time. 

With regard to the depth of burial of the arch, it has been suggested 
in the case of rectangular, buried structures that a depth of cover greater 
them the span length should not be used in estimating the mass of soil moving 
with the structure. This is a logical restriction and is incorporated 
into the derivations in this chapter. 

In deriving the mass concentrations on the model consider first that 
these concentrated masses will eiffect the inertial forces developed at each 
of the joints as the system responds to a blast loading. Second, the magnitude 
and direction of these inertial forces is variable but can be determined at 
any time dxiring the response. Now, instead of dealing with the res\iltant 
inertial forces at the joints, consider instead the horizontal and vertical 
components of these forces. The mass concentrations on the model vdiich are 
derived in the remainder of this section correspond to these horizontal and 
vertical components. This approach is more convenient than dealing with the 
masses corresponding to the resultant inertial forces at the joints. Note 
that this concept is not based upon horizontal and vertical "components" of 
the masses themselves. 

The masses at each joint corresponding to the horizontal and vertical 
components of the inertial forces at the joint may be expressed as: 

j sj aj 



( 5 . 1 . 1 ) 



T6 



\rtiere 



m. = mass concentration at Joint J conresponding to the 
^ horizontal, or the vertical, component of the inertiaO. 
force at joint j. 



m . = soil contribution to m . . 
sj J 

m . = arch contribution to ra.. 
aj J 



As mentioned previously, the mass contribution of the arch at each 



joint is the mass of the arch rib for half a bar length to either side of the 
joint. Hxis may be expressed as: 



\diere 



m 



7 h 
a a 



aj 2g '“j-1 * ^j^ 



( 5 . 1 . 2 ) 



m . 
aj 



8 

L. 

J 



= mass contribution of the arch at joint j 

= wei^t per unit voliane of the material from >rtiich the 
arch is made 

= thickness of the arch rib 
= acceleration due to gravity 

length of the bar to the left of joint j 
length of the bar to the right of joint j. 



The masses of soil associated with the horizontal components of the 
inertial forces at the joints are shown in Figure 5.1. While both of the 
havtnches are free to move either inward or outward. Figure 5.1 is presented, 
for convenience, as though the left-hand havinch, joint 1, were moving inward 
and the ri^t-hand haunch, joint 3 / moving outward. Note that in order for 
the dimensions of the soil blocks not to exceed the "corresponding span" of 
the model, restrictions have been placed on two of the dimensions as follows: 



= (| . c + d^) < (b + d) (5-1-3) 

- <1 * i “ (5.1.1.) 



vrtiere the dimensions are defined in Figure 
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\rtiere 



vhere 



vAiere 



Note also 



vhere 



•vrtiere 



Referring to Figure ^. 1 , one may write: 
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h + d> 




a+c \ /d\ 
2 ^' 2 ^ 




(5.1-5) 



^*°sl^in ° contribution of the soil at joint 1 corresponding 

to the horizontal component of the inertial force acting 
at joint 1 when the joint is moving inward 



".2 = ¥ <^.2 - 



(5.1.6) 



m ^2 “ mass contribution of the soil at joint 2 corresponding to 

the horizontal component of the inertial force acting at 
the crown for both inward and outward movement of the 
crown 



^"*s3 ^out 
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2g 









(5.1.7) 



(m^,)^ ^ = mass contribution of the soil at joint 5 corresponding 
to the horizontal component of the inertial force at 
joint 3 "when the joint is moving outward. 



that: 



(m -. ) . = (m , ) . 

s3 in si in 



( 5 . 1 . 8 ) 





= mass contribution of the soil at joint 3 corresponding 
to the horizonted component of the inertial force at 
joint 3 when the joint is moving inward 



(®i) + = (“!,) . 

si out s;) out 



(m T ) J. 

si out 



(5.1.9) 

= mass contribution of the soil at joint 1 corresponding 
to the horizontal component of the inertial force at 
joint 1 \dien the joint is moving outward. 



The masses of soil associated with the vertical components of the 



inertieil forces at the joints are shown in Figure 5-2. Note that the 
restrictions stated by Equations ( 5 . 1.3) and (5.1.4) also apply here. 



1 

I 
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Referring now to the figure: 

7 



m 



s4 



' i - 1)(¥) - 



where 



m = mass contribution of the soil at joint 1 corresponding to 
the vertical component of the inertial force at joint 1 
for both inward and outward movement of the joint 



m 



s5 



7 d 

— (d ^ - cA) 
g m2 



(5.1.11) 



where 



m c = mass contribution of the soil at joint 2 conresponding to 
^ the vertical component of the inertial force at joint 2 



for both inward and outward movement of the joint 



m 



s6 



= i * ("na - 



a+c 



ml 



)( 






(5.1.12) 



where 



m^g = mass contribution of the soil at joint 5 corresponding to 



the vertical component of the inertial force at joint 5 
for both inward and outward movement of the joint. 



■ |)(¥) * 



Now substituting Equations (5.1.2) and ( 5 . 1.5) through ( 5 . 1.12) 
into Equation ( 5 . 1 . 1 ) : 

+ =ot(J + |)j t 7hL/s 
^"l^out" 2g ['^1111^2^ * ^■aml ■ 

. 4a - t>j a- 

“2 

^”*5 ^out ^out 



m. 



Dl, 



5 

“6 



7 

s 

" 2g 
= 

= 



" (■aal ‘ " ‘'ami ' " "a^^" 

( 5 . 1 . 15 ) 



a+c \/d> 
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Equations (5«1.13) are expressions for the effective mass con- 
centrations at the joints of the model. Ihey are used in the next section in 
deriving expressions for the equivalent masses of the two -degree -of -freedom 
system. 

5.2 Development of Equivalent Mass 

It is now necessary to combine the inertial resistances of the 
arch-soil "structxire" in such a manner as to obtain their resultant effect on 
the response of the model. Again, the expressions which follow are developed 
considering the horizontal displacement components of the havinches as the 
independent variables. 

The technique followed here consists of using the mass concentration 
expressions derived in the preceding section to write the horizontal and 
vertical components of the inertial forces acting on the model at each of the 
joints. From the fact that in the response of any one of the masses the com- 
ponent of its acceleration in a given direction is proportional to the correspond- 
ing displacement component, the accelerations are related to one another using 
the "Displacement Relationships" derived for the model in Section 3*2.2. The 
inertial forces are then applied to the model and their resulteint effect in 
the horizontal direction at the haunches found using the "Load Relationships" 
derived in Section 3*2.1. 

Thvis, using Equations ( 3 . 2 . 7 ) and the principle stated above, one 

• « 

may relate the acceleration components of the masses, h. as: 

J 

h 2 = + I (1 - bc/ad)(ii^ + h^) 

= - (b/a)h^ 

h^ = + I (d/c - b/a)(h^ - h^) 
h6 = + (b/a)h^ 



( 5 * 2 . 1 ) 
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where = horizontal component of the acceleration at Joint 1 

iig = horizontal component of the acceleration at joint 2 
h^ = horizontal component of the acceleration at joint 5 
= vertical cmponent of the acceleration at joint 1 
h^ = vertical component of the acceleration at joint 2 
hg = vertical component of the acceleration at joint 5* 

Note that the acceleration components chosen for reference, i.e. h^ and h^, 
correspond to the horizontal displacement components of the haunches, which 
have alreeidy been chosen as the independent displacement coordinates for the 
dynamic analysis. 

Using Equations (5.1.15) and (5.2.1) it is now possible to write 
the horizontal and vertical components of the inertial forces acting on each 
joint of the model. Since there is nothing to be gained by rewriting here 
the mass concentration expressions given in Equations (5.1.15)^ the notation 
will be simplified by merely indicating the proper mass concentration, raj, to 
use for each inertial force component. Thus, the inertial force components, 
Zj, are: 

^ 

\ ~ ■ Wad)(hj^ + hj) 




= (m^/2)(d/c - b/a)(hj^ - h^) 

Zg = m^(b/a)h^ 

5y applying the inertial force components, Zj, to the model as a 

system of loads, the equivalent inertial force in the horizontal direction 

at the haunches, Z , may be fovind, using Equations (5*2.5)# asi 

J 
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* I (f - 

^ - 1 ^ - 1 (f - - (|)% 



( 5 . 3 - 5 ) 



Now let fflj ^ = equivalent mass of the system in the horizontal 



direction at haunch J. Then: 



Z. = m. h, 
J J 



(5.2.4) 



Substituting Equations (5.2.2) and (5.2.4) into Equations 



( 5 . 2 . 3 ): 



“l,m^l “ ™1^1 (“2^^)(l - bc/ad)^(h^ + h^) 

+ mj^(b/a)^^ + (m^/4)(d/c - b/a)^(h^ - h^) 

= (m2/4)(l - bc/ad)^(h^ + ii^) + m^h^ 

- (m^/4)(d/c - b/a)^(h^ - h^) + m^(b/a)^^ 

Rearranging terms gives: 

“l,m ” |_“l + m|^(b/a)^ + (m^/4)(^ " j 

- j^-(m2/4)(l - bc/ad)^ + (m^/4)(d/c - b/a)^j 



( 5 . 2 . 5 ) 



m. 



),m “ |^(®2/^)(l - bc/ad)^ + + (di^/ 4)(^ - ^)^ + m^(b/a)^j 

- ^ |^-(m2/4)(l - bc/ad)^ + (m^/4)(d/c - b/a)^j 



Note that either (®:j)jjj ^™j^out used for a particular joint in 

Equations ( 5 . 2 . 5 ) as the motion of the system dictates. 

Since the mass effects at the haunches are not separated in 
Equations (5.2.5)^ it is not possible to treat them separately here as was 
done in the derivation of expressions for equivalent soil resistance in 




■ :K 






82 



Section 4.4. This, together with the fact that each of Equations (5.2.5) 
can have either of two possible values, one for inward and one for outward 
movement of the corresponding haunch, makes it necessary at this point to 
consider the possible motions of the model (cf. Section 1.2). There are four 
of these motions as illustrated in Figure 5.3* For convenience, they are 
designated as follows: 

#1 = Asymmetric -Ri^t. This motion is characterized by the 
simultaneous movement of both haunches to the ri^t. 

#2 = Symmetric -In. This motion is characterized by the simultaneous 
inward movement of both haunches. 

^5 = Symmetric -Out. This motion is characterized by the s im\iltaneouB 
outward movement of both haunches. 

y/4 = Asymmetric -Left. This motion is characterized by the simul- 
taneous movement of both haunches to the left. 

Note that in the case of each of the four motions defined above the crown may 
move either upward or downward depending upon the relative velocities and dis- 
placements of the haxinches. 



Ghe notation of Equations (5«2.5) may now be simplified as: 




(5.2.6) 




where 



^l,m = “x (m2/4)(l-bc/ad)^ + mj^(b/a)^ + (m^/4)(d/c-b/a) 
= - (m2/4)(l-bc/ad)^ + (m^/4)(d/c - b/a)^ 



2 





^ = (m2/4)(l-bc/ad) + m^ + ^m^/*tAa/c-D/a; + mg(b/a)^ 

^ = - {m^k){l-hc/eA)^ + (m^/4)(d/c-b/a)^ 




m = index denoting the motion of the system as defined above. 



m = 1,2,5/4. 
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Figure 5*5« Motions of Model 
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CHAPTER VI 

BLAST PRESSURE LOADING 

It is the pxirpose of this chapter to define the loading on the 
dynamic model in terms of the geometry of the arch-soil "structure" and the 
shock wave (hlast overpressure) resulting from the explosion of a nuclear 
weapon. The majority of the "blast parameters are taken from references (l 6 ) 
and (17 )• The shape of the oveiTpressure decay curve behind the shock front 
is taken from reference ( 20 ). Hie assumptions made in selecting these 
parameters are discussed in Section 2 . 4 . 

In defining the load on a "buried structure the shape and char- 
acteristics of the blast wave in soil are of obvious importance. Studies of 
the forces transmitted through the soil cover are being made by several groups 
at the present time. A summary of much of the significant work \diich had been 
done on this problem, together with a convenient numerical technique for solv- 
ing the one -dimensional problem of stress wave transmission throu^ soil, is 
contained in reference (18). From these studies it appears that there is 
little or no significant attenuation in the vertical, free -field overpressures 
within the relatively shallow depths considered in this thesis. Hie results 
of a study of the configuration of the shock front in soil are presented 
herein as Appendix B. 

For derivations in this chapter the blast loeiding is considered to 
be transmitted by tlie soil directly to the arch without modification. Hiis 
is consistent with the approach used in Chapters IV and V in deriving expres- 
sions for the deformation eind inertial resistances of the soil surrounding the 
arch. The assumption is made that the ground surface is flat in the vicinity 
of the arch so that dynamic pressures and reflection phenomena due to surface 
irregularities need not be considered. 
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While the model and the methods of dynamic analysis used in this 
thesis are suitable for studying blast overpressxire loadings due to shock 
waves approaching the structure from any direction, the expressions developed 
in this chapter are derived on the basis of a shock wave approaching from the 
side of the arch, i.e. with the direction of travel of the shock wave normal 
to the longitudinal axis of the arch. This geometry produces maximum eisymmetry 
in the arch loading and thus should lead to maximum lateral response of the 
arch. 

The €LSs\imption that nuclear blast waves may be approximated using 
the idealized ctirves of references (l6), ( 17 ) and (20) is beisic to all of the 
derivations in this chapter. A convenient, but as will be seen not altogether 
essential, assumption is made that the seismic velocity of the soil is practical- 
ly constant between ground zero and the jwint of interest, i.e. the location of 
the structure. 

6.1 Definition of Loading Regions 

The arch itself may receive load at all points on its periphery 
during passage of the shock wave. However, in defining the load on the replace- 
ment model used for the dynamic analysis load is brou^t in only at the joints. 
These joint loadings are simply concentrations of the distributed loading 
sensed by the arch itself. The load concentrations are obtained by defining 
regions in the vicinity of each joint •vdiich will tend to bring load to that 
joint. Then the average blast overpressures acting on these regions are con- 
centrated and applied to the model at the joints. It is more convenient here 
to deal with the horizontal and vertical components of the load at each joint 
than with the resultant forces themselves, just as was the case with the 
inertial forces considered in Chapter V. 
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In the following sections regions tending to bring load to each 
joint, i.e. "loading regions," are defined for various shock wave-soil 
combinations. In general, for shallow-buried structures such as are considered 
in this study, the configuration of the shock front in the soil may be con- 
sidered as a plane perj)endic\ilar to the direction of travel of the shock front. 
It is oriented, in the plane of the arch, at an angle with respect to the 
shock front in air which is dependent upon the relative values of the seismic 
velocity of the soil and the velocity of the shock wave in air at the point 
of interest (cf. Appendix B). Since in this study the abutments are eissumed 
immovable (cf. Section l.d), loading regions tending to bring load primarily 
to the abutments are neglected. 

The three cases considered below are for: 

(1) A vertical shock front in the soil, 

(2) A shock front in the soil which is inclined at an angle 
less than the inclination of the upper bars of the model, and 

(5) A shock front in the soil which is inclined at an angle 
greater than the inclination of the upper bars and less than the inclination 
of the lower bars of the model. 

A series of problems (cf. Chapter VIII and Appendices C and D), 

■vdiich was solved using various angles of inclination for the shock front and 
holding the remaining variables constant, indicated that the lateral response 
of the model was not particxilarly sensitive to variations in this angle. Ih\is, 
the assianption of a constant seismic velocity in the soil between ground zero 
and the location of the structure, while convenient in defining the inclination 
of the shock front in the soil, is not too significant here. Therefore, the 
case of a shock front in soil inclined at an angle between the vertical and 
the inclination of the lower bars of the model is omitted. Also omitted is 
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I the case of a horizontal shock front in soil. While this may be readily 

I 

handled by the techniques developed herein, it is outside the scope of this 
study (cf. Section 2.4). 

I 6.1.1 Case I - Vertical Shock Front in Soil 

From the results of the study presented in Appendix B it can be seen 
that a nearly vertical shock front in the soil results when the velocity of 
the shock front in air at the point of interest and the seismic velocity of 
the soil are nearly equal. Ihis case is illustrated in Figure 6.1. 

Since the configuration of the shock front in soil has been assumed 

I 

I 

as a plane perpendicular to the plane of the arch, it will, in general, strike 
the arch at some point of tangency or at an abutment. The problem of the 
shock wave enveloping the arch itself can be handled quite readily (cf. refer- 
ence ( 21 )). However, for the purposes of this study it is more convenient and 
not significantly less accurate to consider instead the shock wave as enveloping 
the replacement model. 

I 

I 

In interpreting Figure 6.1 consider the shock front in air as moving 
along the ground surface from left to ri^t \diile the shock front in soil, as 
an extension of the shock front in air, moves throu^ the soil. Prior to the 
arrival of the shock front at point Oa', the shock wave is assumed to have no 
influence on the structure. As the shock front moves from Oa' to Ob, the 
shock wave engulfs the "passive" soil failure wedge corresponding to outward 
displacement of the windward haunch, joint 1 (cf. Section 4.2). The average 
overpressvire acting at the ground sxxrface between points Oa' eind Ob is thus 
p^^. As tlie shock front moves from point Oa to point Ob, the shock wave engulfs 
the "active" soil failure wedge corresponding to inward displeicement of the 
windward haxuich, joint 1 (cf. Section 4.l). Hie average overpressvu*e acting 
between points Oa and Ob is thus p^^. 
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As the shock front moves from point 0 to point la, the shock wave 
is assumed to load the left abutment, Joint 0; therefore, this region is not 
considered in the arch loading. The shock front then moves from the midpoint 
of bar 1, point la, to the midpoint of bar 2, point lb, and in so doing is 
assvnned to load joint 1. The average overpressure acting between points la 
and lb is designated p^. 

As the shock front moves on from point lb, the shock wave starts to 
load joint 2; therefore, for convenience, point lb is also labeled point 2a. 

Now the shock wave may be assumed to load joint 2 as the shock front moves 
from point 2a to point 2b. Uie average overpressure acting between points 2a 
and 2b is designated p^. 

It may be noted here that the horizontal component of the loading 
between points 2a and 2 tends to act to the right at joint 2, while the hori- 
zontal component of the loading between points 2 and 2b tends to act to the 
left. For this reason it is convenient to designate the average overpressure 
acting between points 2a and 2 as p^^ and that acting between points 2 and 2b 
as p^^. 

As the shock front moves on from j)oint 2b, it starts to load joint 5. 
Accordingly, it is convenient to label point 2b also as point 3a. Now, the 
average overpressure acting between points 3a and is assumed to load, joint 3 
and is designated p^. 

The average overpressure acting between points 3b and 4 is assumed 
to load the ri^t abutment, joint 4, and is therefore neglected. 

As the shock front moves from point 4a to point 4b, the shock wave 
engulfs the "eictive" soil failure wedge corresponding to the Inward displace- 
ment of the leeward ha\inch, joint 3* aversige overpressure acting between 

points 4a and 4b is thus Pg^* 
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As the shock front moves from point to point 4b the shock vave 
engulfs the "passive" soil failure wedge conresponding to the outward dis- 
placement of the leeward haunch, joint 5* Ihe average overpressure acting 

between points 4a and 4b' is thus p’,. 

s:5 

For purposes of calciolating the "rib-shortening" correction derived 
in Section 5*5^ the average overpressxare acting on the entire arch is needed. 

This is taken as the average overpressirre acting between points 0 and 4 and 
is designated p . For convenience in this case point 0 is also labeled "aa" 

8l 

and point 4 is also labeled "ab." 

Note in the foregoing description that no actual "loading" of the 
arch takes place until the shock front reaches point la. It is assumed that 
no movement of the arch occurs prior to this time. Therefore, the time of 
arrival of the shock front at point la is taken as zero time for this case. 

This assvimption is consistent with the derivation of the "deformation" resistance 
of the soil in Chapter IV. It is not unreasonable ■vrtien one considers the hi^ 
velocity of the shock wave, the relatively small dimensions of the arch-soil 
"structure, " the inertia of the soil loaded prior to this time and the 
relatively small displacements of the structure during the early part of the 
loading history. 

In the derivation of expressions for the loeid concentrations on 
the model (cf. Section 6.2) it is convenient to use the times at idiich the 
shock front arrives at the windward and leeward bomdaries of the various 
loading regions. Zero time for this ctise has already been chosen at point la. 
Hiis is expressed algebraically by Equation (6.1.1). 

t^^ = 0 (6.1.1) 

The remaining time parameters are given by Equations (6.1.2), expressed in 
terms of the velocity of the shock front in air, U, and the geometry of the 
arch-soil "structure" as shown in Figure 6.1. 
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^Oa' 




- [l - 1 


j)/2U 


"'^Oa 


= 


Fn V 

- (b + 2d^cot |^5- + 2 


)/2U 

■ 


"'^Ob 


= 


- (b - 2bdya)/2U 




^Ib 


= 


(b + d)/2U 




^2a 


= 


‘'^Ib 




^2 


= 


(b + 2d)/2U 




'*^2b 


= 


(b + 5d)/2U 








Hb 




V 


= 


(b + 2d)/U 




Ha 




(-b + 2B - 2bd^a)/2U 


^4b 


= 


(-b + 2B + 2d^cot 1 


; . |j)/» 
; ■ ljv« 


Hb' 


= 


(-b + 2B + 2d^cot 1 


^aa 




- b/2U 




*ab 




(-b + 2B)/2U 





( 6 . 1 . 2 ) 



where ^ in general: 




= time of arrival of the shock front at the windweird 
boundary of loading region j 



t ., = time of arrival of the shock front at the leeward 
boundary of loading region J. 



6.1.2 Case II - Shock Front in Soil Inclined at Angle Less Than 
Inclination of Upper Bars of Model 

The results of the study presented in Appendix B indicate that for 
seismic velocities less than the velocity of the shock wave in air at the 
point of interest, the shock front in soil will be inclined at an angle to 
the shock wave in air. The expressions derived in this section are for the 
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case vhere the inclination of the shock wave in soil with respect to the 
horizontal, is less than the angle of inclination of the upper bars of 
the model, i.e. : 

T] < arctan (6.1.3 ) 

This case is illustrated in Figure 6.2. Note that in this case 
the shock front in soil first strikes the model at the crown. Therefore, for 
this case the time of arrival of the shock wave at the crown is taken as zero 
time. 

tg = 0 (6.1.1^) 

The considerations in selecting this time as the reference are the same as 
for the preceding case (cf. Section 6.1.1). 

The loading regions shown in Figure 6.2 are defined in an analogous 
manner to those for the verticeil shock front in soil and therefore are not 
described in detail here. However, there are several differences worthy of 
note, viz: 

1. The "active" and "passive" soil failure wedges are not fully 
enveloped by the shock wave until the shock front in soil reaches the corre- 
sponding abutment. 

2. On the windward side of the model the shock front strikes the 
midpoint of the upper bar before reaching the midpoint of the lower bar. 

However, the midpoints of the bars are not relabeled. This introduces an 
inconsistency in the notation, but not a serious one, since the order in \diich 
the boundaries of a "loading region" are reached by the shock front is obvious 
in any set of computations. 

5. The crown, joint 2, is enveloped by the shock wave from the top 
downward instead of from left to right as was the case in the preceding section. 
Here also the notation is not changed. 
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The average overpressure acting on the eorch for this case is 
taken as the average overpressure acting between points 2 and 4. This 
corresponds to the envelopement of the model by the shock wave as the shock 
front travels from the crown to the leeward abutment. For convenience, 
points 2 and 4 6ure also labeled "aa" and "ab." 

The time parameters defining liie times of arrival of the shock 
front at the boundaries of the loading regions are given in Equations (6.1.5)* 
These may be written readily by referring to Figure 6.2. 



^Oa' 



"Oa 



‘'Ob 



^la 



'lb 



'2a 



'2b 



5a 



5b 



'4a 



'4b 



'4b' 



aa 



ab 



- d /C - (B + 2d cot 

c s o 

- d /C' - (B + 2d cot 

c s o 




(a + c)/C - B/2U 
s 

(a + 2c)/2C' - (b + 2d)/2U 
s 

c/2C' - d/2U 
8 

^Ib 

c/2C^ + d/2U 
^2b 

(a + 2c)/2C' + (b + 2d)/2U 
s 

- d^/c^ + (B - 2bdya)/2U 

(a + c)/C* + b/2U 
6 




(6.1.5) 
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Equations (6.1. 5 ) are written in terras of the apparent vertical 

component of the velocity of the shock front in soil, C'. If the inclination 

s 

of the shock front in soil, q, with respect to the horizontal is obtained from 
Appendix B, Equation (B.J), C may be obtained from Equation (6.1.6). 

C = U tan q (6.1.6) 

5 

6 . 1.5 Case III - Shock Front in Soil Inclined at Angle Intermediate 
Between Inclinations of Upper and Lower Bars of Model 

The expressions derived in this section are for the case where the 

inclination of the shock front in the soil is greater than the inclination 

of the upper bars of the model and less than the inclination of the lower 

bars of the model, i.e. : 

arctan j < q < arctan ^ (6.1. 7 ) 

This is illustrated in Figure 6.5* Note that the shock front strikes the 
windward haunch first. Therefore, the time of arrival of the shock front at 
the windward haunch is taken as zero time for this case. 

t^ = 0 (6.1.8) 

The considerations for the selection of t^ as the reference time are the 
same as those listed in Section 6.1.1 for the selection of the reference time 
for the vertical shock front in soil. 

Ihe loading regions shown in Figure 6.5 are defined in the same 
fashion as those in the preceding sections. However, several differences do 
exist for this case, viz: 

1. The windward haunch is assumed to be loaded as the shock front 
moves from the haxmch to the midixjint of either the upper or lower adjoining 
bar depending upon which occvurs later. In other words, the shock front may 
arrive at point la either before or after it strikes point lb. The order of 
arrival depends upon the relative geometries of the shock front and the model. 
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2. Ihe soil failure wedges on the leeward side of the model are 
assumed to be loaded as the shock front moves from point 4a to either the 
leeward abutment or the intersection of the failure plane and the ground 
surface, vdiichever occurs later. Uie order of arrival depends upon the 
relative geometries of the shock front and the model Just as in the case of 
the windward haxinch discussed in the preceding paragraph. 

5. The average overpressiure acting on the arch in this case is taken 
as the average overpressure acting between points 1 and 4. For convenience 
these points are also labeled for this case as "aa" euid ”ab". This corresponds 
to the envelopment of the model by the shock wave £is the shock front travels 
from the windward haunch to the leeward abutment. 

Ihe time parameters defining the times of arrival of the shock front 
at the bovmdarles of the loading regions are given by Equations (6.1.9)* These 
may be written readily by referring to Figure 6.5* 



^Oa’ 



Oa 



^Ob 



*t 



*t 



la 



lb 



2a 



- (c + d^)/c; - (b + d^cot - 1 )/u 

- (c + d^)/Cg - (b + d^cot + I )/u 

a/C^ - b/u 

a/2C' - b/2U 



- c/2C^ + d/2U 



= t 



2b 






3b 



lb 

c/C + d/U 



- c/2C^ + 3d/ 2U 



= t 



2b 

a/2C + (4d + b)/2U 
8 



(6.1.9) 



* Use either or 

la Id 



whichever is larger. 
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+ (-b + B - bdya)/u 
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0 I 
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+ (-b + B + d cot 

0 1 


r» 

[f-sj 



)/U 
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Equations (6.1.9) are written in terms of the apparent vertical 

component of the velocity of the shock front in soil^ C\ If the inclination 

of the shock front in soil, with respect to the horizontal is obtained from 

Appendix B, Equation (B.jj), then C may be obtained from Equation (6.1.6). 

6 



6.2 Load Concentrations on Model 

The load concentrations on the model are derived in this section by 
considering the average overpressure acting at the ground surface on the load- 
ing regions defined in Section 6.1. Bie approach used is illustrated in 
Figure 6.4. The overpressure is eissumed to have no effect upon the region 

until time t . , the time of arrival of the shock front at the windward boundary, 
ja 

Ja, of the region, J. The average overpressure acting on the region, p ., is 

J 

then assumed to build up linearly as the shock front traverses the region from 
the windward boundary, Ja, to the leeward bovindary, jb. At times greater than 
the time of arrival of the shock front at the leeward boxmdary, tj^, the 



Use either (^ 4 ^)^ ^^4b^2^ whichever is larger. 

Use either (tj^^t )^^ or 'diichever is larger. 
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I average overpressiire acting on the region is eissumed to have a value equal to 
the instantaneous value of the overpressvtre at the center of the region. This 
approach is essentially the same as is used in reference (l6) in defining the 
loading of flat surfaces oriented parallel to the direction of travel of the 
shock wave. Various techniques involving the integration of the oveirpressure- 
time pvilse followed by an averaging of the Impulse over the region of interest 
were studied and found to offer little or no euivantage over the simpler approach 
described above. 

In general, the average overpressure acting on a loading region is 
expressed by Equations (6.2.1) 



t < t . : 
- Ja 



t, < t < t,, : 
Ja - - jb 



Pj(t) = 0 

't - t 

Pj(t) = p (-- i l 



L^Jb ''ja- 






t, +t 



tjb < t < I Pj(t)= P (t - (6.2.1) 



^ (t. +t., )+t, < t 
2 ' ja jb' d - 



: Pj(t)= 0 



where 



p.(t) = average overpressure acting on region J at any time, t 
^ (cf. Figure 6.4) 

p(t) = variation of overpressure behind the shock front with 
time (cf. reference (20)) 



t = time 

t^ = time of arrival of the shock front at the windward 
boundary of region j 

= time of arrival of the shock front at the leeward 
boundary of region, j. 



Ja 



'Jt 



It has been pointed out earlier, in the cases of the inclined shock fronts, 
that because the reference points on the model are labeled from left to ri^t, 
inconsistencies in the notation exist when the shock front approaches a region 
from a different direction. This occurs, for example, in Case II where the 
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shock front strikes point lb before arriving at point la (cf. Figure 6.2). 
This type of situation meiy be handled readily using Equations (6.2.1) by con- 



the designation of the corresponding points on the model. 

Expressions for the concentrated loads acting at the joints of the 

model may now be written in terms of the average overpressures acting on the 

corresponding loading regions, Equations (6.2.1). Here the assumption is 

used that the overpressure acting at the grovuid surface is transmitted vertical- 

ly downward by the soil ais a vertical overpressure in the soil. ' ' Further, 

it is assumed that the horizontal oveirpressure in the soil, p^, bears a con- 

(l5 17 19) 

stant relationship, to the vertical overpressure, p^. ’ 



Note that the following expressions. Equations (6.2.3), are for the horizontal 
and vertical components of the load concentrations at the joints rather than 
for the resultant loads themselves. Hiis is simply for convenience in deriva- 
tion and for use later in determining the equivalent loads on the system 
(cf. Section 6 . 3 ). Referring to Figures 3*2 and 6.1, 6.2 or 6.3^ one may 
write: 





(6.2.2) 



“1 = - 1 

*2 “ ■ I 



2 pr'^2a ^2b 




( 6 . 2 . 3 ) 



= + I P^(b+d) 




^6 = + 2 



Note that since the oveiT)ressure acts downward, it is considered negative in 
accordance with the sign convention used herein (cf. Section 1.2). 
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6.5 Development of Equivalent Load 

The concentrated loads acting at the joints of the model mvist be 
combined in order to obtain their resultant effect upon the tvo -degree -of - 
freedom system. The "Load Relationships" derived in Section 5 . 2 . 1 ^ 

Equations ( 5 . 2 . 5 )^ express the resultant loads in terms of the equivalent 
horizontal loads at the haunches. Ihls is in agreement vith the choice of 
the horizontal displacement components of the haunches as the independent 
variables for the dynamic analysis. 

While Equations ( 6 . 2 . 5 ) and (5.2.5) may be used successively to 
obtain the equivalent loud on the system, it is more convenient to combine 
them and express the equivalent loads in terms of the average overpressvtres 
acting on the various loading regions. Equations (6.5.1) are obtained by sub- 
stituting Equations (6.2.5) into Equations (5.2.5)* 

^1 = ' 2 ^1 - I 



+ I P2d(d/c - b/a) 



(6.5.1) 



- 2 -b/a) + ~ Pjj |Kpj.(a+c) + (b/a)(b+d) 

The notation in the preceding equations may be simplified as shown 
in Equations (6.5.2). 

®i ' - •'A - ¥^ 2 ^ 

q = I Kp^(atc) + I (b/a)(b-Mi) 

■'2*5 

Jj . I a(d/o-b/a) 



- Jj(P2^ - P2t> - JjPj ■. qp. 



(6.5.2) 



where 
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Note that the dimensions of J^, and must be in agreement vith the 
units in "vdiich the overpressure is expressed. 



I 
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Figure 6.2. Loading Regions for Case II 
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Figure 6.^. Generalized Overpressure -Time Diagram 
for IVpical Loading Region 
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CHAPTER VII 
EQUATIONS OF MOTION 

Now that expressions for the equivalent forces, masses and re- 
sistances of the two -degree -of -freedom system have been derived, it is possible 
to develop the equations of motion for the system in terms of these quantities. 
While a simplified notation is used throxighout this chapter, the expressions 
may be expanded readily by inserting the proper expressions from Chapters III 
through VI. As an aid in visualizing the two -degree -of -freedom system, the 
equivalent masses, resistances and loads are considered to exist and act at 
the haunches, i.e. at joints 1 and 5 of the model. This is only a physical 
analogy, but it is convenient in comparing the expressions derived in this 
chapter with those of the preceding chapters. 



7.1 Development of Equations 

The forces acting on one of the two equivalent masses are shown in 
Figure 7*1 in general terms. A general expression for the equations of motion 
may be obtained by applying d'Alembert's Principle and equating the sum of the 
horizontal forces acting on the mass to zero. The result is given in 
Equation (7.I.1). 



where 



m, h . + Q. + Q . = H . 
J J sj J 



m. = equivalent mass of j 
J^ni 

h. = acceleration of mass j 
J 

Q. = equiveilent deformation resistance 
J 

Q^j = equivalent deformation resistance 

H. = equivalent load etc ting on mass j. 
J 



of the arch at J 
of the soil at j 



(7.1.1) 



Since each of the equivalent masses is actually a function of both 
of the independent accelerations (cf. Equations (5.2.6)), the two equations 
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of motion as expressed by Equations (7*l*l) are not in a convenient form for 
integration. However, since Equations ( 5 . 2 . 6 ) are linear in the acceleration 
terms, the equations of motion may be rewritten explicitly for the independent 
acceleration components as shown in the following derivation. 

Equations (7.1.2) are obtained by substituting Equations (5.2.6) 
into Equations ( 7 .I.I). 



4, h, - V h, = H- - (i, - Q - 
l,m 1 l,m 5 1 u si 

-V h^ + 4, h, = H, - Q- - Q , 
5,m 1 3,m 3 3 ^ s5 



( 7 . 1 . 2 ) 



Solving Equations ( 7 . 1.2) simultaneously for h^^ and h^ gives Equations (7.1.3)< 



h, = 



JLiR. 






1 4, 4, -V V 

l,m 3>® 1^® 3^ni 






l,m^,m l,m l,m ^,m l,m 3^™ 



l,m 3^m l,m ^,m 



l.m 



( 7 . 1 . 5 ) 



The notation of Equations ( 7 . 1.3) niay be simplified as shown in Equations 

( 7 . 1 .^). 



(T.l.it) 



vhere 






l,m 4 , 4, -V. V 

' l,m ^,m l,m ^,m 



ItP 



l,m 4, 4, -V, V 

l,m ^,m l,m 






l,m 3>® l,m PfXa 



l,m 



3,m 4, 4., -V V 

^ l,m l,m 3>® 



m = index denoting the motion of the system as defined in 
Section 5.2. m = 1,2,3>^* 
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Equations (7.1.4) may now be integrated slnultaneously in order to 
obtain the velocities and displacements of the system at any time under a 
given set of dynamic loads. 

7.2 Natural Periods of Vibration 

The natural periods of vibration of any system are highly signifi- 
cant parameters in predicting the response of that system to dynamic loads. 

The natiural periods of the model which is used in this study for the dynamic 
analysis may be found quite readily from the equations of motion derived in 
the preceding section. 

Before proceeding with the derivation of expressions for the natural 
periods of vibration of the model, however, it is convenient to express the 
resistance of the arch. Equations (3.4.7)> in a slightly different form. 
Equations (5.4.7) express the equivalent resistance of the arch in terms of 
the resisting moments developed at the hinges of the model. While this is con- 
venient for use in an elas to -plastic anedysis, it is inconvenient for the 
purpose of finding the natural periods of vibration corresponding to the 
initial, elastic response of the system. 

Hie elastic resistance of the model may be expressed in terms of 
the independent displacement coordinates by substituting Equations ( 3 . 2 .II) 
and ( 5 . 3 . 1 ) into Equations (5.4.7). 

( 7 . 2 . 1 ) 

S = K3(-A^h^+A^h^)(q^+q^h3^-q^h^) + Kg(-A2h^+A2h^)(q^+q^h3^-q^h^) 

Equations ( 7 . 2 .I) may be simplified by limiting consideration to small dis- 
placements such that the products of the displacements eire small with respect 
to the remaining terras and may be neglected. The simplified expressions are 
given as Equations ( 7 . 2 . 2 ). 



m 



% “ - (^VVc Vl *■ (Va*McV} 
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( 7 . 2 .^) 



The notation in Equations ( 7 . 2 . 2 ) may be simplified as shovn in Eqxiations 

(7.2.5). 



S = ^ 1^1 - '^ 5^5 

% = - 

where = (A^q^ 

= (Kg/2L^)(5+a/c+b/d+2c/a+bc/ad) 

= (Kg/ 2 L ^)(3 + a/c - b/d - bc/ad) 



(7.2.3) 



The discussion will now proceed with the derivation of expressions 
for the natiiral periods of vibration of the system. One may write the equations 
of motion for the system in free vibration. Equations (7.2.4), directly from 
the general equations of motion for the system. Equations ( 7 .I.I), by elim- 
inating the loading terms. 

m, Ji. + Q + Q = 0 (7.2.4) 

J,m j J sj 

If a sinusoidally varying time function is assumed for the motions of the 

system, the accelerations, ti., and displacements, h., of the system are related 

J J 

as shown in Equations ( 7 . 2 . 5 ). 

where w = circuleur natural frequency of the system for the particxilar 
motion being studied (cf. Section 5.2). 

Equations (7.2.4) may be rewritten in expanded form by inserting Equations 

( 5 . 2 . 6 ) and ( 7 . 2 . 5 ). The result is given as Equations (7.2.6). 
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•M, h + V h, + 0, + Q T = 0 

l,m m 1 l,m m3 ^ si 

2 2 — 

V w h, - 4, w h, + Q, + Q , = 0 

5, mm 1 m 3 j s3 



( 7 . 2 . 6 ) 



Substituting Equations (4.4.14) eind ( 7 . 2 . 3 ) into Equations (7.2.6) gives: 

^1 ■ "l,m“m “3 ■ ‘l"l ‘3“3 ■ '%l'^'“l = ^ 

( 7 . 2 . 7 ) 

.2 . .2 



u (jT h- - w, w h, - T h, + t,h, - (K T/a)h, = 0 
l,m m 1 l,m m3 11 3 3 ' sr ' 1 

00 — . . 

-V w h, + u, w h, + T h, - T h, - (K ,/a)h, = 0 
3>m m 1 3^m m 3 3 1 13 s3 3 



Note that Equation (4.4.15) must be used in lieu of Equation (4.4.14) for 
either of Equations ( 7 . 2 . 7 ) when the motion being considered is such that the 
haunch is displaced outward against the soil rather than inward into the arch. 
This qualification also applies to the remainder of this derivation. Rewriting 
and rearranging Equations ( 7 . 2 . 7 ) gives: 






( 7 . 2 . 8 ) 



Equations (7.2.8) yield solutions different from zero only if their 
determinant vanishes. This permits Equation (7.2.9)^ the frequency equation, 
to be written in terms of the circular natural frequencies of the system. 





- v_. 


-1 


L l,m 3,m 


l,m 


■" 1 






. [(- - <’ 5 >^ 



m 

1 



= 0 



( 7 . 2 . 9 ) 



TSie circular natural frequency for any one of the four motions of 

the system (cf. Section 5 . 2 ) may be obtained by substituting the proper values 

of u, , V, , , and v_ (cf. Equation (5.2.6)) and the proper expressions 

X in ^ ^ HI ^ ^ m 

for soil stiffness, i.e. either K . or K* . (cf. Equations (4.4.14) and (4.4.15)), 

s J S J 
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into Equation (7.2.9) and solving for w . Note that since both K . and K' . 

m sj sj 

are functions of the overpressure restraining the soil, the circular natural 
frequency, hence the natural period of vibration, of the system will vary as 
the shock wave passes over the structure. In general, solution of Equation 
(7.2.9) yields two values of Ihe proper value may be selected by sub- 

stituting the values, in tvim, into Eqviations (7.2.8) and noting the relation- 
ship between the displacements h^^ and h^. Kie proper value is indicated >dien 
the relationship between h^ and h^ agrees with the system motion being 
studied. 



Once the proper circular natviral frequency has been found, it may 
be substituted into Equation (7. 2.10) and the natural period of vibration of 
the system, T^, obtained for the particular system motion eind overpressure 
conditions being studied. 



_ 

T = — 
m u 
m 



(7.2.10) 



While there is no real adveuitage in solving Equation (7.2.9) in 
general terms, a simplified expression for the circular natural frequency can 
be written for those special cases in which the stiffness and mass parameters 
are symmetrical. In this study this condition exists, for example, for the 
"Symmetric -In" and "Symmetric -Out" motions with no overpressure restraining 
movement of the soil. Equation (7.2.II) is obtained by setting M- = = 

JL^ m ^ jjn 



m' 



= 



l,m 



= and in Equation (7.2.9) and solving for the 



m 
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value of corresponding to symmetrical response. 



u = 
m 



+ T * Kjl 



M + V 

m m 



(7-2.ll) 



Note that K'/a may be used in place of K /a in Equation (7.2.II) if the motion 
s s 



being studied so indicates. 



I'l 
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Figure 7-1* Forces Acting on Typical Equivalent Mass 
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CHAPTER VIII 
NUMERICAL STUDIES 

This chapter covers the basis of the numerical studies made in 
connection with this thesis. The range of variables Included in the studies 
is given together with a description of the blast wave properties used. The 
numerical examples are described and, finally, a set of sample cedcvilations 
is included. The resxilts of these studies are presented in tabular and in 
graphical form in Appendices C and D. 

8.1 Range of Variables Studied 

The quantitative investigations made in connection with this thesis 
eu:e designed to Illustrate the behavior of shallow -buried, two -hinged arches 
subjected to nuclear blast loading. Two different arches are vised, one a 
semicircular, reinforced concrete axch, l 6 ' 8 " in diameter, eind the other a 
coriMgated metal arch in the form of a segment of a circle with a rise of 
11 ' 8 1/2" and a span of 25 ' 1 l/ 4 ". The concrete eo'ch hsis a uniform thick- 
ness of 8 " and is reinforced with No. k bars running circumferentially on 
10" centers at mid-thickness. Intermediate grade reinforcing steel is used 
with an assumed yield strength of 50 ksi under dynamic load. The elastic 
modulus of the steel is taken as ^ 0,000 ksi. The average strength of the 
concrete eis measured by the standard 6 " x 12 " test cylinder is taken as 
3.5 ksi at 28 days. Substituting this value for the concrete cylinder 
strength into the appropriate expressions of reference (22) yields a corre- 
sponding elastic modulus for the concrete of 3^820 ksi. The cormgated 
metal eurch is 10 gauge with an elastic modulus of 30,000 ksi and an assumed 
yield strength of 30 ksi under dynamic loading. 
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The arches are studied at depths of cover over the crown veirylng 
from zero to one -half of the arch span. Intermediate depths of cover of 
one-ei^th and one -quarter of the span are used. 

The soil cover is taken as a granular, cohesionless material with 
a vmit weigit of 120 pcf, an angle of interneil friction of 50 ° and a hydro- 
static pressure ratio of horizontal to vertical soil pressures in the "at 
rest" state of stress of 0.5. The friction angle between the soil and the 
arch is assumed as two-thirds of the angle of internal friction of the soil, 
or 20°. The seismic velocity of the soil is assumed as 3>000 feet per second 
for the majority of the numerical work; however, a series of problems is 
Included varying the seismic velocity from 1,000 feet per second to '^,000 
feet per second in order to study the significance of this parameter. 

The blast loading is assxnned to result from a contact sxirface 
burst of a one -megaton nuclear weapon. Overpressures ranging from 50 psi to 
1,000 psi are considered. The shock wave properties are discussed in 
Section 8.2. 

Approximate calculations based upon reference (l6) indicate that 
the overpressures selected occur outside of the crater and therefore are in 
agreement with the loading assumptions discussed in Section 2.4. Structures 
located at the higher overpressures will probably be inside the fireball and 
will also be subjected to heavy doses of nuclear radiation. Again, as dis- 
cussed in Section 2.4, these effects must be considered in any actual design; 
however, they will have little or no effect upon the structurcil response of 
the arch. 

The majority of the numerical solutions are computed assxanlng a 
ratio of horizontal to vertical overpressure in the soil of one-fovurth. In 
addition a series of problems is included in which this ratio is varied from 



one -fourth to one. 
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8.2 Shock Wave Properties 

The shock vave parameters used in this study are taken from ref- 
erences (l6), (17) and (20). A summary of the more significant parameters 



on the duration of the positive phase of the shock wave in the higher over- 
pressure regions. The differences in the values given are not great, however, 
when one considers the actual, variation in all shock wave parameters to be 
expected in the field. Further, small differences in positive phase duration 
are not significant in this study because of the rapid response of the arch 
(cf. Appendix D). However, the rapid response of the arch does make the 
system quite sensitive to the initial shape of the blast wave. Therefore, 
while the positive phase duration for the 1000-psi shock wave may be taken 
as the same as that for a 100-psi shock wave without significantly affecting 
the results, the initial shape of the two shock waves must be defined 
carefully. This latter aspect of the problem is discussed in detail in the 
following paragraphs. 

The results of a study of the configia*ation of the shock front in 
soil are presented in Appendix B. 

The decay with time of the overpressxare behind the shock front is 
discussed in reference (20). Later studies of shock wave characteristics 
are available but are not used here because of their security classification. 
Reference (20) gives the following approximate expressions for the vsiriation 
with time of the overpressvire at a point behind the shock front during the 
passage of the positive phase of the shock wave. 




0 1 "t' < p(t')/PgQ = (1 - ■t’/t^) e“^^^ ^^d^ 



(8.2.1) 



45 psi < p^^ < JOO psi: G = 



■ ^ l+c^(t'/t^) 



(8.2.2) 



where 



t' = elapsed time at the reference point after the shock front 
has passed 

p(t' ) = overpressure at the reference point at time t' 
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p = peak value of the overpressure as the shock front passes 
the reference point 



t^ = duration of the positive phase of the shock vave at the 
reference point 

e = natural logarithmic base - 2.71828... 

^ = - 0.231 + 0.588(pg^/p^) - 0.0352(PgyP^)^; p^^ < I50 psi 
= 0; PgQ > 150 psi 

= (p 7 p^)[0.88 + 0.072(p /P^)]; p_ < 150 psi 



SO 



so 



so 



L- (P3„/P„)IX.67 - 0.01l(p^^/p^)J; > 150 psi 

c^ = 8.71 + 0.l81i3(P3yPo) - 104/[(pg^/p^) + 10] 

P = ambient atmospheric pressure (taken as 14.7 psi for 
numerical vork in this study). 

The overpressure-time relationship for 50 psi is obtained by sub- 
stituting p =50 psi eind t, = I.60 seconds (cf. Table 8.I) into Equations 
(8.2.1) and (8.2.2). The resulting expression is given as Equation (8.2.5 )• 
p(f)/50 = (l-f/l60) e-[°-T05+5.85/(lif))tV'l.60 (8.2.J) 

The computer program used in this study for the solution of the numerical 
examples (cf. Section 8.3) approximates the "actual” overpressvire-time 
cuarves with a series of straight-line segments. These line segments may be 
obtained for the 50"Psi shock wave by plotting Equation (8.2.5) and selecting 
a series of strai^t-line segments which approximate the curve to a reason- 
able degree over the time interval of interest. The coordinates of the end- 
points of the lines selected eire given in Table 8.2. 

The overpressure -time relationship for 100 psi may be obtained by 
substituting p = 100 psi and t, = 1.48 seconds (cf. Table 8.I) into 

SO cl 

Equations (8.2.1) and (8.2.2). The resulting expression is given as 
Equation (8.2.4). 

p(t’)/l00 = (l-t’/l.48) l+3.78(t7l.48)^ ITTO ^ 



Proceeding as in the case of the 50"Psl shock wave gives the lines defined 
by the overpressure-time coordinates listed in Table 8.3 as an approximation 
to the overpressure -time curve. 
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Since the 1,000-psi shock wave falls outside the range in which 
Equations (8.2.1) eind (8.2.2) apply, the overpressure -time relationship for 
this shock wave may be obtained by Interpolating in the curves of references 
(20) and (25). This was done euid a series of strai^t-llne segments selected 
which approximate the "actual" overpressure -time curve reasonably well. Hie 
line segnents selected are defined by the overpressure -time coordinates 
listed in Table 8.4. 

8.^ Description of Numerical Examples 

The properties of the arches selected for the numerical examples 
are given in Section 8.1. Table 8.5 indicates the variations in parameters 
applied to each of these arches. 

The lower vaJ-ue of overpressure to which the reinforced concrete 
arch is subjected is well below the level which should cause buckling or 
general yielding of the arch due to thrust in the arch rib (cf. Section 8.4). 
The hi^er value of 1,000 psi is selected to magnify any tendencies of the 
arch to respond laterally. It is not intended to be below buckling or 
general yielding limitations. As pointed out in Chapter IX, however, the 
arch rib contributes relatively little to the lateral resistance of the 
arch-soil "structure" as compared with the soil contribution. Therefore, 
the 1,000-psl results are indicative of the response of an arch with a 
sufficiently large cross section that buckling and general yielding are no 
longer critical. 

The 50~psl loading on the corioigated metal arch is reasonably 
close to the overpressure causing general yielding due to thmst in the arch 
rib for the assumed soil and loading conditions. In fact, the case of zero 
depth of cover over the crown could very well be ruled out of any practiceil 
consideration on this basis alone; however, it is included in this study for 
comparison purposes. The 100-psi and 1,000-psi loadings eu*e used to give a 
comparison with the corresponding loadings on the reinforced concrete arch. 
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Further, the observation made above In the case of the reinforced concrete 
arch, regarding the relative contributions of the arch and the soil to the 
lateral resistance of the system, is even more applicable here. Therefore, 
it may be stated that if the arch is strengthened sufficiently to withstand 
the thrusts developed in the arch rib, the resjxjnse resvilts of the lOO-psi 
and 1,000-psi studies should be indicative of the response of the strengthened 
sti*ucture . 

Integration of the equations of motion for each of the nvtmerical 
examples was done on the Illiac, the hi^-speed digital computer at the 
University of Illinois. Structural Research Laboratory Code 462^^^^ was 
adapted to perform this integration. 

Tabular comparisons of the results obtained from the examples 
listed in Table 8.5 are presented in Appendix C. Response curves and con- 
figuration diagrams for the arch at various critical times are presented 
in Appendix D. 

8. 4 Sample Calcxilatlons 

Calculations are presented here for the reinforced concrete arch 
with a depth of cover over the crown eqvial to one -fourth of the arch span. 

The arch and soil cover are shown in Figure 8.1. The arch is subjected to a 
lOO-psi shock wave from the explosion of a one megaton nuclear weapon as a 
contact svirface burst. The seismic velocity of the soil is assumed as 5,000 
feet per second and the horizontal overpressure induced in the soil by the 
shock wave is assumed to be one-fourth of the vertical overpress\ire. In 
general, the calculations are based on a one -foot width parallel to the longi- 
tudinal axis of the arch (cf. Section 2.l). 
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8.4.1 Calculation of Peirameters eind Coefficients 

The horizontal displacements of the haunches, h^^ and h^, are 
chosen as the Independent displacements. The remaining displacements are 
related to them by Equations (5.2.7)* 



h^ = I (l-bc/ad)(h^+h^) 

= I (l-(2. 440)^/(5.895 )^](h^+ h^] 
h^ = 0.4l45(h^+ h^) 



(8.4.1) 



^4 = 



(b/a)h^ 

= - (2. 440/5.895 )h^ 



h|^ = - 0.4l4l hj^ 



(8.4.2) 



II5 = I (d/c-b/a)(h^-h^) 

= I (5.893/2.440 - 2 . 440/5.895 )(h^-hj) 

h^ = 1.0006(h^-h^) 
hg = (b/a)h^ 

= (2. 440/5.895 )h. 



(8.4.5) 



hg = 0.4141 h^ 



(8.4.4) 



The total rotations of the joints are related to the independent 
displacements by Equations (5.2.II). 



= I (2/a+l/c+b/ad) 



= I [2/5.893 + 1/2.440 + 2.440/(5.893)^] 
= 0.4098 ft“^ 

Ag = l/c 

= 1/2.440 
= 0.4098 ft"^ 
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= I (l/c-b/ad) 

= I [1/2.440 - 2.440/(5.893)^] 

= 0.1698 ft"^ 



'*'1 " 

= 0.4098 - 0.1698 



(8.4.5) 




tg = - 0.4098 + 0.4098 hj 



(8.4.6) 



'*'5 “ 'S^l"'^^3 

\|r^ = 0.1698 - 0.4098 



(8.4.7) 



The genereil form of the moment-rotation relationship for the model 



hinges is shovn in Figure 8.2. The dimensions on this figure may be found 
using Equation (3.3. !)• Ibe material properties for the arch are listed in 
Section 8.1. Figure 8.3 is a sketch of the cross section of the arch rib. 
Nov, in the derivation which follows, neglect direct stress in the arch rib 
and let: 



h = thickness of the concrete €irch rib 
a 

d^ = effective depth of the concrete arch rib 

b = width of the eirch rib parallel to the longitudinal 
axis of the arch 

J = ratio of the distance, jd , between the resultants of the 
compressive and tensile stresses to the effective depth, d 

& 

k = ratio of the distance, kd , between the extreme fiber 
and the neutral axis to tte effective depth, d^ 

Ig = gross mement of inertia of the eirch rib cross section 

= transformed moment of inertia of the arch rib cross section 
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I. = average moment of inertia of the arch rib cross section 
J 

Ag = area of reinforcing steel in a one -foot vidth of arch 

f* = 28-day strength of a standard 6" x 12" concrete test 
cylinder 

f , = yield strength of the reinforcing steel under dynamic 

loading 

= modulus of elasticity of the reinforcing steel 
Ej = modulus of elasticity of the concrete in the arch 
n = modular ratio for steel and concrete. 

From reference (22); 

n = 5 + 10, 000/ f^ 

= 5 + 10,000/5,500 
= 7.86 



Bien: E. = E /n 

J s' 

= 50 X 10^/7.86 

= 5 *82 X 10 ^ psi 

M = A f , Jd 
y s dy*^ a 

= 0.24(50)(7/8)(4)(i/12) 

=5*50 ft -kips 

I = b (h 9/12 
g a' a' ' 

= 12(8)^12 

= 512 in 

It = + nAgd^(l-k)2 

= y [(5/8)(4)]^ + 7.86(0.24)(U)^(1-5/8)^ 

4 

= 25.5 in^ 



% 
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= I (512 + 25.5) 
= 268.6 in^ 



= (5.82 X 10^)(268.6)/(6.578)(12) 



= 1117.2 ft -kips/ radian 



♦y = 

= 5.50/1117.2 

= 0.005155 radians 

Tlie relationships shown in Figure 8.2 are expressed algebraically in 
Eqxiations ( 8 . 4 . 8 ) 



0 < < .iTy.. Mj = Kg 








( 8 . 4 . 8 ) 



Note that the value of t must be adjusted for any "permanent set" in the 

y 

hinge as the structure responds. Also, Equations ( 8 . 4 , 8 ) are valid for 
rotation in either direction. For convenience, the inequalities at the 
left of Equations ( 8 . 4 . 8 ) are shown for the positive direction of rotation 



only. 



Ihe equiveilent arch resistance may be found using Equations ( 5 . 4 . 7 ). 



q = (l/L^)(a+c) 

cl 

= [1/(6.578)^] [5.895 + 2.440] 



= 0.2049 ft 



q^ = (l/2L^)(d/c-b/a) 

= [1/2(6.578)^] [5.895/2. 440 - 2.440/5.895] 
= 0.0246 ft"^ 



& 

■ 



I 



■f 



1 
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= 2.1+1K)/(6.3T8)^ 

= 0.0600 ft'^ 

= (l/2L^)(d/c+t/a) 

= [1/2(6.578)^] [5.895/2. »+40 + 2.W+O/5.895] 
= 0.0548 ft'^ 






0^ = M^(0.2049 + 0.0246 - 0.0246 h^) 

- M2(0.0600 + 0.0548 - 0.0246 h^) 

Qj = - M^(0.2049 + 0.0246 - 0.0246 h^) 

+ M 2 ( 0.0600 + 0.0246 - 0.0548 h^) 



(8.U.9) 



From the results of the study presented in Appendix A one may 
conclude that the hydrostatic pressure ratios for soil pressures corre- 
sponding to the fully active and the fully passive states of stress in the 
soil may be calculated to a reeisonable degree of approximation using the 
Coxilomb formulas (cf. reference (l4)). 




(8.4.10) 



vhere K. = hydrostatic pressure ratio corresponding to the "Active 

State of Plastic Equilibrivun" in the soil (cf. Section 2.2) 

? = angle of internal friction of the soil 

w = angle of inclination with respect to the vertical of 
the lower bars of the model 

6 ss effective friction angle between the eirch and the soil 
i = slope of the ground surface with respect to the horizontal 



[I 



> 
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K. = 



cos^(50° - 22°29.6*) 



cos^(22°29.6* )cos(20°+22°29.6' ) 



1 + 



\/ 



/ 



s in ( 20‘^+50^ s in 50° 



n2 






= 0.510 




COS^(?+t*j) 




(8.4.11) 



vhere 



Kp = hydrostatic pressure ratio corresponding to the "Passive 

State of Plastic Equilihrium" in the soil (cf. Section 2.2). 






cos^(50°+22°29.6' ) 



cos^(22°29.6’ )cos(22°29 



.6' -20°) j^] 



sin(50 +20 )sin 50 



-- l2 



cos ( 22°29. 6 ' -20° )C08 ( 22°29 . 6 ' ) 



J 



= 5.455 

The soil pressure ratios corresponding to the directions of Pj^ and Pp for an 
"at rest" state of stress in the soil (cf. Section 2.2) may he found using 
Equations (4.1.11) and (4.2.11). 

K^. = \J 1+k^ cos I (w+6)-arctan(l/k )] 

wA O O 

= \/l+(0. 5 )^ cos [ (22°29.6'+20°)-arctan(l/o«5)] 

= 1.051* 

Kqp = \/l+k^ cos [arctan(l/k^) - (<*^-6)j 

= \/l+(0.5)^ cos [arctan(l/0.5 ) - (22° 29. 6' -20°)] 

= 0.545 



The deformation corresponding to full development of the "active" 
state of stress in the soil is taken as an inward rotation of the lower bar 



An error of 1,5'P in Kq. is introduced at this point. It was found to cause 
a change in the calculated response of less than l*^. 



I 



! 



i 
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of the model, of 0.001 radian (cf. Section 4.4). The deformation corre- 

sponding to full development of the "passive" state of stress in the soil is 

then taken as an outward rotation of the lower bar of the model, , as 

’ ty 

given by Equation (4.4.8). 



^ - 



K, 



OP, 



I' = i { 
hy hy\^ - K/ 

^ ^'‘1.051 - 0.510'' 

= 0.005342 radian 

While is not known to this degree of precision, it is convenient to 
retain a sufficient number of significant figures that calculation of the 
slope of the "elastic" line in the soil resistance -deformation relationship 
(cf. Figure 4.12) will give the same resxilt for the build-up of both the 
"active" and the "passive" states of stress in the soil at a given over- 
pressure level as assumed in Section 4.4. 

Expressions for the soil "stiffness" may now be obtained from 
Equations (4.4.10) and (4.4.12). 

7 d^ 

k' = " ■ (KQ^-K^)[cos(w+6)+(b/a)sin(w+6)] 

iy 

0.120(12.500)^(1.051-0.510) r .,^o^ 

= ^(5.S3)(0.00i) 29.6 +20 ) 



+ (|i^)sin(22°29.6’+20°)] 



= 5647 kips/radian 



k" = (5/^3d^)k' 

= [3(3647)/(0.120)(12. 500)] [0.144] 

= 1050 kips/radian-psl 
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Kgj » 5647 - 1050 (8.4.12) 

^sj = 

K' = 5647 - 1050 F’ (8.4.15) 

fc> J S J 

Alternatively, the soil stiffness may be expressed using Equations (4.4.14) 
and (4.4.15). 

K ./a = (k' -k" p .)/a 
sj' 

= (5647-1050 Pgj)/5.895 

ic ./a = 618.9 - 178.2 p . (8.4.14) 

SJ *^sj 



K' ./a = (k'-k" p' ,)/a 
sJ ^sj*" 

k; / a = 618.9 - 178.2 F'. (8.4.15) 

S J B J 

Hie soil resistance -deformation diagram for this problem is shovn in 
Figure 8.4. 

The mass concentrations corresponding to the horizontal and verticeil 
components of the inertial forces acting at the Joints of the model are given 
by Eqviations (5.1.15). Hie geometry of the arch-soil "strectxxre" is shovn 
in Figure 8.1. Hie geometry of the soil masses assumed to move with the 
arch is shown in Figures 5.1 sjid 5*2. The limiting vertical dimensions of 
the soil masses at the haunches are determined using Equations (5*1.5) and 

(5.1.4). 

‘^ml ^ t(a/2)K;-Hi^] < (b+d) 

[(L »| 2 1 ) + 2.440 + 4.167] < (2.440+5.895) 

9.554 > 8.555 
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use : 



use: 



= 8.333 ft 

d„2 = l(c/2)4d^J< 2d 

? 

[(2.1+40/2) + I+.I67] <2(5.893) 
5.387 < 11.786 

d„2 - 5.387 ft. 



' 2i '“‘<5 * l>] 

+ 7 h l/g 



a a 



[8.333(£fi2) . (8.533 - 

. (8.333 . ? . -8g3 . 1 ^ (8.333 )=oot 60°j 

0.150(8)(6.378) 

12(32.2) 



= 0.1782 k-sec^/ft 



(^U = it |H^) - ¥)(|) t d^ cot(5 - |)j 

= afegfr [8-535(24!!^) . (8.333- 

+ (8.333 . 5 - 8?3 t g-'*'^ )(lt|2i) + (8.333)^ oot 30°j 

0.150(8)(6.378) 

12(32.2) 

= 0.3276 k-sec^/ft 









m 




\ 
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= 



7 d 

i- - 1> 



7 h L 
c \ a a 
+ 



g 



_ 0.120(3.893) /c ,87 - + 0.1^0(8)(6.378) 

32.2 ^ 12(32.2) 

= 0.1247 k-sec^/ft 



(m^)in = (m^)^ 

= 0.1782 k-sec^/ft 



^“3 ^out 



^“l^out 

0.3276 k-sec^/ft 



mi = 



!s 

2g 






1 7 h L 



ml 2'' 2 



ml 2 '^ 2 * 



g 



. 2^ [8.555(4^!^) . (8.353- 
. (8.J55 . ^ 9 . . ,. 13P . || K 6. 1 78) 

= 0.1035 k-sec^/ft 



“5 " “2 



= 0.1247 k-sec7ft 



"“6 = “4 



= 0.1035 k-secVft 



•nie equivalent masses for the tvo -degree -of -freedom system can 
now be found using the mass concentrations calculated above and Equations 
(5.2.6). The calculations for these parameters can be shortened consider- 
ably by noting that (103^ = (“5)^0^ ^“l^out ^ ^“3^out^ ™4 ° “6 

V = V . 
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.. ( \ ^ /bv2 ^ /d bx 

in _L ^ ^ 



= 0.1782 * 1^1 . , 0.1055(ft^: 

= 0.5422 k-secVft 



2 



l,m| . , I, ' * T + if <1 • 

. 0.J276 . ^ [i - g^J . o.io„(|^: 



(5.895)‘ 

^ "4^ (|tlS - 

= 0,4916 k-sec^/ft 






%m| 



m 



= 1,2,5A 



“*o X. 2 , , 2 



_ 0. 1247 

5 



1 - ^ ^ (|tli - 



(5.893)“ 



= 0.1054 k-sec^ft 

^5>ml , ^ " ^l,m| ^ . 

'm = 1,5 ' 'm = 5,4 

= 0.4916 k-sec^/ft 

= 2,4 = 1,2 

= 0.5422 k-sec^/l't 



= 1,2,5,4 ^^“L = 1,2,5,4 

= 0.1054 k-sec^/ft 



2 



2.440x^ 



I 
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Bie equivalent masses for the four motions of the system (cf. Figure 5*5) 
may nov be written. 

Asymmetric -Ripht : m^ . = 0.5422-0.1034 h^/h\, 

= 0.4916-0.1054 



(8.4.16) 



Symmetric -In ; 



m. 



1.2 

S,2 



Symmetric -Out : 



“1,5 

S,5 



Asymmetric -Left : m. 



I.'* 



m 



‘5 A 



0.5422-0.1054 h^h^ 

0.5422-0.1054 h./ii, 
1 5 

0.4916-0.1054 h^h^ 
0.4916-0.1054 h^/h^ 

0.4916-0.1054 kJ\ 
0.5422-0.1054 ii^/h^ 



(8.4.17) 



(8.4.18) 



(8.4.19) 



The equivalent loads on the two -degree -of -freedom system are given 
by Equations (6.3.2). 

= I Kp^(a+c) + |(b/a)(b+d) 

= 0.144 |^|(J)(5.895 + 2.440) + |<f7^)(2.440 + 5-893) 

= 0.3985 k/psi 



^ cK (1 - ^) 
2 4 pr ad' 



= J(2.440)(J)(0.144)[1 
«= 0.0182 k/psi 



Jj = I d(d/c - b/a) 



(2.440)7(5.893) ] 



= (5.893/2.440 - 2.44 o/ 5.893)(0.144) 



= 0.8491 k/psi 



(8.4.20) 



= - 0.5985 - 0.0182(p 2^ ■ ^2b^ 0.8491 Pg 

= - 0.0182 (^2^ - Pgb^ ■ °*®^91 ?2 + 0.3985 P3 



Since Equations (8.4.12) throu^ (8.4.15) for equivalent soil stiff- 
nesses and Equations (8.4.20) for equivalent loads are written in terms of 
the average overpressures acting on the loading regions defined in Section 6.1, 
it is necessary to find these overpressxires using Equations (6.2.1) and the 
proper set of time parameters from Section 6.1. Table 8.1 gives the velocity 
of the shock front at the 100-psi overpressure level as 2,950 feet per 
second. This is approximately equal to the ass\imed value of 5»000 feet per 
second for the seismic velocity of the soil. From the results of the study 
presented in Appendix B, one may assume, for the depth considered in this 
example, that the shock front in soil is essentially vertical and a continua- 
tion of the shock front in air. This corresponds to Case I, Section 6.1.1. 

The time parameters for this case are given by Equations (6.I.I) and (6.1.2). 



t 

oa 



t 



= - (b+2d^cot 



jt 

IT 




2U 



= - (2. 440+2(12. 500) cot 30°j/2(2950) 



= - 0.0078 seconds 



t • = ^ (t , + t , ) 

o' 2 oa' Ob' 



= I (-0.0078+0.0015) 



= - 0.0032 seconds 
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- (b-2bdya)/2U 

- [ 2 . 440 - 2 ( 2 . 440 )( 12 . 500 )/ 5 . 895 ]/ 2 ( 2950 ) 
+ 0.0015 seconds 

- (b+ 2 d^cot[^ + |])/ 2 U 

- [ 2 . 440 + 2 ( 12 . 500 )cot 60 °]/ 2 ( 2950 ) 

- 0.0029 seconds 



I (‘oa " U) 

I (-0.0029+0.0015) 
- 0.0008 seconds 



0 



I (‘la ‘ ‘lb> 

I ( 0 + 0 . 0014 ) 

+ 0.0007 seconds 



(b+d)/ 2 U 

( 2 . 440 + 5 . 895 )/ 2 ( 2950 ) 

+ 0.0014 seconds 



lb 

+ 0.0014 seconds 



2 ^^2a '*^2b^ 

I (+0.0014+0.0054) 

+ 0.0024 seconds 
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(b+5d)/2U 

[2.4itO+3(5.895)J/2(2950) 
+ 0.0034 seconds 

2 ^"^2a ^2^ 

I (+0.0014+0.0024) 

+ 0.0019 seconds 

(b+2d)/2U 

[ 2 . 440+ 2 ( 5 . 895 ) ] /2 ( 2950 ) 

+ 0.0024 seconds 

I (+0.0024+0.0034) 

+ 0.0029 seconds 



2b 

+0.0054 seconds 



I " ‘3b> 

I (+0.0034+0.0048) 
+ 0.0041 seconds 



(b+2d)/U 

12.440+2(5.895)J/2950 
+ 0.0048 seconds 



(-b+2B-2bdya)/2U 

[-2. 440+2(16. 66t)-2(2. 440)(i2. 500 )/5.893]/2(2950) 



+ 0.0035 seconds 
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2 

I (+0.0055+0.0077) 

+ 0.0056 seconds 



( -b+2B+2d^cot (J - |])/2U 

[-2. 440+2(16.667 )+2(l2.500)cot 60°]/2(2950) 

+ 0.0077 seconds 



I ) 

I (+0.0035+0.0126) 

+ 0.0080 seconds 



( -b+2B+2d^ cot [ ^ )/2U 

[-2. 440+2(16.667 )+2(l2. 500 )cot 30°j/2(2950) 

+ 0.0126 seconds 

- b/2U 

- 2.440/2(2950) 

- 0.0004 seconds 



I (-0.0004+0.0052) 

+ 0.0024 seconds 



(-b+2B)/2U 

( -2. 440+2( 16 . 667 ) j/2( 2950 ) 



+ 0.0052 seconds 
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The solution of Equations (6.2.1) may be presented conveniently 
in tabular form as shovn in Table 8.6. Since the integration of the equa- 
tions of motion is to be performed on the Illiac, each of the overpressure- 
time functions is approximated by a series of straight-line segments, 
therefore, the overpressure -time coordinates given in Table 8.5 are used for 
the general overpressure -time function, p(t'), with the exception that the 



- ^1b'^a 

maximum value of each overpressure, i.e. * p( “ 2 taken 

directly from the p(t') vs. t' curve itself. This is easily done and improves 
the approximation sli^tly. Thus, Table 8.6 contains the coordinates of the 
end points of the series of strai^t-line segments used to approximate the 
average overpressure acting on each of the loading regions. 

The equations of motion may be written using Equations (7.1.^). 



l,m h- 4, -V V 
’ l,m l,ra ji,m 

, ^ Q.k9l6 

(0.5422)(0.k9l6)-(0.l034)^ 

= 5.1215 ft/k-sec^ 



Itm 



l,m u, V 

' l,m J),m l,m 

, 0.1054 



(0.5422)(0.4916)-(0.1054)^ 
= 0.6565 ft/k-sec^ 



F = It™ 

P 0.1054 



(0.5422)(0.49i6)-(0.1034)^ 
= 0.6565 ft/k-sec^ 
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n 4 -V V 

. ^ 0.5^22 

(0.5^22)(0.4916)-(0.1054)^ 

= 2.1727 ft/k-sec^ 



The remainder of the coefficients are calculated in the same manner. 
Inserting the coefficients into the equations of motion for the four motions 
of the system (cf. Figure 5*5) gives: 



Asymmetric -Rirfit : 

h^ = 5.1215(H^-\-Qg^) + 0.6565(H^-^-Qg^) 

= 0.6565(H3^-Q3^-Qg^) + 2.1727(ii^-Q^-Qgj) 



Symmetric -In : 

h^ = 5.2159(H^-Q^-Qg^) + 0.9717(aj-<^-Qg^) 
h. = 0.9717(H^-Q^-Qgl) + 5.2159(H^-Q5-Qg^) 



Symmetric -Out : 

= 2.1283(H^-Q^-Qg^) + 0.4477(H^-Q5-Qg3) 
= 0.4477(H^-Ql-Qgl) + 2.1283(H^-Q5-Qg^) 



Asymmetric -Left : 

h^ = 2.1727(H^-Q^-Qg^) + 0.6565(H^-Qj-Qg^) 

h^ = 0.6565(H^-Q3^-Qg3^) + 3.12l5(^-^-Qg5) 



(8.4.21) 



(8.4.22) 



(8.4.25) 



(8.4.24) 



8.4.2 Investigation of General Yielding and Buckling of Arch Rib 

Before proceeding with the integration of the equations of motion, 
it is advisable to Investigate the effects of the thrust in the arch rib. 

The average overpressure acting on the arch is given in Table 8.6 as p . 

£l 



' 



r 
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The thrust in the arch, T, caxised by this average oveiT>ressure is given 
approximately by Equation (3.5»l)> ^ich is repeated here for convenience as 
Equation (8.4.25). 



T = I 5; B 



(8.4.25) 



where p^ = effective uniform radial overpressxtre acting on the arch 
B = span of the arch 



Since the value of being used in this problem is one -fourth, the radial 
overpressure acting on the arch will vary from 255» of p at the abutments to 

3 , 

lOOf) of p^ at the crown. Thus, because the arch is semi-circular, this gives 

an average radial overpressure on the arch of 62. 5^^ of the vertical overpressure. 

However, this neglects the thrust induced in the arch as the sides are forced 

outward against the soil during latereJ. response. Therefore, some higher 

proportion of the vertical overpress\ire is appropriate. For the approximate 

CEilculatlon here, p' is taken as three-fourths of p . 

' a ^a 

A measure of the capacity of the arch rib cross section to carry 
thrust loads may be obtained from Equation (8.4.26). 



T = 0.85f A + f , A 
y c c dy s 



(8.4.26) 



where T = thrust corresponding to general yielding of the arch rib 

f^ = 28-day strength of a standard 6"xl2" concrete test cylinder 

A = area of concrete in the arch rib cross section 
c 

A = area of steel in the arch rib cross section 
s 

f^= dynamic yield strength of the reinforcing steel 



Equation (8.4.26) neglects bending moment in the arch rib and is therefore 
approximate; however, because the mode of response being considered here is 
essentially a breathing mode, and further because of the high degree of soil 
restraint, bending in the arch rib is not of prime importance in this 
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investigation. Substituting values into Equation (8.4.26) gives: 

T = 0.85(5. 5)(96. 00-0.24) + 50(0.24) 

= 297 kips 

The average overpressure, p , corresponding to general yielding 

& 

of the arch rib may now be fovind using Equation (8.4.25). 

(p' ) 

\ _ -^a max 

^^aW ■ 0.75 

2T 

= 

0.75B 

2(297)(1000) 

~ 0. 75(16. 667)(144) 

= 550 psi 

Since 550 psi is greater than the maximum value of p in Table 8.6, i.e. 

& 

p = - 97*7 psi, the arch should not fail by genereil yielding. 

A qualification should be made here regarding the use of an 
"average" radial overpressure in this case. Since the crown is subjected to 
an overpressure level greater than the "average" value, local yielding may 
occur \diich could lead to general fedlvire of the arch. This is particularly 
true when the depth of cover over the crown is zero or very small. 

In considering the elastic stability of the arch rib, one must 
keep in mind that buckling is not instantaneous. For the higher yield 
nuclear weapons, however, the duration of the positive phase of the shock 
wave is quite long when compared with the characteristic response time of 
the structvire so that buckling is a possibility v^ich should be investigated. 
If one assumes that the arch rib may be approximated by a bar pinned at the 
crown and at the abutment and restrained by the soil, the eleistic stability 
of the arch rib may be investigated using Equation (8.4.27).^ 



1^8 



vhere 






cr 



cr 



(m^ 






-) 



t in « El 

UEI , 2 s 

l6m^I 



(8.4.27) 



= thrust corresponding to buckling of the arch rib into a 
simple sine curve 

E = modulus of elasticity of the arch material 

I = moment of inertia of the arch rib cross section 

t = developed length of one -half of the arch 

m = number of half sine loops in the buckled configuration 

= soil stiffness, here assumed to act in both directions 
R = radius of the arch. 



Equation (8.4.28) may be obtained from Equation (8.4.27) ty equating the 
critical thrust values for vhich the number of half-waves changes from m 
to m+1. 






ISei 



(8.4.28) 



The minimum soil stiffness corresponding to the maximum thrust 
in the arch is obtained by noting that the maximum value of p in Table 8.6 
occurs at t = 0.0052 seconds. While the configuration of the arch at 
t = 0.0052 seconds is not yet known, a "lower bound" for the soil stiffness 
may be obtained by using the minimum soil restraining pressure acting at 
this time. From Table 8.6 it appeeurs that p^^ has an absolute magnitude, 
i.e. 18 psi, less than all of the other soil restraining pressures at 
t = 0.0052 seconds. 



- P 



S3 



_ ( 0.0052-0.0033 )/^ u) 

~ ^O.Oli^-O.OoSV''^*^'^ 



= - 18 psi 
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The soil stiffness is now determined using Equation (8.4.15). 

ic; ya = 618.9 - 178.2 

iCgj/a = 618.9 - 178.2(-18) 

= 5827 k/ft^ 

= 5827/0.144 
= 26,600 psi 



The number of half sine loops corresponding to these conditions is obtained 
from Equation (8.4.28). 

= llii 

= 26.600(100)^ 

16 (5. 82x10^) (268. 6) 

= 162.05 



m(m+l) = 12.729 
m = 3.103 



The criticed thrust may now be found from Equation (8.4.27) by substituting 
the integral value of ra which gives the lowest value of T^^. The value m = 4 
gives the minimm value in this case. (Note: Differentiating Equation 

(8.4.27) with respect to m, setting dT^^dm = 0, and solving for m gives 
m = 3*568). 



cr 



4EI , 2 
— 5- (m + 

R 



l6m^I 



) 



_ 4(3. 82x10°) (268. 6) 162.03 \ 

(100)^(1000) ~i6 



= 10,725 kips 
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The average overpressure, p^, corresponding to buckling of the arch may now 
be found using Equation (8.4.25). 



2T 

^ - cr 

^^a^cr “ 0.75B 




= 11,900 psi 



This approach, although approximate, indicates that the arch should not 
buckle elastically prior to general yielding of the arch rib itself. Note 
that this is an "elastic" analysis so that the actual buckling of the arch 
in this case is associated with general yielding of the arch rib as described 
in the first part of Section 8.4.2. 



8 . 4.3 Integration 

Since it has now been determined that neither buckling nor general 
yeilding of the arch rib should occur under the given loading, it is possible 
to proceed with the integration of the equations of motion. The results 
shown in Appendices C and D were obtained using a digiteil computer (cf. Sec- 
tion 8 . 5 ). Two of the problems were eilso integrated by hand in order to 
check the computer solutions. 

A cycle of integration is given below to illustrate the operations 
involved. The nvnnerical Integration formulas used are those derived by 
Professor N. M. Newmark and presented in references (4), ( 5 )^ (6) and (7). 
Equations (8.4.29) are expressions for the velocity and displacement of mass j 
at the end of time interval At in terms of the acceleration at the end of the 
interval and the acceleration, velocity and displacement at the beginning of 



the interval. 



■vriiere 



Ihl 




( 8 . 4 . 29 ) 



h . ^ h . ^ h . = acceleration, velocity and displacement, 

j,n j,n j,n respectively, of mass j at time "t = 



^J,n+1^ ^j,n+l 



h. = acceleration, velocity and displacement, 
respectively, of mass J at time t = 



At = time interval. = t , -t 

n+1 n 

f)' = parameter introduced to indicate how much of the acceleration 
at the end of the time interval enters into the relations for 
the displacement at the end of the time interval. 



Since it is desired to use a linear approximation to the acceleration curve 
within each time interval, the value P' = l/6 is used in Equations ( 8 . 4 . 29 ). 



O.OOOl seconds is used. This time interveil was actually used for the computer 
solutions in order to detect the rapid changes in the loading as the shock 
wave traverses the structvure. Solutions were also obtained using time 
intervals of 0.001 seconds for several of the problems ^ich had more slowly 
varying load-time functions in order to check the significance of round-off 
error. These solutions were nearly identical with the solutions obtained 
using At = 0.0001 seconds. 



interval other than the first. Therefore, for this example three time 
intervals have already been integrated, so that at time t^ = O.OOO 3 seconds. 




(8.4.50) 



In the following Illustrative example a time intei^ral. At, of 



For illiistrative purposes it is preferable to consider a time 




i 







iw 



1U2 



the displacements, velocities and accelerations of the equivalent masses 
are: 



h = +0. 597x10'^ ft 
J.^ n 

h- = +3. 927x10 ft/sec 
l,n 

h^ = +26.1*74 ft/sec^ 



h, = +0.084x10"^ ft 

h, = +0.855x10"^ ft/sec 

ii, = +5.568 ft/sec^ 

3,n ^ 



Assume values of the accelerations at the end of the interval, i.e. at 



time “ 0.0004 seconds, as follows: 



\n+l " +55.279 ft/sec ^ S,n+1 " +'^•^20 ft/sec ^ 

The velocities at = 0.0004 seconds may be foxand using the first of 

Equations (8.4.30). 

hi,n+i = (+3. 972x10 ’5) + I (+26. 474+35. 279)(10'^) 

= +7.060x 10 5 ft/sec 

h, , = +1.485x10 "5 ft/sec 
3,n+l 

The displacements at " 0*0004 seconds may be calculated using the 

second of Equations (8.4.30). 

^l,n+l " (+0-597x10 ■^)+(+3. 972x10 ■5)(io‘‘^) 

+(i/3)(26.474)(io'®)+(i/6)(35.279)(10"®) 

= +0.941x10"° ft 

h^ , = +0.198x10"^ ft 

3,n+l 



The rotations of each joint are found using Equations (8.4.5), (8.4.6) and 

(8.4.7). 

= 0.4098 h^ - 0.1698 h^ 

= [0.1*098(0.941) - 0.1698(0.198)] 110"^] 

= +0.352x10“^ radians 
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-6 

1^2 ® -0.305x10 radians 
—A 

'i'j = +0.079x10 radians 



The resisting moments at each hinge are found from Eq\iations (8.4.8). 

«i = A 

= (+-0.352 x 10"^)(1117.2) 

= +0.393x10"^ ft -kips 

= -0.340x10"^ ft kips 

M = +0.088x10 ^ ft kips 



!Hie equivalent lateral resistance of the arch is nov found from Equations 

(8.4.9). 



0^ = M^( 0.2049+0.0246 h^-0.0246 h^)- M2(0.0600+0.03l+8 h^ 
-0.0246 h^) 

= [ +0 . 393x10 '^ ] 1 0 . 2049+0 . 0246 ( +0 . 941x10 ) 

-0. 0246 ( +0.198x10"^)] - [-0.340x10"^] [0.0600 
+0. 0348(+0. 941 x 10‘^ ) -0. 0246(+0. 198x10"^ ) ] 

= +0.101x10 ^ kips 



0^ = -0.038x10 ^ kips 

1 —6 

The displacement h^ = +0.941x10 ft indicates that the windward haunch, 
joint 1, is displaced inward into the arch. From Table 8.6 the average over- 
pressure acting on the "active" soil failure wedge adjacent to joint 1, p^^, 
at t^^^ = 0.0004 seconds is 77.2 psi. 



5sl = * Ciij - 67.9) 



= 77.2 psi 
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The soil stiffness corresponding to this condition is given by Equation 
(8.4.12). 

^sj = 5647 - 1050 

iCgi = 5647 - 1050(-77.2) 

= 84,710 kips/radian 

The displacement h^ = +0.198x10 ft indicates that the leeward haunch, 
joint 5^ is displaced outward against the soil. From Table 8.6 the average 
overpressure acting on the "passive" soil failure wedge adjacent to joint 5, 
p^^, at t^^^ = 0.0004 seconds is zero because the shock front has not yet 
arrived at point 4a (cf. Figure 6.l). The soil stiffness corresponding to 
this condition is given by Equation (8.4.15). 

ic' . = 3647 - 1050 p' , 

sj ^sj 

= 5647 - 1050(0) 

= 5647 kips/radian 

The rotations of the windward and leeward pressure sxirfaces are obtained 
using Equations (4.4.15). 

= h^/a 

+0.941x10“^ 

" ^^95 

= +0.l60xl0~^ radians 
= h^/a 

_ +0.198x10“^ 

5.895 
. -6 

= +0.054x10 radians 



Since these rotations are less than the yield veilues calculated previously. 



1^5 



i.e., t - 0.001 radians and = 0.0053^2 radians, and no "permanent set" 
has yet occurred in the soil, the equivalent soil resistances may be calcu- 
lated as follows: 



= 84,710(+0. 160x10"^) 
= +15. 55x10"^ kips 



Q, 



■s5 



^s5 S 



= 3, 647 (+0.054x10"^) 
= +0.124x10 ^ kips 



The equivalent loads on the system may be determined from 
Equations (8.4.20). Ihe average overpressures loading the structure at 
t^^j^ = 0.0004 seconds are given in Table 8.6. 



■Pi = = 28.4 psi 



= -0.5985 - 0.0182(p2^-P2^) + 0.8491 P2 

= -0.5985(-28.4) 

= + 11.52 kips 

% = o 

The fact that both of the velocities, h^ and at time 

t^^^ = 0.0004 seconds are positive indicates that both havinches are moving 
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toward the right. Thus, the motion of the structure is "Asymmetric -Ri^t" 

(cf. Figure 5*3) and the equations of motion for the system at this time 

are given by Equations (8.4.21). By substituting the effective load and 

resistance parameters calculated in this section into Equations (8.4.21), 

it is now possible to check the asstnned values of the accelerations at time 

t , = 0.0004 seconds. 
n+1 

= 5.1213 [11. 32-(+0. 101x10 ‘^)-(+13. 55x10"^)] 

+0.6565[0-(-0. 038x10 "^ ) - ( + 0 . 124x10 ) ] 

= +35.29 ft/sec ^ 

h, = +7.42 ft/sec^ 

3, n+1 

Since these values are essentially the same as the assumed values, l.e. 
h^ = +35*279 ft/sec^ and h^ = +7.420 ft/sec^, additional trials are 
not required and the velocities and displacements at t = 0.0004 seconds 
which were calculated above are correct, llius, the integration may now 
proceed to the next time interval. The results of the solution of this 
problem up to t = O.O9 seconds are given in Figure D.3. 

8.4.4 Rib Shortening 

The effect of rib shortening is illustrated in Figures D. 24, D.25 
and D.26. Expressions for the "rib -shortening" corrections are obtained using 
Equations (3 .5. 8 )* In "the following derivations an average imiform radial 
pressure on the arch equal to 75^ of the vertical overpressure in the soil 
is assumed (cf. Section 8.4.2). 









m 
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_ 1 BL^^c 

^ 1 ~ 2 [_AE(axi-bc ) 

= . i p' [ 16.667(6.^78)^(2. 440)(12) 

^ ^ U (5. 82x10^) (5. 895^-2. 440^ )J 

= - 11 . 29 x 10 "^ 

6h^ = -8.47x10"^ p 
1 

bh^ = 0 

5 h^ = - 6 h^ 

6h^ = +8.47x10"^ p 
5 a 



5h,. = 






2 Pa 



BL^d 



AE(ad-bc) 



^F' 

2 ^a 



16.667(6.378)^(3. 893)(12) 

(- 8 ( 3 . 82x10° )( 5 . 893 ^- 2 . 440^ )J 



6h,. = 



27.26x10"^ p' 
a 

20.45x10'^ p 



6 h 



5 



= |p;( 



BL^ x 

cAE'' 



= 1 p. [ 16.667(6.378)^(12) ' 

^ ^ U. 440 (8) ( 5 . 82x10^)- 

= 55 x 10 '^ F; 

6h^ = 40.92x10"^ F 
? a 



6 hg = 6 hj^ 

6hg = 20.45x10"^ F^ 



(8.4.31) 

(8.4.52) 



(8.4.53) 



( 8 . 4 . 54 ) 



(8.4.55) 



(8.4.56) 
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The parameters required to plot Figure D. 25 are calc\ilated below 
in order to illustrate the operations involved. The time selected for this 
computation corresponds to the first maximum dovmward displacement of the 
crown. From Figure D.5; this occurs at t = O.OI 98 seconds. The displace- 
ment components at this time may be calculated using the data plotted in 
Figvire D.5 and Equations (8.4.1) through (8.4.4). 

h^ = -5.8x10'^ ft 

hg = +0.4l45(h^+h^) 

= +0.4i45(-5.8+5.0)(10"^) 

= -0.4x10*^ ft 

h^ = +5.0x10'^ ft 

hj^ = -0.4141 h^ 

= -0.4141 (-5.8x10"^) 

= +2.4x10'^ ft 

h^ = -10.8x10"^ ft 

h^ = -f0.4l4l h, 
o 5 

= 40. 41 41 (+5.0x10'^) 

= +2. 1x10 ft 

The average overpressure acting on the arch at t = O.OI 98 seconds 

may be obtained either from Trble 8.6 or from a plot of p (t) vs. t. For 

& 

this example a value of p = -88.1 psi was obtained from a plot of p (t) vs. t. 

Q, & 

This value is now substituted into Equations (8.4.31) through ( 8 . 4 . 56 ) in 
order to obtain the "rib -shortening" corrections at t = O.OI 98 seconds. 



D 
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-f, — 

= -8.47x10 
= -8. 47x10 “^(-88.1) 
= +0.7x10"^ ft 



Shg = 0 




- oh. 



-0.7x10 



ft 



6h|^ = +20.45x10"^ 

= +20. 45x10 ■^(-88.1) 
= -1.8x10"^ ft 



5h^ = +40.92x10"^ p 
j a 

= +40. 92x10 "^(-88.1) 
= -3.6x10"^ ft 



5h^ = 5h^ 

= -1.8x10“^ ft 



The net displacement components, h,, of the joints of the model 

J 

may now he obtained by adding algebraically the "rib -shortening” corrections 
to the displacement components calculated above. 



h . = h . + 5h . 
J J J 



(8.4.57) 



h^ = (-5.8+0.7)(10“^) 

= -5.1x10"^ ft 



hg = (-o.4+o)(io"5) 



= -0.4x10"^ ft 






h 

li 

,'J 







|: 

i 



i 
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= (+ 5 . 0 - 0 . 7 )( 10 “^) 

= +4.5x10"^ ft 



hj^ = (+2.4-1.8)(10"^) 



= +0.6x10"^ ft 



= (-10.8-5.6)(io‘^) 
= -l4.4xl0'^ ft 



= (+ 2 . 1 - 1 . 8 )( 10 "^) 

= +0.5x10"^ ft 




the square root of the sxim of the squares of the net horizontal and vertical 
displacement components at that joint. 



8.4.5 Calcvilation of Natural Periods of Vibration 

As an illustration of the calculation of the natureQ. periods of 
vibration of the system, the natural period of vibration is calculated below 
for the structure as it approaches its first maximxan response at the leeward 
haimch. Because the "stiffness" of the soil varies with the overpressure 
restraining the soil wedges (cf. Section 4.4), the time selected for this 
example is the time corresponding to maximum soil "stiffness" up to this 




(8.4.38) 



= 5.1x10"^ ft 
Ag = 14.4x10"^ ft 
A^ = 4.3x10"^ ft 



instant. 
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From Figure D.J one determines that the motion of the strxictxire 
is "Symmetric -Out" (cf. Figure 5*3) as the system approaches its first maxi- 
mum displacement at the leeward haunch and that this maximum occurs at 
t = 0.0173 seconds. Thus, the haunches are displaced outward and a "passive" 
state of stress is developing in the soil on both sides of the arch. From 
Table 8.6, or from plots of the overpressia*e-time functions, the maximum 
soil restraining pressures acting on the "passive" soil failure wedges during 
the interval between t = 0 and t = 0.0173 seconds are found to occur at 
t = 0.0126 seconds when p^^^ = -89. 0 psi (obtained from plot) and p^^ = -96.4 
psi. 

The soil stiffnesses corresponding to these restraining pressures 
are obtained from Equations (8.4.15). 

K^^a = 618.9 - 178.2 p’j 

ic^^/a = 618.9 - 178. 2( -89.0) 

= 16, 480 k/ft 

K^^/a = 618.9 - 178. 2( -96. 4) 

= 17,800 k/ft 

The arch stiffness parameters required for substitution into the 
frequency equation may be found by referring to Equation (7.2.3). 

Kp 

T = — — (5+a/c+b/d+2c/a+bc/ad) 

- HIT. 2 fc . 5.893 . 2,440 2(2.4^ {2.kkof '\ 

■ 2(6.578)2 r ^ ^ 5:^ (5.393)2j 



=121.2 k/ft 






S 




152 



Kg 

t = — - (5 + a/c - b/d - bc/ad) 

= 1117.2 r 5.895 _ 2.440 _ (2.440)^ 1 

2(6.378)^ L 2TTO 5.895 (5. 895)^-! 

= 66.5 k/ f t 

The equivalent mass parameters for the "Symmetric -Out" motion are 
those appearing in Equation (8.4.18). 

= d , = 0.4916 k-sec^/ft 

V = V = 0.1034 k-sec^/ft 

1^3 

The parameters obtained above may now be substituted into Equation 
(7.2.9) and the resulting equation solved for . 

1 d, d, -V V j - u (t +K' /a) 

L l^m 5^ra l>ni 3>niJ m |_ l,m 1 s3 

+ d, (T,+K',/a) - T,(v'- +v ) 

1 sr ' 3 l>Di 5 >di'J m 

= 0 

[(0.4916)^-(o.1034)^] u)^ - [ 0.4916(121+17,800) 

+ 0.4916(121+16,48o)-66. 5(0. 1034+0. 1054)] 

+ [(121+16,48o)(121+17,80o)-(66.3)^] 

= 0 

0.2510 -16,960 +279-4x10^ = 0 

m m 

^ = 29,000 rad^/sec^ 

3, 1 

= 44,400 rad^/sec^ 
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2 It II 

The value of w corresponding to "Symmetric -Out" motion may be found by 



substituting w, . and w , in t^lm, into either of Equations (7.2.8). 

-L 5^2 

(u, - 'tT-K'-/a)hT + (-V u)^ + T )h, = 0 

l,m m 1 si ' 1 ' l,m m 5 5 

For (1)^ - = 29,000 rad^/sec^: 

[0. 4916(29, 000)-121-l6,l|80]h^+[-0.1054(29,000)+66.3Jh^ = 0 

h^ = -0.80 h- 
5 1 

Since h^ and h^ are of opposite signs, the haunches are displaced in opposite 

II II 2 

directions and the configuration is "symmetric". (Substitution of ^ = 

5 , 2 

p . P 

44,400 Taa/e&c into Equations (7.2.8) gives h^ = +I.I6 h^, indicating that 
both haunches are displaced in the same direction and that the configuration 
is therefore "asyTnmetric"). 

2 2/2 

The correct value of w is thus 29,000 rad /sec from vhich 

m 

w = 170.5 rad/sec. This value may now be substituted into Equation (7. 2.10) 
m 

in order to find the natural period of vibration of the system at this 
particular time in this particular motion. 

T = 2 k/w 
m m 



T, = 2«/l70.3 
5 



= 0.057 seconds 



i 

I 
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TABLE 8.1 

SHOCK WAVE PROPERTIES 



Peak Side -on Overpressvire 


psi 


1000 


100 


50 


Distance from 1 MT Contact 
Surface Burst 


yds 


480 


1100 


1500 


Duration of Positive Phase 


sec 


1.48 


1.48 


1.60 


Velocity of Shock Front in Air 


fps 


8700 


2950 


2240 





OVERPRESSURE -TIME 


TABLE 8.2 

COORDINATES FOR 50-psi 


SHOCK WAVE 




t' 


sec 


0 


0.065 


0.20 


0.40 


0.80 


p(f) 


psi 


50 


40 


26.5 


15.5 


5.0 



TABLE 8.5 

OVERPRESSURE -TIME COORDINATES FOR 100 -psi SHOCK WAVE 



t' sec 0 0.070 0.205 0.595 0.702 

p(t' ) psi 100 60 51»5 16 6 



TABLE 8.4 

OVERPRESSURE -TIME COORDINATES FOR 1000 -psi SHOCK WAVE 



t' sec 0 0.01 0.04 0.12 0.40 

p(t') psi 1000 600 250 100 25 
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TABLE 8.5 

SUMMARY OF NUMERICAL EXAMPLES 



Depth of 
Cover Over 
Crovn 


Reinforced Concrete Arch 


Corrugated Metal Arch 


100 psi 


1000 psi 


50 psi 


100 psi 


1000 psi 


0 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


g- (Span) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


- 


^ (Span) 


( 1 ).( 5 ) 


( 1 ) 


( 1 ) 


( 1 ),( 2 ) 


( 1 ) 


1 (Span) 


( 1 ) 


( 1 ) 


( 1 ) 


( 1 ) 


- 



Notes: (l) = l/k, = 5OOO fps 

(2) = 1/4, 1/2, 5/4, 1 

(5) = 1000, 2400 , 5000 fps 

o 
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TABLE 8.6 

OVERPRESSURE -TIME RELATIONSHIPS 
oading Region Parameter Coordinates 



la-lb 


t=t' +0.0007 


sec 


0 




t' 


sec 


- 




-Pi 


psi 


0 


2a-2b 


t=t' +0.0024 


sec 


0 




t* 


sec 


- 




-P2 


psi 


0 


5a-5b 


t=t '+0.0041 


sec 


0 




t' 


sec 


- 






psi 


0 


2a-2 


t=t' +0.0019 


sec 


0 




t’ 


sec 


- 




-P2a 


psi 


0 


2-2b 


t=t’ +0.0029 


sec 


0 




f 


sec 


- 




-P2b 


psi 


0 


Oa-Ob 


t=t' -0.0008 


sec 


0 




t' 


sec 


0.0008 




-Psl 


psi 


67.9 


OA-Ob 


t=t' -0.0032 


sec 


0 




t' 


sec 


- 




-Psl 


psi 


82.6 


4a-l<-b 


t=t' +0.0056 


sec 


0 




t' 


sec 


- 




-Ps3 


psi 


0 


4a-4b' 


t=t' +0.0080 


sec 


0 




t' 


sec 


- 




■^s3 


psi 


0 


aa-ab 


t=t' +0.0024 


sec 


0 




t' 


sec 


- 




-Pa 


psi 


7 



0.0014 


0.0707 


0.2057 


0.3957 




0.0c 07 


0.0700 


0.2050 


0.5950 




99.4 


60 


51.5 


16 




0.0014 


0.0054 


0.0724 


O.2O74 


0.5974 




0.0010 


0.0700 


0.2050 


0.5950 


0 


99.2 


60 


51.5 


16 


0.0054 


0.0048 


0.0741 


0.2091 


0.5991 


-- 


0.0007 


0.0700 


0. ^050 


0.3950 


0 


99.4 


60 


51.5 


16 


0.0014 


0.0024 


0.0719 


0.2069 


0.3969 


-- 


0.0005 


0.0700 


0.2050 


0.5950 


0 


99.6 


60 


31.5 


16 


0.0024 


0.0034 


0.0729 


0.2079 


0.5979 




0.0005 


0.0700 


0.2050 


0.5950 


0 


99.6 


60 


31.5 


16 


0.0015 


0.0692 


0.2042 


0.5942 




0.0021 


0.0700 


0.2050 


0.3950 




98.5 


60 


51.5 


16 




0.0013 


0.0668 


0.2018 


0.5918 




0.0045 


0.0700 


0.2050 


0.3950 




96.4 


60 


51.5 


16 




0.0055 


0.0077 


0.0756 


0.2106 


0.4006 


-- 


0.0021 


0.0700 


0.2050 


0.5950 


0 


98.5 


60 


31.5 


16 


0.0035 


0.0126 


0.0780 


0.2150 


0.4050 




0.0046 


0.0700 


O.2O5O 


0.5950 


0 


96.4 


60 


31.5 


16 


0.0052 


0.0724 


0.2074 


0.3974 




0.0028 


0.0700 


O.2O5O 


0.5950 




97.7 


60 


31.5 


16 





I 



4 

\ 

3 





12.500 
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Figure 8.1. Reinforced Concrete Arch, 



b/u 




Figure 8.2. Moment-Rotation Diagram for Joints of Nfodel 






1 



I 
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Figure 6.3 • Arch Rib Cross Section 




Figure Q.h. Equiveilent Soil Resistance-Displacement Diagram 
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CHAPTER IX 
SUMMARY 



The resxilts of the nmerical investigation performed in connection 
with this thesis are summeirized in this chapter. Included also are a listing 
of the conclusions which have been drawn on the basis of this investigation 
and some recommendations for further study. The numerical examples are des- 
cribed in Chapter VIII and the resiilts presented in tabular and in graphical 
form in Appendices C and D. In general, the results give an informative 
pictxire of the lateral response of shallow -bxaried arches to nuclear blast 
loadings under a veiriety of conditions. 

9.1 Discussion 

In interpreting the results of the numerical investigation as pre- 
sented in the appendices, one should keep in mind that the relative values of 
the response are more significant than their absolute values. The relative 
values permit comparison of arch behavior eis geometry and loading parameters 
are varied. The absolute values are somevdiat misleading because of their 
dependence upon the assumptions made in the dynamic analysis. In general, 
the calculated displacements are smaller than the displacements ^ich one 
would expect in an actual structvjre. This is true first because of the re- 
placement of the actual structure having continuous flexibility with a dynamic 
model having only two degrees of freedom, and second because of the idealized 
soil conditions assiamed to exist before and during response. For example, 
any lack of compaction in the backfill surroxmding an actual arch will be 
manifested directly as Increased deformation of the arch because of the 
relatively large influence of the resistance mobilized in the soil on the 
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response of the arch. This is discussed in detail belov. Further, 
displacements of an actual stinicture will be affected by differences be- 
tveen the actual blast vave characteristics and those of the idealized 
blast wave used in these calcxilations. This effect could lead to either 
greater or smaller values of displacement. 

Another factor influencing the numerical results is the absence 
of any damping factor in the equations of motion. While damping is not con- 
sidered to have a significant effect upon the magnitude of the lateral re- 
sponse during the initial portion of the response history, it will reduce 
later maxima in all cases short of failure. For this reason, the comparisons 
made in Appendix C are based on "initial", or first, maximum values at each 
joint. 

A comparison of the natural periods of vibration of the structm*es 
is presented in Tables C.l and C.2. In general, as the depth of cover is 
increased, the natural period of vibration of the structure is decreased. 

This indicates that the "stiffness" of the soil tends to be more significant 
than the additional mass \4iich the soil adds to the vibrating system. At 
first glance the natural periods corresponding to zero depth of cover at the 
crown appear to be out of sequence. This is because the change in mass from 
the "no cover" to the "d^ = O" condition is much greater than the changes in 
mass between any two of the succeeding cover geometries listed. In fact, 
because of the limitations imposed on the dimensions of the soil masses by 
Equations ( 5 . 1 . 3 ) and (5.1.4), the change in mass with increasing depth of 
cover eventually becomes zero. On the other hand, the stiffness of the 
stinicture continues to increase with depth becoming, relatively, quite large 
for the case of a soil depth at the crown equal to one-half of the arch span. 
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It should he pointed out here that because of the natxire of the 
assumptions upon vhich the dynamic model is based, the natural periods of 
vibration given in Tables C.l and C.2 for the "Asymmetric" modes are closer 
to the periods of the actual system than the i>eriods for the "Symmetric" 
modes. The symmetric modes are influenced to a significant extent by the 
extensioneil as veil as the flexural characteristics of the arch vhile the 
asymmetric modes are almost entirely flexural, or "bending", modes. 

Table C .2 contains the natursd periods of the systems as they 
approach their first maximiim response at the leeward haunch. Aside from a 
relatively small initial inward movement of the windward haunch and upward 
movement of the crown in certain cases, this is the first maximvim reached 
by the system (cf. Appendix D). In all of the cases studied the system 
motion is "Symmetric -Out" (cf. Figure 5*5) as the system approaches this 
maximum. Since the natioral period is a function, in part, of the soil stiff- 
ness (cf. Equations (7.c;.9) and ( 7 . 2 .IO)), and the soil stiffness in turn is 
a function, in part, of the overpressure restraining the soil failure wedges 
(cf. Equations (4.4.10) and (4.4.12)), it is necessary to specify the par- 
ticular time at \diich the natxiral period is to be calculated so that the over- 
pressure acting on the soil failure wedges may be found. The time selected 
in each case is the time corresponding to maximum restraint of the "passive" 
soil failure wedges as the system approaches its first maximvan response at 
the leeward hamich. 

In Table C .2 it is interesting to note that, as in Table C.l, the 
natural period tends to decrease with increasing depth of cover. Also, the 
natural period decreases with Increasing overpressure. This results from 
the fact that the overpressure acts to increase the soil stiffness but has 
no effect upon the mass of the system. 
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In Tables C.^, C.h, and C.5 the first maximum values of response 
are compared. The tendency of the arch to move initially inward at the 
windward haianch Eind upweurd at the crown in certain cases (cf. Appendix D) has 
not been considered in selecting "first" maxima. As can be seen from the 
response ciurves in Appendix D, the sequence of response is, in general: 

1. The response reaches a maximum (outward) veilue at the leeward 

haunch. 

2. The response reaches a maximum (downward) value at the crown. 

3. Bie response reaches a maximum (outward) value at the windward 

haimch. 

In all cases as the depth of earth cover is increased, the magnitude 
of response is decreased. In general, as the overpressure is incresised, the 
magnitude of response is increased. Ifowever, this increase in the magnitude 
of response is not as great as might be expected. This is caused by the 
increase in soil resistance with increasing overpressure level. An exception 
to this occurs in the case of the corrugated metal arch with no soil cover at 
the crown. Here an increase in the overpressture loading from 50 psi to 100 psi 
results in a reduction in the msignitude of the response. This is due to the 
fact that the mass remains constant with a change in overpress\ire -vdiile both 
the loading and resistance parameters vary with overpress\ire level. Increasing 
the overpressure level to 1000 psi overcomes the "mass effect" and the magni- 
tude of the response Increases as expected. 

It should be noted here that comparing first maximvon in those cases 
having no soil cover at the crown is somewhat misleading. Examination of the 
response curves in Appendix D reveals that the second meixlma tend to become 
very large for certain combinations of structure geometry and loading (cf. 
Figures D. 7, D.15 and D.22). For these cases, the overpressure is not suf- 
ficient to hold the structure down and the crown "pops up." 
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Another point of interest may be illustrated using the data 
presented in Figures C.J>, C.k and C.5* Biis concerns the relative contribu- 
tions of the arch and the soil to the lateral deformation resistance of the 
system. Consider, for example, the reinforced concrete arch with a depth of 
soil cover at the crown equal to one-fourth of the axch span and subjected 
to a 100-psi shock wave. Now consider the configuration of the structure 
at the time of its first maximum response at the leeward haunch, i.e. at 
t = 0.0173 seconds. From Table C.5 the horizontal displacement of the leeward 
haunch, h^, at this time is +0.005^ feet. From Figure D.J>, the horizontal 
displacement of the windward haunch, h^^, at this time is -0.0048 feet. 
Proceeding as in the sample calculations of Section 8.4 gives corresponding 
values of the equivalent arch and soil resistances an: 

0^ = -0.87 kips = -77.09 kips 

Qj = +0.91 kips = +93.^7 kips 

This illustrates the fact that the arch contribution to the lateral response 
of the system is relatively small as compared with the soil contribution. 

This also indicates that the larger displacements obtained in the case of the 
corrugated metal arch may be due more to the larger span of the corrugated 
metal arch rather them to the smaller resistance of the steel arch rib eis 
compared with the concrete arch rib. 

The data contained in Table C.6 show the effect of varying the 
ratio, of the horizontal to vejrbical overpressures induced in the soil 

by the shock wave in air (cf. Figures D.I7 through D. 20). In general, as 

is increased, the magnitude of the response is decreased. Exceptions to 
this are the initial inward movement of the windward haunch and the accompany- 
ing upward movement of the crown. These movements are magnified by increasing 

K . 

pr 
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The effect of the angle of inclination of the shock front in the 
soil was studied and the results of the study presented in Table C.7 (cf. 
Figures D.3> D. 4 and D. 5). Since the angle of Inclination of the shock 
front in soil, t), is a function of the relative values of the velocity of 
the shock wave in air, U, and the seismic velocity of the soil, C , (cf. 
Appendix B), varying n for a constant U has the same effect as a correspond- 
ing variation in C (cf. Equations (B.3))* Table C.7 lists the maxima for 

s 

the solutions in chronological order. In general, increasing the angle of 
inclination of the shock front in soil with respect to the horizontal 
accentuates the initial inward movement of the windward haunch (and accompany- 
ing upward movement of the crown) but has no significant effect on the magni- 
tudes of the other maxima. 

Response curves for all of the cases studied are presented in 
Appendix D, Figures D. 1 throu^ D. 25. The figures are arranged with those 
for the reinforced concrete arch first followed by those for the corrugated 
metal arch. The figures are grouped by overpressure level and are ordered 
within each group by increasing depth of cover. The solutions for p =50 
psi and p^^ = 100 psi are plotted for the first 90 milliseconds after the 
shock wave reaches the structure, those for p = 1000 psi for I8 milliseconds. 
The cu2*ves include the first maximum response at each joint in all cases. 

Configuration diagrams. Figures D. 24 through D.50, are presented 
for two of the solutions. These diagrams are arranged chronologically and 
illustrate the configuration of the structure at each of the first maxima 
(neglecting the small initial movements of the windward haunch and the crown 
as before). In addition, for purposes of comparison. Figure D. 5O shows the 
configuration of the coi*rugated metal arch at the time of the second maxim^^m 
response at the windward haunch. 
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The configxiration diagrams are drawn first neglecting and then 
including rib shortening. The effect of rib shortening is seen to be a 
genereil inward and downward movement of the structure. While rib shortening 
might seem to cause the arch to move away from the soil and thus invalidate 
some of the assumptions up>on which the lateral response analysis is based, 

(17) 

test resiilts have indicated that the soil will deform also and that the 
defomiations of the soil are of the same order as those of the arch. Further, 
the blast loading is transmitted to the arch by the soil so that any movement 

of the arch away from the soil tends to decresise the load on the arch, thus 

counteracting the tendency for the arch and the soil to move apart. Finally, 

the bleist wave in air tends to push the soil down around the arch by the 

very nature of its action. 

9.2 Conclusions 

The numerical results discussed in the preceding section lead to 
the general concliision that well compacted soil surrounding a shallow-buried, 
undergrovmd arch contributes the major portion of the resistance of the arch 
to lateral deformation when subjected to the blast load res\ilting from a 
nuclear explosion. 

If consideration is limited to those cases where the depth of 
cover is sufficient to prevent the crown from moving upward out of the ground, 
e.g. d^ > b/ 8 , and the arch rib cross section is sufficient to withstand the 
thrust developed in the arch without general yielding or buckling, it may be 
concluded that the arch will successfully withstand the applied loading. 
Further, if the following parameters are varied individually in turn, it may 
be concluded that: 

1 . As the depth of cover is increased, the magnitude of the 



response will decrease. 
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2. As the overpressure level is increased, the magnitude of the 
response will increase. 

3 . As the ratio of horizontal to vertical overpressure induced 
in the soil hy the shock wave in air is increased, the magnitudes of the 
initial movements of the windward haunch inward and the crown upward will 
increase; however, the later magnitudes of the lateral response will decreeise. 

4. As the angle of inclination with respect to the horizontal of 
the shock front in soil is increased, the magnitudes of the initial move- 
ments of the windward havinch inward and the crown upward will increase; 
however, the angle of inclination of the shock front in soil has no signifi- 
cant effect upon the magnitude of the later response. 

5 . As the span of the arch is increased, the magnitude of the 
response will increase. 

6. The maximxnn response of the arch will occ^l^ at the crown, 
principally in the vertical direction. 

In general, it may be concluded that the natural i>eriods of 
vibration of the bending inodes of the arch will be decreased by: 

1. Increasing the depth of soil cover, and 

2. Increasing the oveirpressure level. 

9.3 Recommendations 

While the present study enables one to make certain predictions 
regarding the relative behavior of shallow -buried arches under various 
combinations of geometry and nuclear blast loading conditions, it is limited 
by the assumptions ■vrtilch are made in writing and solving the equations of 
motion for the system. This study might logically be extended in one or 
more of the following areas. 
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There appears to be a significant lack of knowledge regarding the 
interaction between the structure and the soil. While classical eaarth pres- 
sure theory gives an indication of the resistance offered by the soil^ this 
aispect of the problem could well be investigated in detail. For example, 
there is a notable lack of experimental data relating to the amount of de- 
formation required to develop the passive state of plastic equilibrium. 

A soil resistance-deformation relationship somewhat more in agree- 
ment with existing experimenteil data might be obtained by considering a more 
rapid shift in the direction of the deformation resistance offered by the 
soil as the system "rebounds." This approach is unattractive from a computa- 
tional standpoint but may bear further investigation. 

A logical extension of this study is a consideration of more than 
two degrees of freedom for the dynamic system. This would give a closer 
approximation to the actual structure and would permit simultaneous considera- 
tion of extensional and bending modes of response. It rai^t also permit 
footing motions to be considered. 

The determination of the overpressure causing elastic buckling of 
the arch rib is an interesting problem. While the calculations in Chapter 
VIII indicate that this condition is not critical, they are approximate. A 
closer estimate of the critlcaJL overpressinre would be helpful. 

A closer approx itaat ion to actual response values might be obtained 
by considering the rise time of the pressure pulse in the soil. This is dif- 
ficult to define at depths of interest in this problem using existing 
information. In fact, the entire problem of the transmission of the shock 
wave throu^ the soil merits careful study. 

A closer approximation to the later maxima of the displacement-time 
relationships may be obtained by including a damping term in the equations 
of motion. Since this terra would appear to be most significant in those 
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cases with little or no soil cover, it shoxxld be kept in mind that depths 
of soil cover approaching zero at the crown may be unattractive from a 
practical standpoint. This is due to the fact that blast resistance is 
only part of the over-all problem of prxjviding protection against the effects 
of nuclear weapons. Thermal and nuclear radiation effects will sdso be 
significant at overpressure levels of interest in this problem. 
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APPENDIX A 

HYDROSTATIC PRESSURE RATIOS DJ SOIL 
SUBJECTED TO NUCLEAR BLAST LOADCTG 

This study is an investigation of the effect of blast overpressure 
on calculated values of the hydrostatic pressure ratios in soil corresponding 
to the "Active State of Plastic Equilibrivun" and the "Passive State of 
Plastic Equilibrium." The calculations are based upon cleissiceil earth pres- 
sure theory using the Coulomb assumptions to describe conditions in the soil 
failure wedges (cf. Section 2.2). A perfectly homogeneoxos, cohesionless soil 
is assumed with the blast overpressure considered as a sxorcharge loading. 

Two solutions for eeuih of the hydrostatic pressxire ratios are 
compared, one obtained using the "Trial Wedge" method and the other obtained 
using an algebraic expression of Coulomb's method. The example chosen for 
these calculations is a semicircular arch, l6' 8" in diameter, with a depth 
of cover over the crown eq\ial to one -fourth of the arch span. This structure 
is illustrated in Figure A. 1. The blsist loading is taken as the 100-psi 
shock wave resulting from a one -megaton nuclear explosion at the ground 
surface. The pressure surfaces are defined based on the dynamic model as 
described in Chapter IV. The soil is assumed to have a \init weight of 120 
pounds per cubic foot and an angle of internal friction of 30 degrees. The 
effective friction angle between the arch and the soil is taken as two-thirds 
of the angle of internal friction of the soil, or 20 degrees. 

The geometry and soil conditions of this problem are the same as 
those of the illustrative example discussed in Section 8.U. Thus, the 
average overpressvire acting on the failure wedges may be obtained from 
Table 8.6. For the purposes of this study, the maximum value of the average 
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overpressxire restraining the "passive" soil failure wedge is used, i.e. 

p' . = -96. h psi. While the average overpressvire restraining the "active" 
s j 

soil failure wedge has a somewhat hi^er maximum veilue, i.e. p , = -98*5 Psi, 

s j 

it is considered at a later time \dien it has decayed to a value of 96.4 psi 
so that both sets of calculations may be made for the same restraining 
pressure. This is sli^tly more convenient. Actually, since this investi- 
gation deals principally with the effect of "large" surcharge loadings on 
calculated veilues of the hydrostatic pressure ratios corresponding to plastic 
equilibrium in soil, any reasonable value of the restraining pressure may be 
used. 



Coulomb's formula for the hydrostatic pressure ratio corresponding 
to the "Active State of Plastic Equilibrium," K^, as derived in reference (l 4 ), 
is given below as Equation (A. l). 



where 



K. = 



cos^C^-to) 



(A.l) 



2 

cos (Jcos(6+u) 



1 +\l 



^in(5+P)sin(T-l) 1 ^ 



'P = angle of internal friction of the soil 

w = angle of inclination with respect to the verticeil of the 
lower bars of the model 

6 = effective friction angle between the arch and the soil 
i = slope of the groxmd surface with respect to the horizontal. 



Substituting the soil properties listed above and the geometry shown in 
Figure A.l into Equation (A.l) gives = O.5IO. 

^ ^ cos^(50°-22°29.6') 

cos2(22“29.6')cos(20W29.6') [1 . I sto(20°.?0°)Bln ?0° f 

L \1 cos(20 +22 29.6' )cos(22 29.6' )J 



= 0.510 



\ 
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Coulomb’s formula for the hydrostatic pressure ratio corresponding 
to the "Passive State of Plastic Equilibrium," Kp, as derived in reference 
(l4), is given below as Equation (A. 2). 




cos 



Wcos(w-6 




(A. 2) 



Substituting values into Equation (A. 2) gives Kp = 5-^33* 






cos^(30°+22°29.6' ) 



cos^(22°29.6' )cos(22°29.6'-20°) 






.2 

cos(22°29.6'-20°)cos(22°29.6' )- 



sin(30°+20°)sin 30° 1 

3°?Q.6' U 



= 3.‘^33 



The results of the graphical solutions for the resiiltant forces on 
the pressure surface using the "Trial Wedge" method are given in Table A. 1. 
The corresponding hydrostatic pressure ratios are calculated below. 

No overpressure acting ; 



^A = r 



2 ^s*^o 



where P^ = resultant force acting on the pressure surface for an active 
state of stress in the soil 



4,800 

I (120) (12. 500)^ 



• _ 

A " 9>375 



= 0.512 
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s o 



where Pp = resultant force acting on the pressure surface for a 
"passive" state of stress in the soil. 

_ 51.050 

9.575 



= 5.512 



With overpressure acting ; 



^A = r 



^ 7 + p .d 

2 so sj o 



^ 92,500 

9, 575+96 . 4 ( 12. 500 ) ( lh4 ) 

= 0.506 






1 P — 

— 7 d + p' ,d 

2 so •^sj o 



619.000 

182,900 

= 5.584 



Hie inclination of the failure surfaces with respect to the 
horizontal are presented in Table A. 2. Hie Inclinations obtained from the 
"Trial Wedge" solutions are compared with the assumed inclinations of + 
for the "active" failure svirface and for the "passive" failure 

surface for the numerical work in this thesis. 

From the comparative results presented in Table A.l, it may be 
concluded that the overpressure loadings have a negligible effect upon the 
hydrostatic pressure ratios in the soil. Further, both the "Trial Wedge" 
method and Equations (A.l) and (A. 2) give essentially the same values of 
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and Finally, in Table A. 2 it may be seen that the assumed inclina- 
tions of the failure siirfaces in the soil are in reasonable agreement with 
the inclinations obtained from the "Trial Wedge" solutions. 
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TABLE A.l 

HYDROSTATIC PRESSURE RATIOS IN SOIL AS COMPUTED 
BY "TRIAL WEDGE" METHOD AND COULOMB FORMULAS 





State of Stress 
in Soil 


Trial Wedge 


Formula 




Without Blast 


With Blast 


Resultant Force on 


Active, P. 


U.80 


92.5 




Pressvu*e Surface, 


kips 


Passive, Pp 


31.05 


619 




H/dro static 


Active, K. 


0.512 


0.506 


0.510 


Pressure Ratio 




Peissive, Kp 


3.512 


5.58U 


3.‘^35 



TABLE A . 2 

INCLINATION WITH RESPECT TO HORIZONTAL OF 
FAILURE SURFACES IN SOIL 



State of Stress 
in Soil 



Trial Wedge Solution Assumed Value 

Without Blast With Blast 



Active 



61 .° 61 .° <5 + 1) - ^0° 
26° 26° " 30° 



Passive 



1 



12 500 
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Figure A.l. Semi.jircular Arch, E = l6*b”, d ^ = v/k 
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APPENDIX B 

SHAPE OF SHOCK FRONT IN SOIL 



The object of this study is the detennination of the configuration 
of the shock front of the air-induced shock wave in soil. Primary interest 
is limited to those depths below the ground surface associated with the 
shallow-biiTied arches studied in this thesis. Knowledge of the shock front 
configuration is needed in order to define the loading on the arch during 
its envelopment by the shock wave (cf. Chapter VI ). 

The approach used in this study is outlined below. Essentially, 
a i>oint of interest is selected and the location of the shock front of the 
air-induced shock wave in soil determined for the time after the explosion 
corresponding to the time of arrival of the shock wave in air at this point 
of interest. It is assumed that the shock wave in air travels with velocity, 
U, which is a continuously varying f\mction of the peak overpressure and 
ambient atmospheric conditions as expressed by Equation (B.I).^^^^ 



U = c \/l + 6p /tP 
o so o 



(B.l) 



where c = ambient sound velocity in air ahead of the shock front 
° (assumed as II30 feet per second) 

p^^ = peak air overpressxire at the shock front 

P = ambient atmospheric pressure ahead of tlie shock front 
(assimied as l4.7 psi). 

Tlae shock wave in soil is assumed to travel at the seismic velocity of the 

soil, C . This study considers a semi-infinite, homogeneous, soil mass with 
s 

no surface irregularities and with a constant seismic velocity. Further, 
the seismic velocity of the soil is assumed not to vary with depth, direction 
or with the overpressure level in the shock wave. 
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The shock wave tends to travel symmetrically outward frxam Grovind 
Zero as a roughly spherical segment of ever-increasing radius. Therefore, 
this problem may be considered in two dimensions and the shock front assumed 
to extend indefinitely in the third dimension. At distances from Ground 
Zero "vdiich are of interest in this problem, the radius of curvature of the 
shock front in air is quite large, so that this assumption does not lead to 
significant error. 

The calculations in this appendix are based upon a one -megaton, 
contact surface burst. The pressure -distance relationship is taken from 
reference (l7)t The point of interest is selected at a distance from Groxmd 
Zero such that the peak value of the overpressure at the shock front is 
100 psi. Seismic velocities of 1,000, 2,000, ^,000 and 6,000 feet per 
second are considered. 

The procedure followed in the calculations may be divided into 
six steps as follows: 

1. Select the point of interest. 

2. Select a series of check points, 1, 2, between the 

point of interest and Ground Zero. Plot these points to scale along a hori- 
zontal line representing the ground surface. 

3. Calculate the time interval, required for the shock wave 

in air to travel between each of the check points and the point of interest. 

4. Cedculate the distance traveled by the shock front in soil, 
d^^, at the seismic velocity of the soil, C^, during each time interval 
found in Step 3* 

d = Ut) C (B.2) 

sn n s ' ' 

5. Using each check point, n, selected in Step 2 as a center, 

swing a semicirculeir arc in the "soil" with a radius equal to the value of 

d calculated for that check point in Step 4. 
sn 



I8l 



6. A curve tangent to the arcs dravn in Step 5 is thus an 
approximation of the shock front in soil at the instant the shock front 
reaches the point of interest selected in Step 1, 

Ihe results of this study are presented in Figures B.l through 

B.4. 



While this study is quite limited, it may be concluded tentatively 

that for the relatively shallow depths considered in this thesis, the shock 

front in soil may be considered as a plane perpendicular to the direction of 

travel of the shock wave and inclined at an angle, q, with respect to the 

horizontal as given by Equations (B.3)* Tliese equations may be written with 

the aid of Figure B. 5* 

C /U < 1: i| = arcsin C /u 

" " " (B.5) 

Cg/U >1: q = n /2 



where = seismic velocity of the soil 

U = velocity of the shock wave in air at the point of interest. 
While the shock front in soil at depths close to the ground surface for 

C /u > 1 remains essentially vertical, the fact that it tends to "run ahead" 

s 

of the shock front in air may tend to incresise the "rise time" of the 
overpressure -time function in the soil at a particular point of interest. 
Since an increase in the rise time tends to reduce the dynamic response of 
a system to a given impulse, this effect is not considered in the numerical 
work of this thesis (cf. Section 2.4). 
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Figure B.2. Shock Front in Soil: P = 100 psi, C = 2,000 fps 
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Figure B.5* Shock Front in Soil: p = 100 psi, C = ^,000 fps 



18U 



<D 

CD 



O CD - - 




Figure B.4. Shock Front in Soil: p =100 psi, C = 6,000 fps 
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Figure B.5* Inclination of Shock Front in Soil 
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APPENDIX C 

TABULAR COMPARISONS OF NUMERICAL RESULTS 

This appendix contains a series of tables in which comparisons 
are made of the various numerical solutions. In general, quantities are 
expressed as follows: 

Displacements, h,, in feet 
J 

Overpressvures, p, in psi 

Cover depth at crown, d^, in proportion of arch span, B 

Time, t, in seconds 

Natural periods of vibration, T^, in seconds. 

The sign convention used in these tables is the same as that stated 
in Section 1.2, i.e. displacements to the ri^t and upward are considered 
positive. The shock wave is considered to engulf the structure from left 
to right. 

A summary of the mmierical examples considered is contained in 



Table 8 . 5 . 
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NATURAL PERIODS OF VIBRATION WITH NO OVERPRESSURE 
ACTING ON SOIL FAILURE WEDGES 

(seconds ) 



Soil Cover 


Reinforced 


Concrete 


Arch 


Cornigated Metal 


Arch 


Asymmetric 


Symmetric 
In Out 


Asymmetric 


Symmetric 
In Out 


None 


0.147 


0.115 


0.115 


0.484 


0.589 


0.589 


d = 0 
c 


0.165 


0.156 


0.155 


0.175 


0.254 


0.189 


= B/8 


0.144 


0.179 


0.155 


0.185 


0.255 


0.214 


= B /4 


0.129 


0.172 


0.148 


0.162 


0.255 


0.199 


= B/2 


0.092 


0.156 


0.125 


0.111 


0.175 


0.154 




I 

( 

« 
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TABLE C.2 

NATURAL PERIODS OF VIBRATION AT FIRST MAXIMUM 
RESPONSE OF LEEWARD HAUNCH 

( seconds ) 

(Motion: "Symmetric -Out” in all cases) 



Depth of 
Cover Over 
Crovm 


Reinforced Concrete Arch 


Corrugated Metal Arch 


100 psi 


1000 psi 


50 psi 


100 psi 


1000 psi 


0 


0.054 


0.012 


0.060 


0.045 


0.015 


b/8 


0.058 


0.015 


0.075 


0.045 


(1) 


b/u 


0.057 


0.015 


0.075 


0.054 


0.019 


b/2 


0.052 


0.011 


0.062 


0.046 


(1) 



(l) Not calculated. 
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TABLE C .5 

FIRST MAXIMUM RESPONSE AT WINEWARD HAUNCH 
(Neglect small initieil inward movement) 



Cover Para- . 
Depth meter' 
at 

Crown 


Reinforced Concrete Arch 
(Span = 16 '8") 


Corrugated Metal Arch 
(Span = 25' 1 1/4”) 


100 psi 


1000 psi 


50 psi 


100 psi 1000 psi 


0 




-0.0155 


-0.0157 


-0.0191 


-0.0158 -0.0199 




t 


0.0214 


0.0068 


0.0585 


0.0254 0.0087 


b/8 


^1 


-0.0087 


-0.0092 


-0.0121 


-0.0152 (2) 




t 


0.0256 


0.0075 


0.0450 


0.0556 


b /4 


"^1 


-0.0065 


-0.0068 


- 0.0084 


-0.0095 -0.0101 




t 


0.0229 


0.0074 


0.0450 


0.0557 0.0110 


b/2 


^1 


- 0.0040 


-0.0045 


-0.0051 


-0.0057 (2) 




t 


0.0204 


0.0065 


0.0593 


0.0296 



(1) = horizontal displacement of windward haunch, in feet, 
t «= time corresponding to h^, in seconds. 

(2) Not computed. 
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TABLE C .4 

FIRST MAXIMUM RESPONSE AT CROWN 
(Neglect small initial upward movement) 



Cover 

Depth 

at 

Crown 


Para- 

meter^ 


Reinforced Concrete Arch 
(Span = 16 '8”) 


Corrugated Ifetal Arch 
(Span = 25' 1 1/4") 


100 psi 1000 psi 


50 psi 


100 psi 


1000 psi 


0 




-0.0257 -0.0281 


-0.0544 


-0.0275 


- 0.0441 




t 


0.0186 0.0055 


0.0525 


0.0215 


0.0069 


b/8 




-0.0151 -0.0175 


-0.0205 


-0.0218 


(2) 




t 


0.0205 0.0060 


0.0590 


0.0285 




B /4 




-0.0108 -0.0125 


- 0.0141 


-0.0151 


-0.0171 




t 


0.0198 0.0060 


0.0588 


0.0285 


0.0085 


b/2 




-0.0069 -0.0075 


-0.0085 


-0.0095 


( 2 ) 




t 


0.0174 0.0055 


0.0559 


0.0252 





(1) = vertical displacement of crown, in feet, 
t = time corresponding to in seconds 

(2) Not computed. 
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TABLE C .5 

FIRST MAXIMUM RESPONSE AT LEEWARD HAUNCH 



1 Cover 
Depth 
at 

Crown 


Para- 
meter^^ ^ 


Reinforced Concrete Arch 
(Span = 16 *8") 


Corrugated Metal Arch 
(Span = 25' 1 1/4’’) 






100 psi 


1000 psi 


50 psi 100 psi 1000 psi 


0 




+0.0119 


+0.0163 


+0.0181 +0.0162 +0.0268 




t 


0.0165 


0.0050 


0.0293 0.0195 


0.0065 


b/8 




+0.0075 


+0.0094 


+0.0102 +0.0108 


(2) 




t 


0.0181 


0.0055 


0.0341 0.0247 




B/k 




+0.0054 


+0.0066 


+0.0070 +0.0076 +0.0092 




t 


0.0173 


0.0053 


0.0337 0.0246 


0.0075 


b/2 




+0.0037 


+0.0039 


+ 0.0044 +0.0049 


( 2 ) 




t 


0.0154 


0.0045 


0.0295 0.0220 





(1) = hori 2 X»nteLL displeuiement of leeward haunch^ in feet, 
t = time corresponding to h^, in seconds. 

(2) Not computed. 
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TABLE C.6 



EFFECT OF HORIZONTAL-TO -VERTICAL O' 
RATIO ON LATERAL RESPONSE 



(First maximum values) 



'RESSURE 



(2) 

Parameter 











l/b 


1/2 


5/4 


1 


(+hj 
1 max 


0.0007 


0.0015 


0.0022 


0.0054 


Time 


0.0068 


0.0086 


0.0104 


0.0121 


(+hj 
3 max 


0.0076 


0.0057 


0.0045 


0.0054 


Time 


0.0246 


0.0221 


0.0192 


0.0164 


(-bJ 

5 max 


-0.0151 


-0.0106 


-0.0060 


-0.0016 


Time 


0.0285 


0.0279 


0.0268 


0.0258 


(-h^ ) 

1 max 


-0.0095 


-0.0076 


-0.0059 


-0.0044 


Time 


0.0557 


0.0540 


0.0542 


0.0542 



(1) Based on Corrugated Metal Arch, = 100 psi, d^ = fi/h. 

(2) = horizontal movement of windward haunch , in feet 
h^ = horizontal movement of leeward haunch , in feet 
h^ = vertical movement of crown, in feet 

Time, in seconds 
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TABLE C .7 

EFFECT OF INCLINATION OF SHOCK FRONT IN SOIL 
ON LATERAL RESPONSE^ 

(2) 

(First maximum values' ') 



Inclination of Shock 
Front in Soil 

Seismic Velocity in Soil 


90°00' 
5000 fps 


54 ° 22 ' 
2400 fps 


18 ° 44 ' 
1000 fps 


(+hj 

5 max 


+0.0001 


+0.0001 


( 5 ) 


Time 


0,0051 


0.0022 




(+h^ ) 

1 max 


+ 0.0004 


+0.0002 


( 4 ) 


Time 


0.0045 


0,0056 




(+hJ 
5 max 


+0.0054 


+0.0052 


+0.0061 


Time 


0.0175 


0.0165 


0.0157 


(-hj 

5 max 


-0.0108 


-0.0105 


-0.0112 


Time 


0.0198 


0.0189 


0.0174 


(-h,) 

1 max 


-0.0065 


-0.0062 


-0.0058 


Time 


0.0229 


0.0220 


0.0210 


(1) Based on Reinforced 


Concrete Arch, 


d = b/ 4 , p = 

C ' ^ so 


100 psi 



(2) = horizontal movement of windward haunch, in feet 

hj = horizontal movement of leeward haunch, in feet 
h^ = vertical movement of crown, in feet. 

(5) No initial upward movement of crown. 

( 4 ) No initial inward movement of windward haunch. 



APPEITDIX D 



GRAPHICAL PRESENTATION OF NUMERICAL RESULTS 

The results of the numerical vork done in connection vith this 

thesis are presented in this appendix in graphical form. Displacement-time 

"response" curves are included for each of the problems solved. Three 

curves are plotted for each problem as follows: 

hj^ = horizontal displacement of the windward haunch 

h^ = horizonteil displacement of the leeward haunch 

he = vertical displacement of the crown 
I? 

The sign convention used in this appendix is the same as that stated in 
Section 1.2, i.e. displacements to the right and upward are considered 
positive. In addition to the response curves, the configuration of the 
dynamic model is illustrated for two of the problems at several significant 
times in their response histories. 

The problems are described in Chapter VIII and summarized in 
Figure 8 . 5 . A set of sample calculations for one of the problems is pre- 
sented in Section 8 . 4 . The significance of the resvilts is discussed in 



Chapter IX. 
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Figure D.jJ. Response Curves for Reinforced Concrete Arch: = 100 psi, d^ = b/U, Cg = 5>000 fps 









■i 



+ 0 030 



198 




y ‘ lU8LUaOD|dS|Q 



Figure Response Curves for Reinforced Concrete Arch: p = 100 psi, = B/4, C = 2,400 fps 
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Figure D.6. Response Curves for Reinforced Concrete Arch: p = 100 psi, ~ 
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Figure D.7. Response Curves for Reinforced Concrete Arch: p = 1,000 psi 
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Figure D. 9 . Response Curves for Reinforced Concrete Arch: p = 1,000 psi, d = b/4 
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Figure D.IO. Response Curves for Reinforced Concrete Arch: p = 1,000 psi, d = b/2 
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Figure D.ll. Response Curves for Corrugated Metal Arch: p = 50 pci, d = 0 
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Figure D. 15. Response Curves for Corrugated Metal Arch: p = 50 psi> d = b/4 
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Figure D.15* Response Curves for Corrugated Metal Arch: p = 100 psi; d = 0 
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Figure D.I 7 . Response Curves for Corrugated Metal Arch: p = 100 psi, d = b/U, K = l/k 
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Figure D. 19. Response Curves for Corrugated Metal Arch: p = 100 psi, d = b/4, K =3/^ 
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Figure D. 21. Response Curves for Corrugated Metal Arch: p = 100 psi, d = b/2 
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Figure D.22. Response Curves for Corrugated Metal Arch: = 1,000 psi 
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Figlire D.25* Response Curves for Corrugated MetaJL Arch: p = 1,000 psi, 
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Figure D.29. Configuration Diagram for Cornigated Metal Arch: 
p = 50 PBI, d = b/ 4, t = 0.0450 sec., 
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